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
. 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,
, where
is the superconducting gap at
and
[
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
. However, the influence of impurities may reduce
[
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
; (ii) the system should be in a clean limit,
, where
l is the mean-free path of the quasiparticles and
is the superconducting coherence length at
and
[
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)
GaCl
; (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)
Cu(NCS)
[
24,
25]; (iii) evidence for phase transition within the superconducting phase obtained from the local electronic spin polarization in
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
phase diagram that are consistent with the FFLO phase obtained from rf-penetration depth measurements [
30] for organic sample
-(BEDT-TTF)
Cu(NCS)
; (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)
Cu(NCS)
and
-(ET)
SF
CH
CF
SO
[
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)
FeCl
[
34]; (
vi) measurements of the temperature and angular dependencies of
in pnictide superconductor LiFeAs [
35] and KFe
As
[
36] as well as a dip in the interlayer resistance in parallel fields observed in the compound
-(BETS)
Fe
Ga
Cl
for
[
37] (for all types of salts, (
,
,
) 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)
Cu(NCS)
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
of the organic superconductor (TMTSF)
ClO
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
kOe is observed,
kOe [
40]). Moreover, an unusual in-plane anisotropy of
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
, 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
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
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
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.,
, but sufficiently large to make the mean field treatment justified, namely the coupling between layers,
t, satisfies the condition
[
63,
64]. Here,
is the critical temperature of the system at
. 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
(
is a coordinate in
xy-plane), i.e.,
, 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]:
Here, we use the notation
where
with
,
is the Zeeman energy,
is the in-plane Fermi velocity,
is the impurity scattering time, and
. 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
where
with
, the azimuthal angle of
(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
where
is defined as self-consistency relation,
where
is the pairing constant and the brackets denote averaging over
and
,
We assume that the temperature unit is so chosen that the Boltzmann constant
. For the case of an
s-wave pairing interaction,
.
The FFLO state appears as a phase with a modulated order parameter with a wave vector
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
The simple exponential solution
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,
. The solution of Equation (
2) might be written as
Because of the form for
in Equation (
8), one can write the order parameter
as
From symmetry considerations, it follows that
. Substituting Equations (
8) and (
9) into Equation (
2), one gets the following system of coupled equations [
72]:
where
and
. Here, we have taken into account that
and introduced the notation
. This hierarchy of coupled equations is broken on the level of
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
.
3. Angular Dependence of the Upper Critical Field
From Equation (
10) we can obtain function
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
where the renormalized variables
and
are introduced,
is the temperature of the onset of the superconductivity in the pure Pauli regime,
, and
is the angle the
q vector makes from the
x-axis. The
and the modulation vector, maximizing
, are illustrated in
Figure 1. In the case of an
s-wave superconductor (obtained from Equation (
15) by substitution
)
For a
d-wave superconductor, the paramagnetic upper critical field is never smaller than
for an
s-wave superconductor [
71], as shown in
Figure 4. The optimal direction of the modulation vector in a
d-wave superconductor is
for
and
for
with
.
and
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
and
. Therefore, the order of the magnitude of the magnetic field required to observe the FFLO state can be found from the relation
. Taking this into account, one obtains
, where
a is the unit cell in the
xy-plane. Therefore,
. Due to the assumption
, one has
. This condition allows us to hold only the terms up to the first harmonics in Equations (
8) and (
9). Then, one simply gets
The solution of this system of equations is
where the following notations are used
Here, the expression for
can acquire an analytical form
where we have introduced the following notations
,
,
,
,
and
where
k is the cycling index with
.
In our numerical investigations, we restrict ourselves to the parameters of the compound
-(BEDT-TTF)
Cu(NCS)
. The Maki parameter,
[
24], the interlayer coupling is
K (
K) [
74],
(
),
[
75] and the Fermi velocity
[
76]. We have chosen the value
[
74]. Introducing the dimensionless Fermi velocity parameter,
, this value of
corresponds to
. The interlayer distance is
. 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 makes from the x-axis: at high temperatures of the FFLO phase, , while, in the low-temperature phase 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
as a function of
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
makes from the
x-axis. Due to the four-fold symmetry of the
pairing, we plot
curves only for
. The regions
and
are shadowed. In the low-temperature region, the modulation vector
is along the
x-axis, while at
−
, and it describes
in the conventional phase. The central domain illustrates
, when the modulation wave vector
makes the angle
with the
x-axis. The re-orientational transition at
is accompanied by a strong change of the orbital correction, which should result in the sharp step on the
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,
, calculated at
, when the modulation vector is fixed by the symmetry of pairing to the direction of maxima of the order parameter, and at
and
, when it is fixed to the nodes of
. 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
, the transition from the conventional phase to the FFLO modulated phase, accompanied by a rotation of the anisotropy principal axis by the angle
. The second transition occurs at
, when the direction of the FFLO modulation vector rotates by
, 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
dependence at
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
directions and investigate the anisotropy of the superconducting onset in the temperature range
.
Figure 7 displays the orbital-motion induced normalized correction of the transition temperature,
, as a function of normalized temperature,
, 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
, considered in this section, while the solid lines are the solutions with
, discussed in the next section.
As it was intuitively expected, the orbital effects reduce the superconducting onset temperature,
in both cases. While increasing the applied magnetic field,
first decreases in most cases until the tricritical point,
, is reached [
77]. At
the function
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,
may exhibit an upturn and
may approach the paramagnetic limit,
, when
H increases. In contrast, for small
an increase of the magnetic field leads to a decrease of
. For intermediate angles,
can be a non-monotonic function of the field strength. One can infer that the strong field-direction dependence of the superconducting onset temperature,
, 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
the behavior of the upper critical field is isotropic. When increasing
H above
, a strong in-plane anisotropy of
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
and fixed along the direction of the modulation vector. Even in the presence of the “easy axis” along the direction
induced by the crystal field, the four-fold degeneracy of the direction of the FFLO modulation vector restores again for
. 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)
ClO
samples as a realization of the FFLO state with the modulation vector close to the
-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
, where
d is the interlayer distance and
. 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
coincides, i.e., in the resonance with the period of FFLO modulation
(for
,
). The upper critical field at
observed in
-(BEDT-TTF)
Cu(NCS)
is
T. Since in this compound the interlayer distance
nm [
78],
7–9 nm this leads to
nm and
70–90 nm. In (TMTSF)
ClO
T,
nm,
nm, resulting in
nm and
nm. In
-(BETS)
FeO
, the lower
T,
nm,
nm, resulting in
nm and
nm. For these compounds
at
. The FFLO modulation appears only at
with a wave vector
growing from
to
with decreasing temperature. Therefore, if the condition
is satisfied at
, then, at some finite temperature
T, the resonance condition,
, 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)
FeCl
[
80]).
From the system of Equations (
10) and (13), it is seen that, if
then the averaged Equations (
10) and (12) for
and
show that
is of the same order as
. 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
in a second-order approximation on the small parameter
where, additionally to Equation (
20), the following are introduced:
Taking into account the fact that the critical temperature when accounting for the orbital effects,
, is close to
Equations (
25) and (
26) can be written as the following system of coupled equations:
where
The solution of the system of Equations (
29) and (
30) is given as
where, for
the values are chosen that maximize the critical temperature. Usually, the second harmonic of the order parameter,
, can be neglected because the developed theory is valid up to the second order with respect to
, and because
≲
. In this case, the solution of the system is just
, as it was shown in the previous section. However, if
(the resonance condition), the term
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
). Consequently,
becomes on the order of
, and one has to consider both equations in the system on equal footing [
81].
Figure 7 displays the change of the transition temperature,
, 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
. However, there are essential differences in the curves’ behavior for the angles
. These differences are induced by the just discussed resonance effect. The curve
for
exhibits only one cusp at
.
is close to zero and this vicinity induces a wide pseudo-resonance peak on the curve
for
. 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,
, when accounting for the second harmonics of the order parameter,
. 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
, where
and
are the inter-chain distances along the
y- and
z-axis, respectively. We assume that the couplings between chains are small, i.e.,
and
, but sufficiently large to stabilize the superconducting long-range order and to make the mean field treatment justified,
,
[
63]. Here,
is the critical temperature of the system at
. In quasi-1D superconductors, the orbital effect is extremely weak for the magnetic field applied along the chains. Indeed, near
, the parallel (along the
x-axis) upper critical field is
, where
and
while the perpendicular upper critical field (i.e., along the
y-axis) is
with
We see that
and the orbital effect for the parallel component of the magnetic field is weakened by the factor
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
and the corresponding vector potential may be chosen as
,
,
. Then, the Eilenberger equations reduces to the following linearized Ginzburg–Landau equation (for details, see the
Appendix)
where
with
,
,
and
in the uniform superconducting state. This equation contains the magnetic-field induced potential, which is periodic in real space with the period
. Equation (
35) leads to the functional obtained in the limit of
as [
85]
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
, 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,
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
and
, respectively, that provides the expression for the modulation wave vector in the paramagnetic limit,
. This wave vector increases from 0 at the tricritical point to a very large value at
.
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,
, 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
, where
is a slowly varying function with
. This implies that the orbital effect weakly modifies the exact one-dimensional solution, which is justified for
We insert this solution into the functional (
36), perform the average over the period
as
, and taking into account that
, obtain
In the adopted approximation, the functional
depends only on the space varying phase,
. Since we are interested in the behavior of the system in the vicinity of the resonance
, we can write the functional in a simpler form (normalized by
, note that
)
where
is the detuning, or the relative misfit between two periodicities in the system. The phase
is the shift of the FFLO modulation relative to the minima in the potential
. The state
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
we can obtain the ground state minimizing
, which is given by the solutions of the exactly integrable sine-Gordon equation for the phase
where
is the effective interlayer coupling parameter. In the absence of this coupling,
, Equation (
40) has a solution
, which describes an unperturbed one-harmonic incommensurate phase,
. On the other hand, Equation (
40) has a trivial solution,
, that corresponds to a commensurate structure,
. For a finite
Equation (
40) should describe some inhomogeneous distribution of the phase of the order parameter. The first integral of this equation is given by
where
and
is the constant of integration. Then, the exact solution of Equation (
41), expressed in terms of the Jacobi amplitude, is
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
, which favors
and the derivative part that favors
. For
it reduces to
describing a domain wall, which separates two commensurate regions, one with phase
, and the other with
. 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
given by Equation (
42), the normalized energy of the state is expressed in terms of complete elliptic integrals
E an
K as
Minimization of the energy with respect to
leads to the following equation for
where
. Here, the parameter
varies in the range
as long as
v varies in the range
. The variation of the parameter
results in a drastic changes in the behavior of the phase
.
Figure 12 displays the phase for several values of parameter
. As
is increased a plateau section appears within the period of soliton lattice
or
. 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
One can see that, as
parameter
increases its value. In the limit
the period of the soliton lattice diverges
and the modulation wave vector is negligible. For
there are no real valued solutions of Equation (
44). Consequently, there exists a limiting value of
, beyond which the phase is always locked, and the commensurate phase remains stable as
decreases further. This means that the modulation vector coincides with the modulation induced by the field. Therefore, at
we have a lock-in of the FFLO modulation to the field induced potential, i.e., incommensurate/commensurate type phase transition. This corresponds to
and the FFLO modulation is described by
.
Figure 13 illustrates the behavior of the spatially varied order parameter in the vicinity of the incommensurate–commensurate (I–C) phase transition when
. 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
we reduce the quasi-2D case to the just described quasi-1D case, Label (
36), but with the renormalized coefficients:
,
,
,
. 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.