Next Article in Journal
Numerical Solution of Maxwell Equations for S-Wave Superconductors
Next Article in Special Issue
Crystal Growth Techniques for Layered Superconductors
Previous Article in Journal
Correlated Disorder in Myelinated Axons Orientational Geometry and Structure
Previous Article in Special Issue
Superconductivity in 122-Type Pnictides without Iron
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

In Search of Unambiguous Evidence of the Fulde–Ferrell–Larkin–Ovchinnikov State in Quasi-Low Dimensional Superconductors

by
Mihail D. Croitoru
1,* and
Alexandre I. Buzdin
1,2,*
1
Université Bordeaux, CNRS, LOMA, UMR 5798, F-33405 Talence, France
2
Institut Universitaire de France, F-75231 Paris, France
*
Authors to whom correspondence should be addressed.
Condens. Matter 2017, 2(3), 30; https://doi.org/10.3390/condmat2030030
Submission received: 9 August 2017 / Revised: 1 September 2017 / Accepted: 4 September 2017 / Published: 11 September 2017
(This article belongs to the Special Issue Layered Superconductors)

Abstract

:
In layered conductors with a sufficiently weak interlayer coupling in-plane magnetic field cause only small diamagnetic currents and the orbital depairing is strongly suppressed. Therefore, the Zeeman effect predominantly governs the spin-singlet superconductivity making the formation of the spatially modulated Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase possible in such materials. Despite decades of strenuous effort, this state still remains a profound mystery. In the last several years, however, there have been observed several hints indicating the experimental realization of the FFLO state in organic layered superconductors. The emergence of the FFLO phase has been demonstrated mainly based on thermodynamic quantities or microscopically with spin polarization distribution that exhibit anomalies within the superconducting state in the presence of the in-plane magnetic field. However, the direct observation of superconducting order parameter modulation is so far missing. Recently, there have been proposed theoretically several hallmark signatures for FFLO phase, which are a direct consequence of its main feature, the spatial modulation of the order parameter, and hence can provide incontrovertible evidence of FFLO. In this article, a review of these signatures and the underlying theoretical framework is given with the purpose to summarize the results obtained so far, omitting duplications, and to emphasize the ideas and physics behind them.

1. Introduction

High magnetic fields destroy superconductivity and restore the normal conducting state [1]. The underlying physics are based on two effects. The first one is orbital in nature whereby vortices penetrate into a superconductor and the suppression of superconductivity by a magnetic field H is caused by an increase of diamagnetic energy beyond the critical value H c 2 o r b = Φ 0 / 2 π ξ 2 . The second one is the Zeeman effect, which breaks apart the paired electrons in a spin-singlet state if the field is larger than the Pauli paramagnetic limit, μ B H P = Δ 0 / 2 , where Δ 0 is the superconducting gap at T = 0 and H = 0 [2,3]. In most type-II superconductors, the orbital effect plays the dominant role in suppression of superconductivity. However, the upper critical field becomes high when the orbital motion is suppressed like in quasi-1D or quasi-2D layered superconductors for magnetic field applied parallel to the most conducting layers. Under this condition, the destruction of spin-singlet state of Cooper pairs may become the dominant mechanism of the superconducting state suppression, with the upper critical field determined mostly by the Pauli paramagnetic limit [4].
Due to the paramagnetism of conduction electrons, in a magnetic field, the conduction band is splitted and the normal quasiparticles have separate spin-up and spin-down Fermi surfaces (illustrated in Figure 1). They are displaced due to the Zeeman energy. The Fermi surface mismatch should lead to a pair breaking: in the conventional Bardeen–Cooper–Schrieffer (BCS) pairing regime, the total pair momentum is zero, a pairing state, Cooper pair, can be formed between any opposite parts of the Fermi surface. When magnetic field is above the Pauli limit, it is difficult to form pairs with zero total momentum because the difference in displaced Fermi energies becomes larger than the pairing gap. To overcome this, the Zeeman splitting causes a nonzero momentum of the Cooper pairs, resulting in oscillations of the superconducting order parameter, at the expense of the pairs being formed only between some restricted parts of the Zeeman splitted Fermi surfaces [5]. These parts depend on the anisotropy of the Fermi surface and determine the stabilization of the Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) state (in addition to U(1)-gauge symmetry, the spatial symmetry is spontaneously broken). The existence of inhomogeneous superconducting state was predicted back in the 1960s by Fulde and Ferrell and Larkin and Ovchinnikov [6,7], who pointed out that the stability of the superconducting phase in a magnetic field can be increased beyond the Pauli paramagnetic limit by pairing electrons with different momenta. The FFLO emergence temperature cannot exceed T * 0 . 56 T c 0 . However, the influence of impurities may reduce T * [8]. Physically spatial modulation of the order parameter allows for arranging the polarized quasiparticles in the vicinity of the nodes of the order parameter, and hence compensating for the loss of the superconducting condensation energy by the gain in spin polarization [5]. An FFLO state consists of regions of positive and negative pairing amplitude separated by domain walls where the magnetization is accumulated—a form of phase separation with alternating superconducting and polarized normal regions.
Conditions for the stabilization of the FFLO phase are rather stringent [9], namely (i) the orbital pair breaking effect should be sufficiently weaker than the Pauli paramagnetic limit, the Maki parameter α M 2 H c 2 / H P 1 . 8 ; (ii) the system should be in a clean limit, l ξ 0 = v F / π Δ 0 , where l is the mean-free path of the quasiparticles and ξ 0 is the superconducting coherence length at T = 0 and H = 0 [8,10,11,12,13,14]. Moreover, the highly anisotropic Fermi surface [5,15,16,17] favors the FFLO phase formation. This is an inherent property of layered conductors that exhibit a highly anisotropic structure and hence have features of systems with reduced dimensionality [18].
However, despite intense theoretical and experimental efforts since its prediction in the 1960s, the compelling evidence for nonhomogeneous state in superconductors has been provided only recently. There have been observed several signatures indicating the experimental realization of the FFLO state in organic superconductors for in-plane external magnetic field when the flux penetrates between the layers in the form of Josephson vortices, thus limiting orbital suppression [19,20,21]. In particular, (i) the anomaly in the thermal conductivity [22], and fine details in the phase diagram obtained by the tunnel diode oscillator (rf)-penetration depth measurements and pulsed field techniques [23] for the clean organic sample λ -(BETS) 2 GaCl 4 ; (ii) the magnetic torque evidence for the appearance of an additional first-order phase transition line within the superconducting phase in the in-plane high field regime for organic sample κ -(BEDT-TTF) 2 Cu(NCS) 2 [24,25]; (iii) evidence for phase transition within the superconducting phase obtained from the local electronic spin polarization in 13 C NMR spectroscopy experiment [26], NMR detection of spin-polarized quasiparticles [27] forming the Andreev bound states spatially localized in the nodes of the order parameter [5,28,29], and phase transitions in the H T phase diagram that are consistent with the FFLO phase obtained from rf-penetration depth measurements [30] for organic sample κ -(BEDT-TTF) 2 Cu(NCS) 2 ; (iv) NMR relaxation rate evidence for an additional phase transition line [31] and the clear upturn beyond the Pauli limit in the magnetic-field and angular-dependent high-resolution specific-heat measurements for the organic materials κ -(BEDT-TTF) 2 Cu(NCS) 2 and β -(ET) 2 SF 5 CH 2 CF 2 SO 3 [32,33], (v) magnetic torque evidence for the tricritical point between the FFLO, homogeneous superconducting, and paramagnetic metallic phases in the 2D magnetic-field-induced organic superconductor λ -(BETS) 2 FeCl 4 [34]; (vi) measurements of the temperature and angular dependencies of H c 2 in pnictide superconductor LiFeAs [35] and KFe 2 As 2 [36] as well as a dip in the interlayer resistance in parallel fields observed in the compound λ -(BETS) 2 Fe 1 x Ga x Cl 4 for x = 0 . 37 [37] (for all types of salts, ( x = 1 . 0 , x = 0 . 37 , x = 0 ) signs for the FFLO state have been observed [23,38]) probably suggest that these systems may have realized the FFLO phase. Very recently, the first magnetocaloric and calorimetric observations of a magnetic-field-induced first order phase transition within a superconducting state to the paramagnetic state with higher entropy in κ -(BEDT-TTF) 2 Cu(NCS) 2 molecular compound [39] has been reported. These results provide a strong thermodynamic evidence of the FFLO superconducting phase.
Similarly, quasi-1D superconducting compounds have been studied extensively in order to find out if this class of materials can exhibit the inhomogeneous state. The field-amplitude and field-angle dependence of the superconducting transition temperature T c H of the organic superconductor (TMTSF) 2 ClO 4 in magnetic field applied along the conduction planes have been reported [40,41]. The authors observed an upturn of the curve of the upper critical field at low temperatures (An enhancement of almost two times over H P 27 kOe is observed, H c 2 50 kOe [40]). Moreover, an unusual in-plane anisotropy of H C 2 in the high-field regime was observed. Both observations were interpreted as an evidence of the FFLO state formation [42,43,44].
The search for FFLO has also covered other low-dimensional conductors [45]. In particular, a possible realization of the FFLO phase in ultrathin crystalline Al films has been recently reported [46]. In superconductor-ferromagnet (S/F) bilayers, a quasi-one-dimensional FFLO-like state can be realized by Cooper pairs migrating from the superconductor into the ferromagnet. Due to the exchange splitting in the ferromagnet, the Cooper pair gains a nonzero momentum, resulting in an oscillating pairing wave function [47]. The observation of a π -junction predicted in [48] can be considered as the first direct evidence of the FFLO-like state [49,50].
Besides the quasi-low dimensional organic materials, heavy-fermion systems were shown to be favorable to the formation of the spatially modulated state [51,52,53]. The heavy fermion superconductor, CeCoIn 5 , appears to meet all the strict requirements placed on the existence of the FFLO state. In a strong magnetic field, the superconducting transition is of the first order at low temperatures for both directions of the magnetic field. The first order transition at the upper critical field indicates that the Pauli paramagnetic effect dominates over the orbital effect for both field directions. A second anomaly is the heat capacity that has been observed within the superconducting phase. However, the phase found between the normal metal and the homogeneous superconductor show features expected rather from the antiferromagnetically ordered phase than from FFLO (for example, thermodynamic experiments showed an entropy decrease while shifting from the homogeneous phase to that one [54]). Nevertheless, recent measurements revealed an additional order intertwined with another one, probably an FFLO state [55,56].
Besides condensed matter physics, the FFLO phase has been predicted to occur also in other systems, such as cold atomic gases [57], quantum chromodynamics, nuclear physics, and astrophysics [58,59].
Various theoretical studies of the FFLO phase were performed in order to provide unambiguous ways to detect the manifestation of the FFLO phase. For instance, a tool based on the Josephson effect was suggested, since the application of the magnetic fields to a junction recovers Josephson current that is suppressed between a conventional and the FFLO superconductors [60]. FFLO states in multiband superconductors have been analyzed as well. The analysis shows that the s ± multiband pairing facilitates the FFLO transition and results in the upward curvature of the upper critical field in the intermediate temperature range [61]. The multiband s ± pairing with several “shallow” bands with bottoms close to the FS can result in the FFLO state caused by a small shift of the chemical potential. This may enable one to tune the FFLO transition. In [62], it has been shown that the FFLO state is most stable when the Fermi level is in the proximity of one of the system energy subbands.
In this work, we provide an overview of our recent investigations on possible tools for exhibiting the FFLO phase based on its main feature, namely, the spatial modulations of the order parameter.

2. General Settings

In this section, we consider quasi-2D superconducting materials and model a layered superconductor as a stack of conducting layers (see Figure 2). The single-electron spectrum of such system is taken as follows
ξ p = p x 2 2 m x + p y 2 2 m y + 2 t cos p z d μ ,
where the effective mass approximation is used for description of the in-plane charge carriers motion and the tight-binding approximation for the motion along the z-direction, perpendicular to the conducting layers. The Fermi surface acquires the shape of a corrugated cylinder with an elliptical cross section (see Figure 3). We assume that the corrugation of the Fermi surface, because of the coupling between layers, is small, i.e., t T c 0 , but sufficiently large to make the mean field treatment justified, namely the coupling between layers, t, satisfies the condition T c 0 2 / E F t T c 0 [63,64]. Here, T c 0 is the critical temperature of the system at H = 0 . In this model, an effective mass anisotropy can be taken into account by the scaling transformation of the coordinates and orbital magnetic field [65].
We choose the magnetic field to be parallel to the conducting planes and with a gauge for which the vector potential A = H × r ( r = ( x , y , 0 ) is a coordinate in xy-plane), i.e., A z = x H sin α + y H cos α , where α is the angle between the applied field, with amplitude H, and x-axis. Assuming that the vector potential varies slowly at the interlayer distances and taking into account that the system is near the second-order phase transition, we can employ the linearized Eilenberger equation for a layered superconductor in the presence of the parallel magnetic field. It is the momentum representation with respect to the coordinate z [1,66,67]:
Ω n + 1 2 Π ^ f ω n , r , p z , k z = Δ p ^ ( r , k z ) + f ω n , r , p z , k z 2 τ sign ( ω n ) .
Here, we use the notation
Π ^ v F · + 4 i t sin ( p z d ) sin ( Q · r k z 2 d ) ,
where Q = ( π d H / ϕ 0 ) [ sin α , ( m x / m y ) 1 / 2 cos α , 0 ] with ϕ 0 = π c / e , h = μ B H is the Zeeman energy, v F = v F p ^ is the in-plane Fermi velocity, τ is the impurity scattering time, and Ω n ω n i h + sign ( ω n ) / 2 τ . In this work, we consider both s-wave and d-wave pairing interaction. For the case of d-wave pairing, the interaction may be represented in the form
V d p , p = λ 2 γ p ^ γ ( p ^ ) ,
where γ d x 2 y 2 p ^ = p ^ x 2 p ^ y 2 = cos 2 φ with φ , the azimuthal angle of p (angle between the momentum in the crystalline xy-plane and x-axis). In this expression, factor 2 is introduced for normalization. Then, the order parameter may be written as
Δ p ^ r = Δ d r γ p ^ ,
where Δ d r is defined as self-consistency relation,
1 λ Δ d r , k z = 2 π T ω > 0 γ p ^ f ω p ^ , r , p z , k z ,
where λ is the pairing constant and the brackets denote averaging over p z and p ^ = n ,
π d π d d d p z 2 π 0 2 π d φ 2 π .
We assume that the temperature unit is so chosen that the Boltzmann constant k B = 1 . For the case of an s-wave pairing interaction, γ p ^ = 1 .
The FFLO state appears as a phase with a modulated order parameter with a wave vector q whose direction is determined by the crystal field effects [68,69,70] and the pairing symmetry [71]. The orientation of the FFLO modulation vector is arbitrary in a pure Pauli limited s-wave superconductor in the case of a Fermi surface with an elliptical cross section, which can be mapped by scaling transformation to an isotropic case [65]. For a d-wave superconductor, the symmetry of superconducting pairing fixes the directions of the FFLO modulation vector even in the absence of the crystal field effects [71]. If the orbital effects are essential, the actual direction will be determined by the interplay between the anisotropy of pairing and the crystal field effects. For the triclinic symmetry of the organic superconductors, the FFLO modulation is pinned in a certain direction. The magnitude and the direction of the FFLO modulation vector are determined by the condition of the maximum value of the critical field
H c P 2 = max H c 2 q , ϕ .
The simple exponential solution f ω n p , r , p z exp i qr is no more valid in quasi-2D superconductors. Therefore, to describe correctly the angular dependence of the upper critical field in the FFLO phase, we have to account for the orbital effects, which add the higher harmonics in FFLO modulation, q ± m Q . The solution of Equation (2) might be written as
f ω p ^ , r , p z = e i qr m e i m Qr f m ω n , p ^ , p z .
Because of the form for f ω p ^ , r , p z in Equation (8), one can write the order parameter Δ r as
Δ r = e i qr m e i 2 m Qr Δ 2 m .
From symmetry considerations, it follows that Δ 2 m = Δ 2 m . Substituting Equations (8) and (9) into Equation (2), one gets the following system of coupled equations [72]:
L n q f 0 + t ˜ f 1 t ˜ f 1 = Δ 0 ,
L n q ± Q f ± 1 ± t ˜ f 0 t ˜ f ± 2 = 0 ,
L n q ± 2 Q f ± 2 ± t ˜ f ± 1 t ˜ f ± 3 = Δ ± 2 ,
L n q ± 3 Q f ± 3 ± t ˜ f ± 2 = 0 ,
where L n q = Ω n + i v F q / 2 and t ˜ = t sin p z d . Here, we have taken into account that Δ ± 2 m + 1 = 0 and introduced the notation f m f m ω n , n p , p z . This hierarchy of coupled equations is broken on the level of f ± 3 in order to obtain symmetric equations for the first two harmonics of the order parameter up to the second order with respect to the small parameter t / T c 0 .

3. Angular Dependence of the Upper Critical Field

From Equation (10) we can obtain function H c 2 q , ϕ that should be minimized in order to obtain the absolute value and the direction of the modulation vector. In the Pauli limit, when neglecting the orbital motion, we obtain
ln T c 0 T c P = π T c P n 1 ω n 2 γ 2 p ^ L n q F ( h ˜ c P , q ˜ c P ) ,
where the renormalized variables h ˜ = h / 2 π T and q ˜ = q / 2 π T are introduced,
F ( h ˜ , q ˜ ) π T n 1 ω n 2 g 1 v 2 g 1 v 2 Ω n 2 g 1 v g 1 v + 2 Ω n 2 cos 4 ϕ .
T c P is the temperature of the onset of the superconductivity in the pure Pauli regime, g 1 v q 2 v F 2 + 4 Ω n 2 , and ϕ is the angle the q vector makes from the x-axis. The H c P 2 and the modulation vector, maximizing H c P 2 , are illustrated in Figure 1. In the case of an s-wave superconductor (obtained from Equation (15) by substitution ϕ = π / 8 )
F ( h ˜ , q ˜ ) π T n 1 ω n 2 g 1 v .
For a d-wave superconductor, the paramagnetic upper critical field is never smaller than H c P 2 for an s-wave superconductor [71], as shown in Figure 4. The optimal direction of the modulation vector in a d-wave superconductor is ϕ = ± π / 4 , ± 3 π / 4 for T * * < T < T * and ϕ = 0 , ± π / 2 , π for 0 < T < T * * with T * * 0 . 056 T c 0 . H c P 2 and q in s-wave pairing symmetry is shown in the Figure 4a.
Let us now account for the orbital effects. As it has been seen above [73], in the absence of orbital effect, the FFLO state only appears at T < T * 0 . 56 T c 0 and H > H * 1 . 06 T c 0 / μ B . Therefore, the order of the magnitude of the magnetic field required to observe the FFLO state can be found from the relation μ B H T c 0 . Taking this into account, one obtains v F Q v F e d T c 0 / μ B c ( d / a ) T c 0 , where a is the unit cell in the xy-plane. Therefore, v F Q T c 0 . Due to the assumption t T c 0 v F Q , one has t T c 0 v F Q . This condition allows us to hold only the terms up to the first harmonics in Equations (8) and (9). Then, one simply gets
L n q f 0 + t ˜ f 1 t ˜ f 1 = Δ 0 ,
L n q ± Q f ± 1 ± t ˜ f 0 = 0 .
The solution of this system of equations is
T c P T c T c a A t 2 ,
where the following notations are used
a = π T n , ξ = ± T n q , q , ξ Q T = T c P ,
T n g 1 , g 2 , g 3 = 1 2 0 2 π d φ 2 π 2 cos 2 2 φ L n g 1 L n g 2 L n g 3 ,
A 1 h T c P T c P h = 1 1 h ˜ F ( h ˜ , q ˜ ) h ˜ T = T c P .
Here, the expression for T n g 1 , g 2 , g 3 can acquire an analytical form
T n g 1 , g 2 , g 3 = 8 t 2 k = 1 N k 4 a p k 2 p k + 1 p k + 2 + i g k 2 g k g k + g k g k + 1 g k g k + 1 + × 1 + g k 4 g k g k + 2 g k g k + 2 + 2 a g k v ,
where we have introduced the following notations g k v g k 2 v F 2 + 4 Ω n 2 , g k = g k , x i g k , y , g 1 q , g 2 q ± Q , g 3 q ± 2 Q and
g k ± i 2 Ω n ± g k v g k v F ,
where k is the cycling index with k = 1 , 2 , 3 .
In our numerical investigations, we restrict ourselves to the parameters of the compound κ -(BEDT-TTF) 2 Cu(NCS) 2 . The Maki parameter, α 8 [24], the interlayer coupling is t = 1 . 5 K ( 1 . 8 K) [74], t / T c 0 = 0 . 16 ( 0 . 2 ), Δ 0 = 2 . 8 k T c 0 [75] and the Fermi velocity v F = 5 . 0 10 . 0 × 10 4 m / s [76]. We have chosen the value v F = 7 . 5 × 10 4 m / s [74]. Introducing the dimensionless Fermi velocity parameter, η = v F π d / ϕ 0 μ B , this value of v F corresponds to η = 2 . 55 . The interlayer distance is d = 1 . 62 nm . The summation over the Matsubara frequencies was performed numerically.
In the following for the case of a d-wave superconductor, we consider two situations. In the first one, the symmetry of superconducting pairing fixes the direction of the FFLO modulation wave vector, characterized by ϕ , the angle that q makes from the x-axis: at high temperatures of the FFLO phase, ϕ = ± π / 4 , ± 3 π / 4 , while, in the low-temperature phase ϕ = 0 , ± π / 2 , π thus making it four-fold degenerate. In the second case, the crystal field effect leaves this degeneracy.
Figure 5 presents the orbital-motion induced normalized correction of Δ T c = T c T c P as a function of T c P / T c 0 for different directions of the applied field when the symmetry of superconducting d-wave pairing fixes the direction of the FFLO modulation wave vector. The direction of the external field is measured by α , the angle the applied field H makes from the x-axis. Due to the four-fold symmetry of the d x 2 y 2 pairing, we plot Δ T c T c P / T c 0 curves only for α = 0 , π / 9 , π / 4 . The regions T c P < T * * and T c P > T * are shadowed. In the low-temperature region, the modulation vector q is along the x-axis, while at T c P > T * q = 0 , and it describes Δ T c in the conventional phase. The central domain illustrates Δ T c , when the modulation wave vector q makes the angle ϕ = π / 4 with the x-axis. The re-orientational transition at T c P < T * * is accompanied by a strong change of the orbital correction, which should result in the sharp step on the T c H curve.
The existence of the nodes in the order parameter results in particular features of the anisotropy of the superconducting onset temperature induced by the spatially modulated FFLO phase (in the case of d-wave only, when the symmetry of superconducting pairing fixes the direction of the FFLO modulation wave vector). Figure 6 shows the magnetic field angular dependence of the normalized superconducting transition temperature, T c α / T c P , calculated at T c P / T c 0 0 . 03 , 0 . 056 , when the modulation vector is fixed by the symmetry of pairing to the direction of maxima of the order parameter, and at T c P / T c 0 0 . 057 , 0 . 075 , 0 . 15 and 0 . 4 , when it is fixed to the nodes of Δ p ^ r . In the polar plot, the direction of each point seen from the origin corresponds to the magnetic field direction and the distance from the origin corresponds to the normalized critical temperature. One can see that the field-angle dependence of the onset of superconductivity makes evident two transitions: one at T c P / T c 0 = 0 . 56 , the transition from the conventional phase to the FFLO modulated phase, accompanied by a rotation of the anisotropy principal axis by the angle π / 4 . The second transition occurs at T c P / T c 0 0 . 056 , when the direction of the FFLO modulation vector rotates by ϕ = π / 4 , which is reflected in a change of the overall anisotropy of the onset of superconductivity. The transition from the high temperature FFLO phase to the low temperature phase is accompanied by a strong decrease of the orbital effect that makes T c α / T c P dependence at T < T * * much closer to the pure paramagnetic limit.
In real systems, the crystal field effect is inevitably present and influences the pinning direction of the FFLO modulation [68,69]. As already mentioned previously, the orientation of the FFLO modulation vector is arbitrary in the pure Pauli limited regime for an elliptic Fermi surface (s-wave case). The crystal field introduces deviations from the ellipticity and pins the FFLO modulation vector in a certain direction. In the case of a d-wave superconductor and weak crystal effect, it can lift the four-fold degeneracy of the direction of FFLO modulation, making it two-fold. To illustrate this, we pin the FFLO modulation vector along the ϕ = π / 4 , 5 π / 4 directions and investigate the anisotropy of the superconducting onset in the temperature range T * * < T . Figure 7 displays the orbital-motion induced normalized correction of the transition temperature, Δ T c , as a function of normalized temperature, T c P / T c 0 , for several orientations of the applied field (for both s-wave and d-wave cases). Here, the crystal field effect is sufficient to break the four-fold degeneracy to fix the direction of the FFLO modulation vector. The dashed lines illustrate the result for Δ ± 2 = 0 , considered in this section, while the solid lines are the solutions with Δ ± 2 0 , discussed in the next section.
As it was intuitively expected, the orbital effects reduce the superconducting onset temperature, Δ T c < 0 in both cases. While increasing the applied magnetic field, Δ T c first decreases in most cases until the tricritical point, H * , is reached [77]. At H > H * , the function Δ T c ( H ) strongly depends on the in-plane effective mass anisotropy (not considered here and investigated in [42] and angle α (and this is very different from the conventional phase). For some angles, Δ T c may exhibit an upturn and T c may approach the paramagnetic limit, T c P , when H increases. In contrast, for small α , an increase of the magnetic field leads to a decrease of Δ T c . For intermediate angles, Δ T c can be a non-monotonic function of the field strength. One can infer that the strong field-direction dependence of the superconducting onset temperature, T c α , appears at high magnetic fields when the FFLO state develops, while it is absent at low fields.
The change in the anisotropy of the superconducting onset temperature that is induced by the FFLO phase is particularly visible in Figure 8 and Figure 9, where the same plot is given as in the Figure 6 but for the s-wave and d-wave cases. For magnetic fields below H * , the behavior of the upper critical field is isotropic. When increasing H above H * , a strong in-plane anisotropy of H c 2 develops, which remains and becomes essentially pronounced at high fields. In particular, this behavior develops with external magnetic field for the case of the isotropic in-plane Fermi surface. The maximum transition temperature is for the magnetic field orientation perpendicular to the direction of the FFLO modulation vector along the x-axis. In Figure 9, we see that the shape of the field-angle dependence of the onset of superconductivity in the high temperature FFLO phase is similar to that obtained for an s-wave superconductor [42]. However, the principal axis of the plot is not vertical but tilted by π / 4 and fixed along the direction of the modulation vector. Even in the presence of the “easy axis” along the direction π / 4 induced by the crystal field, the four-fold degeneracy of the direction of the FFLO modulation vector restores again for T < T * * . Therefore, anisotropy in the low temperature FFLO phase is presumably the same as that discussed for Figure 6.
The results discussed in this section show that, in the FFLO phase, the in-plane upper critical field anisotropy in quasi-low-dimensional superconductors is settled by the interplay between the modulation and magnetic field wave vectors. The superconducting onset temperature is maximal for the field oriented perpendicular to the FFLO modulation vector. The change of the anisotropy of the critical field as well as of its fine structure may give important information about the FFLO state and unambiguously prove its existence. Our calculations support the interpretation of the experimentally observed in-plane anisotropy of the onset of superconductivity in (TMTSF) 2 ClO 4 samples as a realization of the FFLO state with the modulation vector close to the b * -axis [40].

4. Resonance In-Plane Magnetic Field Effect

We know that the vector potential of the parallel magnetic field results in a modulation of the interlayer coupling with the period λ H = 2 ϕ 0 / d H , where d is the interlayer distance and ϕ 0 = π c / e . In this section, we show that the period of this modulation may interfere with the in-plane FFLO modulation leading to the anomalies in the critical field behavior. The strongest effect occurs when the magnetic length λ H coincides, i.e., in the resonance with the period of FFLO modulation λ FFLO (for T = 0 , λ FFLO 0 = π v F / Δ 0 = π 2 ξ 0 ). The upper critical field at T = 0 observed in κ -(BEDT-TTF) 2 Cu(NCS) 2 is H c 2 30 T. Since in this compound the interlayer distance d = 1 . 62 nm [78], ξ 0 = 7–9 nm this leads to λ H = 85 nm and λ FFLO 0 = 70–90 nm. In (TMTSF) 2 ClO 4 H c 2 0 = 5 T, d = 1 . 31 nm, ξ 0 = 45 nm, resulting in λ H = 630 nm and λ FFLO 0 = 444 nm. In λ -(BETS) 2 FeO 4 , the lower H c 2 0 = 18 T, d = 1 . 85 nm, ξ 0 = 8 . 5 nm, resulting in λ H = 130 nm and λ FFLO 0 = 85 nm. For these compounds λ FFLO 0 < λ H 0 at T = 0 . The FFLO modulation appears only at T < T * 0 . 56 T c 0 with a wave vector q growing from q = 0 to q 0 = q T = 0 = 2 Δ 0 / v F with decreasing temperature. Therefore, if the condition λ F F L O 0 < λ H 0 is satisfied at T = 0 , then, at some finite temperature T, the resonance condition, λ FFLO T = λ H T , should be realized. It corresponds to the situation when the strongly overlapping Josephson vortices form the rectangular lattice with its centers just above the nodes of the order parameter (previously, such a mechanism of the pinning of Josephson vortices by the nodes of FFLO modulation was suggested in [79] (see Figure 10) and observed in organic superconductor λ -(BETS) 2 FeCl 4 [80]).
From the system of Equations (10) and (13), it is seen that, if L n q = L n q ± 2 Q , then the averaged Equations (10) and (12) for Δ 0 and Δ ± 2 show that Δ ± 2 is of the same order as Δ 0 . To account for such degenerate or resonance situations, Equations (12) and (13) are included in our consideration. Making use of the self-consistency relation, one obtains for T < T * in a second-order approximation on the small parameter t / T c 0
Δ 0 1 λ π T c n 2 γ 2 p ^ L n q p ^ + t 2 a = Δ ± 2 t 2 c ± ,
Δ 2 1 λ π T c n 2 γ 2 p ^ L n q ± 2 Q p ^ + t 2 b ± = Δ 0 t 2 c ± ,
where, additionally to Equation (20), the following are introduced:
b ± = π T n , ξ = ± T n q ± 2 Q , q ± 2 Q , q ± 2 Q + ξ Q T = T c P ,
c ± = π T n T n q , q ± Q , q ± 2 Q T = T c P .
Taking into account the fact that the critical temperature when accounting for the orbital effects, T c , is close to T c P , Equations (25) and (26) can be written as the following system of coupled equations:
Δ 0 T c P T c A T c + t 2 a = Δ ± 2 t 2 c ± ,
Δ ± 2 T c P T c T c B ± + t 2 b ± + δ ± = Δ 0 t 2 c ± ,
where
B ± A h ˜ Λ ± h ˜ h ˜ T = T c P A h ˜ Λ ± h ˜ h ˜ T = T c P ,
Λ ± = π T n 2 γ 2 p ^ L n q p ^ 2 γ 2 p ^ L n q ± 2 Q p ^ ,
δ ± = Λ ± T = T c P .
The solution of the system of Equations (29) and (30) is given as
T c P T c T c a A + b ± B ± t 2 + B ± δ ± 2 + 1 2 a A + b ± B ± t 2 + B ± δ ± 2 4 A B ± a b ± c ± 2 t 4 + a t 2 δ ± 1 2 ,
where, for ± , the values are chosen that maximize the critical temperature. Usually, the second harmonic of the order parameter, Δ ± 2 , can be neglected because the developed theory is valid up to the second order with respect to t / T c 0 , and because t T c 0 v F Q . In this case, the solution of the system is just T c P T c / A T c = t 2 a , as it was shown in the previous section. However, if δ ± = 0 (the resonance condition), the term T c P T c / T c B ± + t 2 b ± + δ ± in the l.h.s. of Equation (30) is of the same order as the corresponding term in Equation (29) (up to the second order with respect to t / T c 0 ). Consequently, Δ ± 2 becomes on the order of Δ 0 , and one has to consider both equations in the system on equal footing [81].
Figure 7 displays the change of the transition temperature, Δ T c = T c T c P , due to the orbital effects of the applied magnetic field as a function of the reduced temperature for different field orientations. The solid curves display the results of the second harmonic correction, and they are the solutions of the full system of these equations. One can see (more pronounced in the s-wave case) that two types of curves (dashed and solid) almost coincide for angles α = 0 , 20 , 45 . However, there are essential differences in the curves’ behavior for the angles α = 70 , 90 . These differences are induced by the just discussed resonance effect. The curve Δ T c P for α = 70 exhibits only one cusp at T c P / T c 0 0 . 2 . T c P / T c 0 0 . 2 δ ξ is close to zero and this vicinity induces a wide pseudo-resonance peak on the curve Δ T c P for α = 70 . Our results show that the resonance contribution can become on the order of the “orbital” effect itself. At exact resonances, the suppression of the critical temperature due to the orbital effect becomes negligible. The same conclusion is valid for the case of a d-wave superconductor.
The resonance effect results in particular features of the anisotropy of the superconducting onset temperature induced by the spatially modulated FFLO phase. The thick curves in Figure 8 and Figure 9 show the magnetic field angular dependence of the normalized superconducting transition temperature, T c α / T c P , when accounting for the second harmonics of the order parameter, Δ ± 2 0 . We see that, in addition to the overall anisotropy induced by the FFLO modulation and studied in [42], cusps develop for certain directions of the applied field, when the resonance conditions are realized. These cusps are the result of the resonant interplay between the FFLO wave vector and the magnetic wave vector, when the orbital effects of the field are taken in the second-order approximation.
In this section, we have shown that, in layered superconductors under the applied in-plane magnetic field in the FFLO phase, the resonance between the modulation wave vector and the vector potential may lead to anomalous cusps in the field-direction dependence of the upper critical field. Therefore, we suggest that observation of characteristic cusps in the anisotropy of the onset of superconductivity may be direct evidence for the appearance of the FFLO phase in layered superconductors.

5. FFLO Lock-In Effect

In this section, we keep studying the effects of orbital contribution on the FFLO modulation, targeting unambiguous proof of the FFLO phase. We demonstrate the locking phenomena—the result of the interplay between two competing length scales: the period of the magnetic field potential and the characteristic length associated with the FFLO modulation [82]. The locking occurs when the modulation wave vector remains equal to the period (or matches m periods) of the potential through some range of the magnetic field values. To study effects of competing periodicities in the system, we use the Frenkel–Kontorova (FK) model [83,84].
In the rest part of the article for the purpose of simplicity, we consider a quasi-1D conductor with the following electron spectrum E p = p x 2 / 2 m x + 2 t y cos p y d y + 2 t z cos p z d z , where d y and d z are the inter-chain distances along the y- and z-axis, respectively. We assume that the couplings between chains are small, i.e., t z < T c 0 and t y < T c 0 , but sufficiently large to stabilize the superconducting long-range order and to make the mean field treatment justified, T c 0 2 / E F t z , T c 0 2 / E F t y [63]. Here, T c 0 is the critical temperature of the system at H = 0 . In quasi-1D superconductors, the orbital effect is extremely weak for the magnetic field applied along the chains. Indeed, near T c 0 , the parallel (along the x-axis) upper critical field is H c 2 x = ϕ 0 2 π ξ y ξ z T c 0 T T c 0 , where ξ y = d y t y T c 0 and ξ z = d z t z T c 0 , while the perpendicular upper critical field (i.e., along the y-axis) is H c 2 y = ϕ 0 2 π ξ x ξ z T c 0 T T c 0 , with ξ x v F T c 0 . We see that H c 2 x H c 2 y and the orbital effect for the parallel component of the magnetic field is weakened by the factor t y T c 0 d y ξ x 1 . Therefore, in our analysis of the orbital effect, it is enough to take into account only the perpendicular component of the magnetic field because it provides the dominant contribution to the orbital effect. For the purpose of clarity, we consider the magnetic field to be aligned in the xy-plane, making the angle θ with the x-axis. The perpendicular component of the magnetic field H y = H sin θ , and the corresponding vector potential may be chosen as A z = x H y = x H sin θ , A y = 0 , A x = 0 . Then, the Eilenberger equations reduces to the following linearized Ginzburg–Landau equation (for details, see the Appendix)
α Δ x = β x 2 Δ x δ x 4 Δ x 48 t 4 ν sin 4 χ A Δ x t 2 2 μ 1 7 cos 2 χ A x χ A 2 Δ x + 2 t 2 μ 3 sin 2 χ A x χ A x Δ x + 3 sin 2 χ A x 2 Δ x 8 t 2 γ sin 2 χ A Δ x ,
where Q = π d z H y / ϕ 0 = π d z H sin θ / ϕ 0 with ϕ 0 = π c / e , t t z , α ln T c T c 0 K 1 T c K 1 0 T c = T c T c P A T c and α < 0 in the uniform superconducting state. This equation contains the magnetic-field induced potential, which is periodic in real space with the period λ H = 2 ϕ 0 / d z H y . Equation (35) leads to the functional obtained in the limit of t < v F Q as [85]
F s n = 1 L x d x α Δ x 2 + β x Δ x 2 + δ x 2 Δ x 2 + 1 2 t 2 μ 1 7 cos 2 Q x Q 2 Δ x 2 + 6 t 2 μ sin 2 Q x x Δ x 2 + 4 t 2 γ 1 cos 2 Q x Δ x 2 .
The coefficients found in this expression are provided in the Appendix [86].
It is known [87] that, without the orbital contribution in 1D, the solution of the Ginzburg–Landau (GL) equation is the Jacobi elliptic sine function Δ x = v F ξ k k sn x ξ k , k , expressed in terms of the modulus k, which is determined upon minimizing the free energy at the fixed external field parameter h. This solution in the vicinity of the phase transition line takes the simple form, Δ 0 cos ( q x ) with q the absolute value of the modulation wave vector. The linearized Ginzburg–Landau functional (36) and Equation (35) in the paramagnetic limit (when neglecting the orbital contribution) reduce to
F s n P M = 1 L x d x α Δ x 2 + β x Δ x 2 + δ x 2 Δ x 2 ,
and α Δ x = β x 2 Δ x δ x 4 Δ x , respectively, that provides the expression for the modulation wave vector in the paramagnetic limit, q 2 = β / 2 δ = 2 K 3 T c P / K 5 T c P v F 2 > 0 . This wave vector increases from 0 at the tricritical point to a very large value at T 0 .
So far, we have assumed that the absolute value of the FFLO wave vector might be determined in the paramagnetic limit only. However, the orbital contribution, in its turn, may influence the absolute value of the FFLO modulation vector. Making use of Equations (29) and (30), we may optimize its solution with respect to the vector q. The result of the calculations are depicted in Figure 11, which illustrates the modulus of the FFLO wave vector versus the reduced temperature. One can see that the orbital effect on the modulation vector is weak except for the region in the close proximity to the resonance, q Q , where it may strongly influence the FFLO structure [88,89]. In the vicinity of the resonance, where the unperturbed q and Q curves intersect, an interesting lock-in effect appears: two wave vectors hybridize.
To investigate in this locking effect in detail, we apply the theory developed by Dzyaloshinskii [90,91,92] and look for the solution of Equation (35) in the following form Δ 0 cos Q x + φ x , where φ ( x ) is a slowly varying function with φ ( x ) Q . This implies that the orbital effect weakly modifies the exact one-dimensional solution, which is justified for T c 0 T / T c 0 t / T c 0 2 . We insert this solution into the functional (36), perform the average over the period Q 1 as = Q 2 π π / Q π / Q d x , and taking into account that φ x Q x , obtain
F s n Q = Δ 0 2 L x d x α 2 + 1 2 β Q + φ 2 + 1 2 δ Q + φ 4 + t 2 γ 2 cos 2 φ x + 1 8 t 2 μ Q 2 2 7 cos 2 φ x + 6 8 t 2 μ 2 cos 2 φ x Q + φ 2 .
In the adopted approximation, the functional F s n depends only on the space varying phase, φ x . Since we are interested in the behavior of the system in the vicinity of the resonance q Q , we can write the functional in a simpler form (normalized by π K 5 T c Δ 0 2 , note that K 5 T c < 0 )
F s n φ Q = Δ 0 2 L x d x α 2 + 1 8 v F 4 Q 2 φ δ q 2 + 5 8 t 2 v F 2 Q 2 11 16 t 2 v F 2 Q 2 cos 2 φ x ,
where δ q Q q is the detuning, or the relative misfit between two periodicities in the system. The phase φ x is the shift of the FFLO modulation relative to the minima in the potential cos 2 φ x . The state φ x = 0 is the commensurate phase between Q and q. Equation (39) is a continuum limit approximation of the Frenkel–Kontorova Hamiltonian introduced by Frank and Van der Merwe [93].
Taking into account the periodicity of the function φ x , we can obtain the ground state minimizing F s n φ Q , which is given by the solutions of the exactly integrable sine-Gordon equation for the phase φ x
2 2 φ x 2 + v sin 2 φ x = 0 ,
where v 11 t 2 / v F 2 is the effective interlayer coupling parameter. In the absence of this coupling, v = 0 , Equation (40) has a solution φ x = C x , which describes an unperturbed one-harmonic incommensurate phase, Δ 0 cos Q + C x . On the other hand, Equation (40) has a trivial solution, φ = 0 , that corresponds to a commensurate structure, Δ 0 cos Q x . For a finite v , Equation (40) should describe some inhomogeneous distribution of the phase of the order parameter. The first integral of this equation is given by
φ x 2 = v κ 2 1 κ 2 sin 2 φ x ,
where κ 2 2 v / ε + v and ε is the constant of integration. Then, the exact solution of Equation (41), expressed in terms of the Jacobi amplitude, is
φ x = am v κ x , κ 2 ,
where κ is the constant of integration and must be found from the energy minimum of the system. The solution is the regularly spaced solitons, a soliton lattice. The soliton lattice is a compromise between term cos 2 φ x , which favors φ x = const and the derivative part that favors φ x = C x . For κ 1 , it reduces to φ x = 2 tan 1 tanh v 2 κ x describing a domain wall, which separates two commensurate regions, one with phase φ = π / 2 , and the other with φ = π / 2 . Within this wall (soliton), a π change of the phase occurs. The length of the soliton can be a small meaning the fast change of the phase.
With the solution for φ x given by Equation (42), the normalized energy of the state is expressed in terms of complete elliptic integrals E an K as
F s n φ = δ q 4 π 2 v κ K κ 2 v 16 1 2 κ 2 + 4 E κ 2 κ 2 K κ 2 .
Minimization of the energy with respect to κ leads to the following equation for κ :
κ E κ 2 = v v c ζ ,
where v c = π 2 δ q 2 / 4 . Here, the parameter κ varies in the range 0 κ 1 as long as v varies in the range 0 v v c . The variation of the parameter κ results in a drastic changes in the behavior of the phase φ x . Figure 12 displays the phase for several values of parameter ζ 0 . 2 , 0 . 6 , 0 . 8 , 0 . 9 , 0 . 99 , 0 . 999999 . As ζ is increased a plateau section appears within the period of soliton lattice L = 2 K κ 2 κ / v or δ q L = ( 2 / π ) 2 E κ 2 K κ 2 . With increasing ζ , not only the period of the soliton lattice increases but also the relative fraction of the plateau width. On this plateau, the phase remains almost constant, but strongly changes beyond plateau at the end of the periods. Thus, the system can be represented as a periodic structure of domains of the locked phase separated by solitons. Equation (43) can be expressed differently via the distance L between the domain walls, namely
F s n φ = π 4 δ q L v 16 1 2 κ 2 + 4 κ 2 E κ 2 K κ 2 .
One can see that, as δ q 0 , parameter ζ increases its value. In the limit ζ 1 , the period of the soliton lattice diverges
lim ζ 1 L = 2 κ v ln 4 1 κ 2 ,
and the modulation wave vector is negligible. For ζ 1 , there are no real valued solutions of Equation (44). Consequently, there exists a limiting value of δ q c 2 11 t / π v F , beyond which the phase is always locked, and the commensurate phase remains stable as δ q decreases further. This means that the modulation vector coincides with the modulation induced by the field. Therefore, at δ q c , we have a lock-in of the FFLO modulation to the field induced potential, i.e., incommensurate/commensurate type phase transition. This corresponds to L and the FFLO modulation is described by Δ 0 cos Q x . Figure 13 illustrates the behavior of the spatially varied order parameter in the vicinity of the incommensurate–commensurate (I–C) phase transition when ζ 1 . One can see the domain walls (phase solitons) in the spatial distribution of the order parameter.
The theory developed is applicable not only to the quasi-1D systems but to the quasi-2D materials as well. Since the orientation of the FFLO modulation vector is determined by the crystal field effects [68,70] and/or the pairing symmetry [71], when the in-plane magnetic field is perpendicular to the FFLO wave vector q we reduce the quasi-2D case to the just described quasi-1D case, Label (36), but with the renormalized coefficients: β π 8 K 3 T v F 2 , δ 3 π 128 K 5 T v F 4 , μ π 4 K 5 T v F 2 , γ π 4 K 3 T . The coefficients γ , ν have the same values.
In this section, we have shown that, besides the anomalous cusps in the temperature and angular dependencies of the in-plane critical field, which can be considered as markers for the FFLO phase, a new marker is defined, which is based on the fact that the proximity to the resonance condition is characterized by the field-induced phase transition from the incommensurate to the commensurate structure (C–I) of the spatial modulation of the order parameter with respect to the magnetic wave vector structure. The periodic potential tends to “lock-in” the FFLO modulation into a commensurate configuration; in the incommensurate phase near the C–I transition, the periodicity of the phase solitons or phase-domain walls depends continuously on the temperature and external field strength.

6. Conclusions

In the past few decades, the quest to make evident the FFLO phase has stimulated a lot of activities in this branch of the field of superconductivity and considerably deepened our understanding of the physics behind. The FFLO concept is more than 50 years old and now we are probably at the beginning of its most interesting period, experimental evidence. The developed theoretical framework, both described or discussed/cited or even not mentioned in this review, are very timely to support the burst of the inhomogeneous phase experimental observation [39].
In this article, we have reviewed our contribution to the research field of inhomogeneous superconductivity and, in particular, we have shown that
-
The FFLO modulation strongly interferes with the orbital effect and provides the main source of the in-plane critical field anisotropy. The change of the anisotropy of the critical field as well as of its fine structure may give important information about the FFLO state and unambiguously prove its existence.
-
As soon as the vector potential of the applied magnetic field is commensurate with the wave vector of the FFLO phase, the resonance peaks appear in the field-direction dependence of the onset of superconductivity.
-
At the resonance, the interplay between the orbital and paramagnetic effects may result in an FFLO lock-in effect.
These effects can open up new possibilities to an unambiguous evidence spatially modulated superconducting phase in quasi-low-dimensional conductors.

Acknowledgments

We acknowledge the support by the Agence Nationale de la Recherche (ANR) and the Deutsche Forschungsgemeinschaft (German Research Foundation, DFG) grant ”Fermi-NESt”.

Author Contributions

Both authors contributed equally to this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of the Extended Nonlinear Lowerence–Doniach Equation

We start from the Eilenberger equations for a quasi-2D superconductor [1]
Ω n + 1 2 Π ^ Q f ω n p , r , k z , p z = Δ r , k z g ω n p , r , k z , p z ,
Ω n 1 2 Π ^ Q f ω n p , r , k z , p z = Δ * r , k z g ω n p , r , k z , p z ,
1 2 Π ^ 0 g ω n p , r , k z , p z = Δ * r , k z f ω n p , r , k z , p z f ω n p , r , k z , p z Δ r , k z ,
where Ω n ω n + 1 / 2 τ i h sign ( ω n ) and the definition of Π ^ is provided in Equation (3) or
Π ^ = v F x , y + i 4 t sin p z d sin χ A
with χ A Qr k z 2 d . We have assumed that the vector potential varies slowly at the inter-plane distances (this assumption means that we neglect the diamagnetic screening currents and take the magnetic field as uniform and given by the external field, H). These equations depend on k z , which takes into account the dependence of the order parameter on the center-of-mass coordinate, z + z / 2 . The order parameter is defined self-consistently as
1 λ Δ r , k z = 2 π T Re ω > 0 f ω n p , r , k z , p z ,
where λ is the pairing constant and the brackets definition is in Equation (7).
The derivation is performed iteratively. It is convenient to write the procedure as the following iterative scheme:
f ω k + 1 n p , r , k z , p z = Δ Ω n g ω k 1 2 Ω n Π ^ f ω k ,
f ω k + 1 n p , r , k z , p z = Δ * Ω n g ω k + 1 2 Ω n Π ^ f ω k .
The zero order approximation (the absence of the superconducting state) results in g ω 0 n p , r , k z , p z = 1 and f ω 0 n p , r , k z , p z = 0 . Then, for the first order approximation, the iterative scheme produces f ω 1 n p , r , k z , p z = Δ / Ω n , f ω 1 n p , r , k z , p z = Δ * / Ω n . For the iterative scheme, we need the iterative expression for g ω k n p , r , k z , p z function. If we write it, for example for g ω 5 , as g ω 5 = g ˜ ω 0 + λ g ˜ ω 1 + λ ω 2 g ˜ ω 2 + λ ω 3 g ˜ ω 2 + λ 4 g ˜ ω 4 + λ 5 g ˜ ω 5 and f ω 5 = f ˜ ω 0 + λ f ˜ ω 1 + λ 2 f ˜ ω 2 + λ 3 f ˜ ω 3 + λ 4 f ˜ ω 4 + λ 5 f ˜ ω 5 , then making use of the normalization condition
g ω 2 + f ω f ω = 1 ,
we obtain the following system of coupled expressions for the k-th iteration
g ˜ ω 0 n p , r , k z , p z = 1 ,
i = 0 k g ˜ ω i g ˜ ω k i = i = 0 k f ˜ ω k i f ˜ ω i .
The second iteration of the iteration scheme leads to
f ω 2 n p , r , k z , p z = Δ r , k z Ω n 1 2 Ω n 2 Π ^ Δ r , k z ,
f ω 2 n p , r , k z , p z = Δ * r , k z Ω n + 1 2 Ω n 2 Π ^ Δ * r , k z .
Therefore, here, f ˜ ω 2 = Π ^ Δ r , k z / 2 Ω n 2 , f ˜ ω 2 = Π ^ Δ * r , k z / 2 Ω n 2 . From the normalization condition, we can also obtain:
g ˜ ω 1 , 2 + 2 g ˜ ω 0 g ˜ ω 2 = f ˜ ω 2 f ˜ ω 0 f ˜ ω 1 f ˜ ω 1 f ˜ ω 0 f ˜ ω 2 ,
2 g ˜ ω 0 g ˜ ω 2 = f ˜ ω 1 f ˜ ω 1 ,
which gives g ˜ ω 2 = f ˜ ω 1 f ˜ ω 1 / 2 = f ω 1 f ω 1 / 2 and g ˜ ω 1 = 0 . Therefore, the normal Green function g ω 2 n p , r , k z , p z acquires the following form
g ω 2 n p , r , k z , p z = g ˜ ω 0 + g ˜ ω 1 + g ˜ ω 2 = 1 f ω 1 f ω 1 2 = 1 Δ r , k z 2 2 Ω n 2 .
Following similar steps, then performing averaging over p y and p z and Fermi surface v F x , we obtain, on the 5th iteration for f ω 5 f ω 5 n p , r , k z , p z ,
f ω 5 = Δ r , k z Ω n Δ x Δ r , k z 2 2 Ω n 3 + 1 4 Ω n 3 Π ^ 2 Δ r , k z + 3 Δ r , k z Δ r , k z 4 8 Ω n 5 + 1 16 Ω n 5 Π ^ 4 Δ r , k z 4 8 Ω n 5 Δ r , k z 2 Π ^ 2 Δ r , k z 1 8 Ω n 5 Δ 2 r , k z Π ^ 2 Δ * r , k z 3 8 Ω n 5 Δ * r , k z Π ^ Δ r , k z Π ^ Δ r , k z 2 8 Ω n 5 Δ r , k z Π ^ Δ r , k z Π ^ Δ * r , k z .
Substituting the found averages into the self-consistency relation given by Δ 5 r , k z / λ = π T n f ω 5 n p , r , k z , p z and using the standard regularization rule
1 λ = ln T T c 0 + 2 π T ω n > 0 1 ω n ,
we obtain the following extended version of the Ginzburg–Landau equation (for clarity, we dropped down k z -dependence), if, by Π ^ , one understands its 3D version,
0 = Δ r π K 1 T K 1 0 T ln T T c 0 + π K 3 T 4 Π ^ 2 Δ r π K 3 T 2 Δ r Δ r 2 + 3 π K 5 T 8 Δ r Δ r 4 + π K 5 T 16 Π ^ 4 Δ r π K 5 T 8 4 Δ r 2 Π ^ 2 Δ r + Δ 2 r Π ^ 2 Δ * r + 3 Δ * r Π ^ Δ r Π ^ Δ r + 2 Δ r Π ^ Δ r Π ^ Δ * r ,
where we have introduced the following notations [86] K m T = 2 T ω > 0 ω n i μ B H m and K m 0 T = 2 T ω > 0 ω n m . This result is similar to the one obtained by Houzet and Mineev [13]. After the averaging procedure, we can write (neglecting the k z dependence)
Δ r ln T T c 0 = Δ r π K 1 T K 1 0 T + π K 3 T 4 2 Δ r π K 3 T 2 Δ r Δ r 2 2 t 2 π K 3 T sin 2 χ A Δ r + 3 π K 5 T 8 Δ r Δ r 4 + π K 5 T 16 4 Δ r + t 2 π K 5 T 4 1 7 cos 2 χ A χ A 2 Δ r + 10 t 2 π K 5 T sin 2 χ A Δ r Δ r 2 t 2 π K 5 T 3 sin 2 χ A χ A Δ r + sin 2 χ A 2 χ A Δ r + 3 sin 2 χ A 2 Δ r π K 5 T 2 Δ r 2 2 Δ r π K 5 T 8 Δ 2 r 2 Δ * r 3 π K 5 T 8 Δ * r Δ r 2 π K 5 T 4 Δ r Δ r Δ * r + 6 t 4 π K 5 T sin 4 χ A Δ r .
Here, the derivative signs mean the following 2 Δ r = v F x 2 x 2 + v F y 2 y 2 Δ r , 4 Δ r = v F x 4 x 4 + v F x 2 v F y 2 x 2 y 2 + v F y 4 y 4 Δ r , χ A Δ r = v F x 2 x χ A x Δ r + v F y 2 y χ A y Δ r and χ A 2 = v F x 2 x χ A 2 + v F y 2 y χ A 2 . In a quasi-1D case, this equation reads as
α Δ x = β x 2 Δ x δ x 4 Δ x 2 γ Δ x Δ x 2 3 ν Δ x Δ x 4 + μ Δ x 2 x 2 Δ x + μ 4 Δ 2 x x 2 Δ * x + 3 μ 4 Δ * x x Δ x 2 + μ 2 Δ x x Δ x x Δ * x t 2 μ 2 1 7 × cos 2 χ A x χ A 2 Δ x 80 t 2 ν sin 2 χ A × Δ x Δ x 2 + 2 t 2 μ 3 sin 2 χ A x χ A x Δ x + sin 2 χ A x 2 χ A Δ x + 3 sin 2 χ A x 2 Δ x 8 t 2 γ sin 2 χ A Δ x 48 t 4 ν sin 4 χ A Δ x ,
where α ln T T c 0 π K 1 T K 1 0 T , β π 4 K 3 T v F x 2 , δ π 16 K 5 T v F x 4 , γ π 4 K 3 T , μ π 2 K 5 T v F x 2 , ν π 8 K 5 T . Here, the derivative signs mean the following 2 Δ x = v F x 2 x 2 Δ x , 4 Δ x = v F x 4 x 4 Δ x , χ A Δ x = v F x 2 x χ A x Δ x and χ A 2 = v F x 2 x χ A 2 . The linearized version of this equation,
α Δ x = β x 2 Δ x δ x 4 Δ x 48 t 4 ν sin 4 χ A Δ x t 2 2 μ 1 7 cos 2 χ A x χ A 2 Δ x + + 2 t 2 μ 3 sin 2 χ A x χ A x Δ x + 3 sin 2 χ A x 2 Δ x 8 t 2 γ sin 2 χ A Δ x
is used in Section 5 of the review.

References

  1. Kopnin, N.B. Theory of Nonequilibrium Superconductivity; Clarendon Press: Oxford, UK, 2001. [Google Scholar]
  2. Clogston, A.K. Upper Limit for the Critical Field in Hard Superconductors. Phys. Rev. Lett. 1962, 9, 266–267. [Google Scholar] [CrossRef]
  3. Chandrasekhar, B.S. A Note on the Maximum Critical Field of High-Field Superconductors. Appl. Phys. Lett. 1962, 1, 7. [Google Scholar] [CrossRef]
  4. Zwicknagl, G.; Wosnitza, J. Breaking Translational Invariance by Population Imbalance: The Fulde–Ferrell–Larkin–Ovchinnikov States. In BCS: 50 Years; Cooper, L.N., Feldman, D., Eds.; World Scientific Publishing: Singapore, 2011. [Google Scholar]
  5. Matsuda, Y.; Shimahara, H. Fulde–Ferrell–Larkin–Ovchinnikov State in Heavy Fermion Superconductors. J. Phys. Soc. Jpn. 2007, 76, 051005. [Google Scholar] [CrossRef]
  6. Larkin, A.I.; Ovchinnikov, Y.N. Inhomogeneous State of Superconductors. Sov. Phys. JETP 1965, 20, 762–769. [Google Scholar]
  7. Fulde, P.; Ferrell, R.A. Superconductivity in a Strong Spin-Exchange Field. Phys. Rev. 1964, 135, A550. [Google Scholar] [CrossRef]
  8. Agterberg, D.F.; Yang, K. The effect of impurities on Fulde–Ferrell–Larkin–Ovchinnikov superconductors. J. Phys. Condens. Matter 2001, 13, 9259. [Google Scholar] [CrossRef]
  9. Gruenberg, L.W.; Gunther, L. Fulde–Ferrell Effect in Type-II Superconductors. Phys. Rev. Lett. 1966, 16, 996. [Google Scholar] [CrossRef]
  10. Aslamazov, L.G. Influence of Impurities on the Existence of an Inhomogeneous State in a Ferromagnetic Superconductor. Sov. Phys. JETP 1969, 28, 773. [Google Scholar]
  11. Takada, S. Superconductivity in a Molecular Field. II: Stability of Fulde–Ferrel Phase. Prog. Theor. Phys. 1970, 43, 27–38. [Google Scholar] [CrossRef]
  12. Adachi, H.; Ikeda, R. Effects of Pauli paramagnetism on the superconducting vortex phase diagram in strong fields. Phys. Rev. B 2003, 68, 184510. [Google Scholar] [CrossRef]
  13. Houzet, M.; Mineev, V.P. Interplay of paramagnetic, orbital, and impurity effects on the phase transition of a normal metal to the superconducting state. Phys. Rev. B 2006, 74, 144522. [Google Scholar] [CrossRef]
  14. Cui, Q.; Yang, K. Fulde–Ferrell–Larkin–Ovchinnikov state in disordered s-wave superconductors. Phys. Rev. B 2008, 78, 054501. [Google Scholar] [CrossRef]
  15. Shimahara, H. Structure of the Fulde–Ferrell–Larkin–Ovchinnikov State in Two-Dimensional Superconductors. J. Phys. Soc. Jpn. 1998, 67, 736–739. [Google Scholar] [CrossRef]
  16. Shimahara, H. Enhancement of the Upper Critical Field Due to a Fermi-Surface Effect in Quasi-Two-Dimensional Superconductors in Parallel Magnetic Fields. J. Phys. Soc. Jpn. 1998, 68, 3069–3073. [Google Scholar] [CrossRef]
  17. Shimahara, H.; Moriwake, K. Temperature Dependence of the Upper Critical Field of Anisotropic Singlet Superconductivity in a Square Lattice Tight-Binding Model in Parallel Magnetic Fields. J. Phys. Soc. Jpn. 2002, 71, 1234–1237. [Google Scholar] [CrossRef]
  18. Buzdin, A.I.; Bulaevskiĭ, L.N. Organic superconductors. Sov. Phys. Usp. 1984, 27, 830. [Google Scholar] [CrossRef]
  19. Beyer, R.; Wosnitza, J. Emerging evidence for FFLO states in layered organic superconductors. Low Temp. Phys. 2013, 39, 225–231. [Google Scholar] [CrossRef]
  20. Agosta, C.C.; Bishop-Van Horn, L.; Newman, M. The Signature of Inhomogeneous Superconductivity. Low Temp. Phys. 2016, 185, 220–229. [Google Scholar] [CrossRef]
  21. Piazza, F.; Zwerger, W.; Strack, P. FFLO strange metal and quantum criticality in two dimensions: Theory and application to organic superconductors. Phys. Rev. B 2016, 93, 085112. [Google Scholar] [CrossRef]
  22. Tanatar, M.A.; Ishiguro, T.; Tanaka, H.; Kobayashi, H. Magnetic field–temperature phase diagram of the quasi-two-dimensional organic superconductor λ-(BETS)2GaCl4 studied via thermal conductivity. Phys. Rev. B 2002, 66, 134503. [Google Scholar] [CrossRef]
  23. Coniglio, W.A.; Winter, L.E.; Cho, K.; Agosta, C.C.; Fravel, B.; Montgomery, L.K. Superconducting phase diagram and FFLO signature in λ-(BETS)2GaCl4 from rf penetration depth measurements. Phys. Rev. B 2011, 83, 224507. [Google Scholar] [CrossRef]
  24. Bergk, B.; Demuer, A.; Sheikin, I.; Wang, Y.; Wosnitza, J.; Nakazawa, Y.; Lortz, R. Magnetic torque evidence for the Fulde–Ferrell–Larkin–Ovchinnikov state in the layered organic superconductor κ-(BEDT-TTF)2Cu (NCS)2. Phys. Rev. B 2011, 83, 064506. [Google Scholar] [CrossRef]
  25. Tsuchiya, S.; Yamada, J.; Sugii, K.; Graf, D.; Brooks, J.S.; Terashima, T.; Uji, S. Phase Boundary in a Superconducting State of κ-(BEDT-TTF)2Cu(NCS)2: Evidence of the Fulde–Ferrell–Larkin–Ovchinnikov Phase. J. Phys. Soc. Jpn. 2015, 84, 034703. [Google Scholar] [CrossRef]
  26. Wright, J.A.; Green, E.; Kuhns, P.; Reyes, A.; Brooks, J.; Schlueter, J.; Kato, R.; Yamamoto, H.; Kobayashi, M.; Brown, S.E. Zeeman-Driven Phase Transition within the Superconducting State of κ-(BEDT-TTF)2Cu(NCS)2. Phys. Rev. Lett. 2011, 107, 087002. [Google Scholar] [CrossRef] [PubMed]
  27. Mayaffre, H.; Krämer, S.; Horvatić, M.; Berthier, C.; Miyagawa, K.; Kanoda, K.; Mitrović, V.F. Evidence of Andreev bound states as a hallmark of the FFLO phase in κ-(BEDT-TTF)2Cu(NCS)2. Nat. Phys. 2014, 10, 928–932. [Google Scholar] [CrossRef]
  28. Buzdin, A.I.; Tugushev, V.V. Phase diagrams of electronic and superconducting transitions to soliton lattice states. Zh. Eksp. Teor. Fiz. 1983, 85, 735–745. [Google Scholar]
  29. Vorontsov, A.B.; Sauls, J.A.; Graf, M.J. Phase diagram and spectroscopy of Fulde–Ferrell–Larkin–Ovchinnikov states of two-dimensional d-wave superconductors. Phys. Rev. B 2005, 72, 184501. [Google Scholar] [CrossRef]
  30. Agosta, C.C.; Jin, J.; Coniglio, W.A.; Smith, B.E.; Cho, K.; Stroe, I.; Martin, C.; Tozer, S.W.; Murphy, T.P.; Palm, E.C.; et al. Experimental and semiempirical method to determine the Pauli-limiting field in quasi-two-dimensional superconductors as applied to κ-(BEDT-TTF)2Cu(NCS)2: Strong evidence of a FFLO state. Phys. Rev. B 2012, 85, 214514. [Google Scholar] [CrossRef]
  31. Koutroulakis, G.; Kühne, H.; Schlueter, J.A.; Wosnitza, J.; Brown, S.E. Microscopic study of the Fulde–Ferrell–Larkin–Ovchinnikov state in an all-organic superconductor. Phys. Rev. Lett. 2016, 116, 067003. [Google Scholar] [CrossRef] [PubMed]
  32. Lortz, R.; Wang, Y.; Demuer, A.; Böttger, P.H.M.; Bergk, B.; Zwicknagl, G.; Nakazawa, Y.; Wosnitza, J. Calorimetric Evidence for a Fulde-Ferrell-Larkin-Ovchinnikov Superconducting State in the Layered Organic Superconductorκ-(BEDT-TTF)2Cu(NCS)2. Phys. Rev. Lett. 2007, 99, 187002. [Google Scholar] [CrossRef] [PubMed]
  33. Beyer, R.; Bergk, B.; Yasin, S.; Schlueter, J.A.; Wosnitza, J. Angle-Dependent Evolution of the Fulde–Ferrell–Larkin–Ovchinnikov State in an Organic Superconductor. Phys. Rev. Lett. 2012, 109, 027003. [Google Scholar] [CrossRef] [PubMed]
  34. Uji, S.; Kodama, K.; Sugii, K.; Terashima, T.; Takahide, Y.; Kurita, N.; Tsuchiya, S.; Kimata, M.; Kobayashi, A.; Zhou, B.; et al. Magnetic torque studies on FFLO phase in magnetic-field-induced organic superconductor λ-(BETS)2FeCl. Phys. Rev. B 2012, 85, 174530. [Google Scholar] [CrossRef]
  35. Cho, K.; Kim, H.; Tanatar, M.A.; Song, Y.J.; Kwon, Y.S.; Coniglio, W.A.; Agosta, C.C.; Gurevich, A.; Prozorov, R. Anisotropic upper critical field and possible Fulde–Ferrel–Larkin–Ovchinnikov state in the stoichiometric pnictide superconductor LiFeAs. Phys. Rev. B 2011, 83, 060502(R). [Google Scholar] [CrossRef]
  36. Zocco, D.A.; Grube, K.; Eilers, F.; Wolf, T.; Löhneysen, H.v. Pauli-Limited Multiband Superconductivity in KFe2As2. Phys. Rev. Lett. 2013, 111, 057007. [Google Scholar] [CrossRef] [PubMed]
  37. Uji, S.; Terashima, T.; Konoike, T.; Yamaguchi, T.; Yasuzuka, S.; Kobayashi, A.; Zhou, B. Internal field effect on vortex states in the layered organic superconductor λ-(BETS)2Fe1-xGaxCl4 (x=0.37). Phys. Rev. B 2017, 95, 165133. [Google Scholar] [CrossRef]
  38. Uji, S.; Kodama, K.; Sugii, K.; Terashima, T.; Yamaguchi, T.; Kurita, N.; Tsuchiya, S.; Konoike, T.; Kimata, M.; Kobayashi, A.; et al. Vortex Dynamics and Diamagnetic Torque Signals in Two Dimensional Organic Superconductor λ-(BETS)2GaCl4. J. Phys. Soc. Jpn. 2015, 84, 104709. [Google Scholar] [CrossRef]
  39. Agosta, C.C.; Fortune, N.A.; Hannahs, S.T.; Gu, S.; Liang, L.; Park, J.-H.; Schleuter, J.A. Calorimetric Measurements of Magnetic-Field-Induced Inhomogeneous Superconductivity above the Paramagnetic Limit. Phys. Rev. Lett. 2017, 118, 267001. [Google Scholar] [CrossRef] [PubMed]
  40. Yonezawa, S.; Kusaba, S.; Maeno, Y.; Auban-Senzier, P.; Pasquier, C.; Jérome, D. Anomalous In-Plane Anisotropy of the Onset of Superconductivity in (TMTSF)2ClO4. Phys. Rev. Lett. 2008, 100, 117002. [Google Scholar] [CrossRef] [PubMed]
  41. Yonezawa, S.; Kusaba, S.; Maeno, Y.; Auban-Senzier, P.; Pasquier, C.; Jérome, D. Magnetic-Field Variations of the Pair-Breaking Effects of Superconductivity in (TMTSF)2ClO4. J. Phys. Soc. Jpn. 2008, 77, 054712. [Google Scholar] [CrossRef]
  42. Croitoru, M.D.; Houzet, M.; Buzdin, A.I. In-Plane Magnetic Field Anisotropy of the Fulde–Ferrell–Larkin–Ovchinnikov State in Layered Superconductors. Phys. Rev. Lett. 2012, 108, 207005. [Google Scholar] [CrossRef] [PubMed]
  43. Croitoru, M.D.; Houzet, M.; Buzdin, A.I. Field-Direction Dependence of the Upper Critical Field in Organic Superconductors. J. Supercond. Nov. Magn. 2012, 25, 1283–1287. [Google Scholar] [CrossRef]
  44. Croitoru, M.D.; Buzdin, A.I. Orientational Effect of the in-Plane Magnetic Field on the FFLO Modulation in Layered Superconductors. J. Supercond. Nov. Magn. 2013, 26, 1657–1661. [Google Scholar] [CrossRef]
  45. Tedrow, P.M.; Meservey, R. Spin-Paramagnetic Effects in Superconducting Aluminum Films. Phys. Rev. B 2013, 8, 5098. [Google Scholar] [CrossRef]
  46. Adams, P.W.; Nam, H.; Shih, C.K.; Catelani, G. Zeeman-limited superconductivity in crystalline Al films. Phys. Rev. B 2017, 95, 094520. [Google Scholar] [CrossRef]
  47. Buzdin, A.I. Proximity effects in superconductor-ferromagnet heterostructures. Rev. Mod. Phys. 2005, 77, 935. [Google Scholar] [CrossRef]
  48. Buzdin, A.I.; Bulaevskii, L.N.; Panyukov, S.V. Critical-current oscillations as a function of the exchange field and thickness of the ferromagnetic metal (F) in an S-F-S Josephson junction. JETP Lett. 1982, 35, 178–180. [Google Scholar]
  49. Ryazanov, V.V.; Oboznov, V.A.; Rusanov, A.Y.; Veretennikov, A.V.; Golubov, A.A.; Aarts, J. Coupling of two superconductors through a ferromagnet: Evidence for a π junction. Phys. Rev. Lett. 2001, 86, 2427. [Google Scholar] [CrossRef] [PubMed]
  50. Oboznov, V.A.; Bol’ginov, V.V.; Feofanov, A.K.; Ryazanov, V.V.; Buzdin, A.I. Thickness dependence of the Josephson ground states of superconductor-ferromagnet-superconductor junctions. Phys. Rev. Lett. 2006, 96, 197003. [Google Scholar] [CrossRef] [PubMed]
  51. Bianchi, A.; Movshovich, R.; Capan, C.; Pagliuso, P.G.; Sarrao, J.L. Possible Fulde–Ferrell– Larkin–Ovchinnikov Superconducting State in CeCoIn5. Phys. Rev. Lett. 2003, 91, 187004. [Google Scholar] [CrossRef] [PubMed]
  52. Martin, C.; Agosta, C.C.; Tozer, S.W.; Radovan, H.A.; Palm, E.C.; Murphy, T.P.; Sarrao, J.L. Evidence for the Fulde–Ferrell–Larkin–Ovchinnikov state in CeCoIn5 from penetration depth measurements. Phys. Rev. B 2005, 71, 020503(R). [Google Scholar] [CrossRef]
  53. Mercure, J.-F.; Bangura, A.F.; Xu, X.; Wakeham, N.; Carrington, A.; Walmsley, P.; Greenblatt, M.; Hussey, N.E. Upper critical magnetic field far above the paramagnetic pair-breaking limit of superconducting one-dimensional Li0.9Mo6O17 single crystals. Phys. Rev. Lett. 2012, 108, 187003. [Google Scholar] [CrossRef] [PubMed]
  54. Tokiwa, Y.; Bauer, E.D.; Gegenwart, P. Quasiparticle entropy in the high-field superconducting phase of CeCoIn5. Phys. Rev. Lett. 2012, 109, 116402. [Google Scholar] [CrossRef] [PubMed]
  55. Kim, D.Y.; Lin, S.-Z.; Weickert, F.; Kenzelmann, M.; Bauer, E.D.; Ronning, F.; Thompson, J.D.; Movshovich, R. Intertwined Orders in Heavy-Fermion Superconductor CeCoIn5. Phys. Rev. X 2016, 6, 041059. [Google Scholar] [CrossRef]
  56. Hatakeyama, Y.; Ikeda, R. Antiferromagnetic Order Oriented by Fulde–Ferrell–Larkin–Ovchinnikov Superconducting Order. Phys. Rev. B 2015, 91, 094504. [Google Scholar] [CrossRef]
  57. Liao, Y.; Rittner, A.S.C.; Paprotta, T.; Li, W.; Partridge, G.B.; Hulet, R.G.; Baur, S.K.; Mueller, E.J. Spin-imbalance in a one-dimensional Fermi gas. Nature 2010, 467, 567–569. [Google Scholar] [CrossRef] [PubMed]
  58. Casalbuoni, R.; Nardulli, G. Inhomogeneous superconductivity in condensed matter and QCD. Rev. Mod. Phys. 2004, 76, 263. [Google Scholar] [CrossRef]
  59. Alford, M.; Bowers, J.A.; Rajagopal, K. Crystalline color superconductivity. Phys. Rev. D 2001, 63, 074016. [Google Scholar] [CrossRef]
  60. Yang, K.; Agterberg, D.F. Josephson effect in Fulde–Ferrell–Larkin–Ovchinnikov superconductors. Phys. Rev. Lett. 2000, 84, 4970. [Google Scholar] [CrossRef] [PubMed]
  61. Gurevich, A. Upper critical field and the Fulde–Ferrel–Larkin–Ovchinnikov transition in multiband superconductors. Phys. Rev. B 2010, 82, 184504. [Google Scholar] [CrossRef]
  62. Dutta, S.; Mueller, E.J. Dimensional crossover in a spin-imbalanced Fermi gas. Phys. Rev. B 2016, 94, 063627. [Google Scholar] [CrossRef]
  63. Tsuzuki, T. On the long-range order in superconducting intercalated layer compounds. J. Low Temp. Phys. 1972, 9, 525–538. [Google Scholar] [CrossRef]
  64. Dzyaloshinskii, I.E.; Kats, E.I. Theory of superconductivity of quasi one-dimensional structures. Zh. Eksp. Teor. Fiz. 1968, 55, 2373–2375. [Google Scholar]
  65. Brison, J.P.; Keller, N.; Vernière, A.; Lejay, P.; Schmidt, L.; Buzdin, A.; Flouquet, J.; Julian, S.R.; Lonzarich, G.G. Anisotropy of the upper critical field in URu2Si2 and FFLO state in antiferromagnetic superconductors. Physica C 1995, 250, 128–138. [Google Scholar] [CrossRef]
  66. Eilenberger, G. Transformation of Gorkov’s equation for type II superconductors into transport-like equations. Z. Phys. A HadroN. Nucl. 1968, 214, 195–213. [Google Scholar] [CrossRef]
  67. Swidzinsky, A.V. Spatially Inhomogeneous Problems in the Theory of Superconductivity; Izdatel’stvo Nauka: Moscow, Russia, 1982. (In Russian) [Google Scholar]
  68. Denisov, D.; Buzdin, A.; Shimahara, H. Types of Fulde–Ferrell–Larkin–Ovchinnikov states induced by anisotropy effects. Phys. Rev. B 2009, 79, 064506. [Google Scholar] [CrossRef]
  69. Dao, V.H.; Denisov, D.; Buzdin, A.; Brison, J.-P. Role of crystal anisotropy on the vortex state in the Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase. Phys. Rev. B 2013, 87, 174509. [Google Scholar] [CrossRef]
  70. Buzdin, A.; Matsuda, Y.; Shibauchi, T. FFLO state in thin superconducting films. Europhys. Lett. 2007, 80, 67004. [Google Scholar] [CrossRef]
  71. Maki, H.; Won, H. The sine-wave like d-wave superconductivity in high magnetic fields. Czechoslov. J. Phys. 1996, 46, 1035. [Google Scholar] [CrossRef]
  72. Croitoru, M.D.; Buzdin, A.I. The Fulde–Ferrell–Larkin–Ovchinnikov state in layered d-wave superconductors: in-plane anisotropy and resonance effects in the angular dependence of the upper critical field. J. Phys. Condens. Matter 2013, 25, 125702. [Google Scholar] [CrossRef] [PubMed]
  73. Saint-James, D.; Thomas, E.J.; Sarma, G. Type II Superconductivity; Pergamon Press: Oxford, UK, 1969. [Google Scholar]
  74. Singleton, J.; Goddard, P.A.; Ardavan, A.; Harrison, N.; Blundell, S. J.; Schlueter, J. A.; Kini, A. M. Test for Interlayer Coherence in a Quasi-Two-Dimensional Superconductor. Phys. Rev. Lett. 2002, 88, 037001. [Google Scholar] [CrossRef] [PubMed]
  75. Müller, J.; Lang, M.; Helfrich, R.; Steglich, F.; Sasaki, T. High-resolution ac-calorimetry studies of the quasi-two-dimensional organic superconductor κ-(BEDT-TTF)2Cu(NCS)2. Phys. Rev. B 2002, 65, 140509. [Google Scholar] [CrossRef]
  76. Izawa, K.; Yamaguchi, H.; Sasaki, T.; Matsuda, Y. Superconducting Gap Structure of κ-(BEDT-TTF)2Cu(NCS)2 Probed by Thermal Conductivity Tensor. Phys. Rev. Lett. 2001, 88, 027002. [Google Scholar] [CrossRef] [PubMed]
  77. Croitoru, M.D.; Buzdin, A.I. Extended Lawrence–Doniach model: The temperature evolution of the in-plane magnetic field anisotropy. Phys. Rev. B 2012, 86, 224508. [Google Scholar] [CrossRef]
  78. Yoshimura, M.; Shigekawa, H.; Nejoh, H.; Saito, G.; Saito, Y.; Kawazu, A. Electronic structure of the organic superconductor κ-(BEDT-TTF)2Cu(NCS)2 [where BEDT-TTF is bis(ethylenedithio)tetrathiafulvalene] studied by scanning tunneling microscopy. Phys. Rev. B 1991, 43, 13590. [Google Scholar] [CrossRef]
  79. Bulaevskii, L.; Buzdin, A.; Maley, M. Intrinsic Pinning of Vortices as a Direct Probe of the Nonuniform Larkin–Ovchinnikov–Fulde–Ferrell State in Layered Superconductors. Phys. Rev. Lett. 2003, 90, 067003. [Google Scholar] [CrossRef] [PubMed]
  80. Uji, S.; Terashima, T.; Nishimura, M.; Takahide, Y.; Konoike, T.; Enomoto, K.; Cui, H.; Kobayashi, H.; Kobayashi, A.; Tanaka, H. et al.; et al. Vortex Dynamics and the Fulde–Ferrell–Larkin–Ovchinnikov State in a Magnetic-Field-Induced Organic Superconductor. Phys. Rev. Lett. 2006, 97, 157001. [Google Scholar] [CrossRef] [PubMed]
  81. Croitoru, M.D.; Buzdin, A.I. Resonance in-plane magnetic field effect as a means to reveal the Fulde–Ferrell–Larkin–Ovchinnikov state in layered superconductors. Phys. Rev. B 2012, 86, 064507. [Google Scholar] [CrossRef]
  82. Bak, P. Commensurate phases, incommensurate phases and the devil’s staircase. Rep. Prog. Phys. 1982, 45, 587. [Google Scholar] [CrossRef]
  83. Kontorova, T.A.; Frenkel, Y.I. On the theory of plastic deformation and twinning. II. Zh. Eksp. Teor. Fiz. 1938, 8, 1340–1348. [Google Scholar]
  84. Frenkel, Y.I. Introduction into the Theory of Metals; Nauka: Leningrad, Russia, 1972. (In Russian) [Google Scholar]
  85. Croitoru, M.D.; Buzdin, A.I. Soliton formation in the FFLO phase. Phys. Rev. B 2016, 94, 214512. [Google Scholar] [CrossRef]
  86. Buzdin, A.I.; Kachkachi, H. Generalized Ginzburg–Landau theory for nonuniform FFLO superconductors. Phys. Lett. A 1982, 225, 341–348. [Google Scholar] [CrossRef]
  87. Buzdin, A.I.; Polonskii, S.V. Nonuniform state in quasi-1D superconductors. Zh. Eksp. Teor. Fiz. 1987, 93, 747–761. [Google Scholar]
  88. Croitoru, M.D.; Buzdin, A.I. Peculiarities of the orbital effect in the Fulde–Ferrell–Larkin–Ovchinnikov state in quasi-one-dimensional superconductors. Phys. Rev. B 2014, 89, 224506. [Google Scholar] [CrossRef]
  89. Croitoru, M.D.; Buzdin, A.I. FFLO-wave-vector Lock-in Effect in Quasi-1D Superconductors. J Supercond. Nov. Magn. 2015, 28, 1305–1308. [Google Scholar] [CrossRef]
  90. Dzyaloshinskii, I.E. Theory of helicoidal structures in antiferromagnets. III. Sov. Phys. JETP 1964, 20, 665. [Google Scholar]
  91. Izyumov, Y.A. Modulated, or long-periodic, magnetic structures of crystals. Sov. Phys. Usp. 1984, 27, 845–867. [Google Scholar] [CrossRef]
  92. Izyumov, Y.A. Neutron Diffraction in Long-Periodic Structures; Energoatomizdat: Moscow, Russia, 1987. [Google Scholar]
  93. Frank, F.C.; Van der Merwe, J.H. One-Dimensional Dislocations. I. Static Theory. Proc. Roy. Soc. A 1949, 198, 205–216. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the paramagnetism of conduction electrons. Zeeman effect.
Figure 1. Schematic representation of the paramagnetism of conduction electrons. Zeeman effect.
Condensedmatter 02 00030 g001
Figure 2. Scheme of the quasi-2D layered metal.
Figure 2. Scheme of the quasi-2D layered metal.
Condensedmatter 02 00030 g002
Figure 3. The Fermi surface of the layered metal in the form of a corrugated cylinder.
Figure 3. The Fermi surface of the layered metal in the form of a corrugated cylinder.
Condensedmatter 02 00030 g003
Figure 4. The upper critical field H c P 2 and the absolute value of the FFLO modulation vector q as a function of T c 0 for s-wave (left panel) and d-wave (right panel) superconductors in the pure Pauli limit.
Figure 4. The upper critical field H c P 2 and the absolute value of the FFLO modulation vector q as a function of T c 0 for s-wave (left panel) and d-wave (right panel) superconductors in the pure Pauli limit.
Condensedmatter 02 00030 g004
Figure 5. Contribution of the orbital effect as a function of T c P / T c 0 for several field directions when η = 2 . 6 . The regions showing the results obtained when q is along the x-axis and for the conventional phase are highlighted in yellow, while the white region is for the case when q vector makes angle ϕ = π / 4 with the x-axis.
Figure 5. Contribution of the orbital effect as a function of T c P / T c 0 for several field directions when η = 2 . 6 . The regions showing the results obtained when q is along the x-axis and for the conventional phase are highlighted in yellow, while the white region is for the case when q vector makes angle ϕ = π / 4 with the x-axis.
Condensedmatter 02 00030 g005
Figure 6. Normalized superconducting transition temperature, T c α / T c P as a function of the α for several values of T c P / T c 0 and η = 2 . 6 . The angle between the q-vector and x-axis is ϕ = π / 4 (first and second panels) and ϕ = 0 (last panel). The direction of q-vector is fixed by the symmetry of the order parameter.
Figure 6. Normalized superconducting transition temperature, T c α / T c P as a function of the α for several values of T c P / T c 0 and η = 2 . 6 . The angle between the q-vector and x-axis is ϕ = π / 4 (first and second panels) and ϕ = 0 (last panel). The direction of q-vector is fixed by the symmetry of the order parameter.
Condensedmatter 02 00030 g006
Figure 7. Contribution of the orbital effect as a function of T c P / T c 0 for several field directions when η = 2 . 6 . The case of d-wave with q vector making angle ϕ = π / 4 with the x-axis as well as the s-wave superconductor are shown.
Figure 7. Contribution of the orbital effect as a function of T c P / T c 0 for several field directions when η = 2 . 6 . The case of d-wave with q vector making angle ϕ = π / 4 with the x-axis as well as the s-wave superconductor are shown.
Condensedmatter 02 00030 g007
Figure 8. The same as in Figure 6; however, the direction of q-vector is influenced by the crystal field.
Figure 8. The same as in Figure 6; however, the direction of q-vector is influenced by the crystal field.
Condensedmatter 02 00030 g008
Figure 9. The same as in Figure 6; however, the direction of q-vector is influenced by the crystal field. Here is the s-wave case. Thick lines are for Δ 2 0 ; thin lines are for Δ 2 = 0 .
Figure 9. The same as in Figure 6; however, the direction of q-vector is influenced by the crystal field. Here is the s-wave case. Thick lines are for Δ 2 0 ; thin lines are for Δ 2 = 0 .
Condensedmatter 02 00030 g009
Figure 10. Schematic picture of the FFLO modulated order parameter and the oscillations of the interlayer coupling, explaining the physical origin of the resonances.
Figure 10. Schematic picture of the FFLO modulated order parameter and the oscillations of the interlayer coupling, explaining the physical origin of the resonances.
Condensedmatter 02 00030 g010
Figure 11. The absolute value of the FFLO modulation vector q versus of T c P / T c 0 for an s-wave quasi-1D superconductor when accounting for the orbital effects within the first iteration. Here, t c / T c = 0 . 2 . The solid line takes into account the influence of resonance, and the short dashed line shows the field vector Q.
Figure 11. The absolute value of the FFLO modulation vector q versus of T c P / T c 0 for an s-wave quasi-1D superconductor when accounting for the orbital effects within the first iteration. Here, t c / T c = 0 . 2 . The solid line takes into account the influence of resonance, and the short dashed line shows the field vector Q.
Condensedmatter 02 00030 g011
Figure 12. The phase shift of the FFLO modulation for several values of the anisotropy parameter ζ . The straight line corresponds to an unperturbed incommensurate structure.
Figure 12. The phase shift of the FFLO modulation for several values of the anisotropy parameter ζ . The straight line corresponds to an unperturbed incommensurate structure.
Condensedmatter 02 00030 g012
Figure 13. The spatial distribution of the order parameter in the vicinity of the I–C phase transition for Q = 0 . 2 and for several values of the anisotropy parameter ζ (a) and Q = 0 . 7 (b).
Figure 13. The spatial distribution of the order parameter in the vicinity of the I–C phase transition for Q = 0 . 2 and for several values of the anisotropy parameter ζ (a) and Q = 0 . 7 (b).
Condensedmatter 02 00030 g013

Share and Cite

MDPI and ACS Style

Croitoru, M.D.; Buzdin, A.I. In Search of Unambiguous Evidence of the Fulde–Ferrell–Larkin–Ovchinnikov State in Quasi-Low Dimensional Superconductors. Condens. Matter 2017, 2, 30. https://doi.org/10.3390/condmat2030030

AMA Style

Croitoru MD, Buzdin AI. In Search of Unambiguous Evidence of the Fulde–Ferrell–Larkin–Ovchinnikov State in Quasi-Low Dimensional Superconductors. Condensed Matter. 2017; 2(3):30. https://doi.org/10.3390/condmat2030030

Chicago/Turabian Style

Croitoru, Mihail D., and Alexandre I. Buzdin. 2017. "In Search of Unambiguous Evidence of the Fulde–Ferrell–Larkin–Ovchinnikov State in Quasi-Low Dimensional Superconductors" Condensed Matter 2, no. 3: 30. https://doi.org/10.3390/condmat2030030

APA Style

Croitoru, M. D., & Buzdin, A. I. (2017). In Search of Unambiguous Evidence of the Fulde–Ferrell–Larkin–Ovchinnikov State in Quasi-Low Dimensional Superconductors. Condensed Matter, 2(3), 30. https://doi.org/10.3390/condmat2030030

Article Metrics

Back to TopTop