## BCS Stripe Phase in Coupled Bilayer Superconductors

Uddalok Nag, 1 Jonathan Schirmer, 1, 2 Enrico Rossi, 2 C.-X. Liu, 1 and J. K. Jain 1

1 Department of Physics, 104 Davey Lab, The Pennsylvania State University, University Park, Pennsylvania 16802, USA 2 Department of Physics, William &amp; Mary, Williamsburg, Virginia 23187, USA

As a signature of competing correlations, stripes occur in a variety of strongly correlated systems, such as high temperature superconductors (SCs) and quantum Hall effect. We study a double layer SC in the presence of a parallel magnetic field B within the Bogoliubov-de Gennes framework. We find that for low B the system remains in the 'Bardeen-Cooper-Schrieffer (BCS) phase' with a spatially uniform gap, but with increasing B , a transition occurs into a phase which contains stripes of the BCS phase separated by regions where the interlayer phase difference rotates by 2 π due to the presence of inter-layer vortices. This stripe phase is predicted to manifest through oscillations in the amplitude of the SC gap and an alternating pattern of supercurrents. We will comment on the relation to previous works based on the Landau-Ginzburg formalism as well as on the possible experimental realization and signature of this phase.

Introduction -Stripes occur in many strongly correlated systems and are a signature of competing orders. In the quantum Hall effect they occur at a filling of the form ν = n +1 2 where the / ν = n and ν = n +1 phases compete [1, 2], and in high temperature superconductors (SCs) they arise because the holes in an antiferromagnetic background want to cluster together but a competing phase separation has a prohibitive Coulomb energy cost [3]. Here we predict that in a bilayer SC in the presence of a magnetic field parallel to the layers, a state consisting of stripes of the Bardeen-Cooper-Schrieffer (BCS) phase is stabilized in some parameter range.

Liu considered a bilayer Ising SC model in the presence of a parallel magnetic field using the Ginzburg-Landau (GL) approach [4], which is valid near the critical temperature. With increasing magnetic field, a phase transition was predicted from the BCS phase into a phase with the gap function ∆ ± ( x ) = | ∆ | e ± iqx , where ± labels two SC layers. As | ∆ is spatially uniform and each SC | layer has a single q value, we refer to this phase as the 'layer helical phase' below. We note that the layer helical phase is also an example of the Fulde-Ferrell-LarkinOvchinnikov (FFLO) phase [5-8]. In a system with time reversal and spatial inversion symmetry, which ensures that for every electron with the momentum k on the Fermi surface there exists another with the momentum -k and opposite spin, Cooper pairs with zero momentum are stabilized, as assumed in the BCS theory. But if any of these symmetries are broken, for example the time reversal symmetry (due to a magnetic field) or spatial inversion symmetry (due to antisymmetric spin-orbit coupling), the resulting Cooper pairs can have non-zero center of mass momentum q , producing the FFLO superconducting phase described by a gap parameter of the form ∆( r ) = ∆( ) | r | e i q r · . Even though the FFLO phase was predicted six decades ago, its realization has been challenging and it has been reported in very few condensed matter [9, 10] and atomic systems [11-13].

SCs [14, 15], characterized by an upturn in the in-plane upper critical magnetic field when lowering temperature, as well as by broken translational and rotational symmetries [16]. Given the experimental feasibility of this remarkable physics, we theoretically evaluate the phase diagram of a bilayer SC in the presence of an in-plane B by obtaining self-consistent solutions of the Bogoliubovde Gennes (BdG) equations for electrons on a lattice subject to an attractive Hubbard U interaction. The BdG formalism, which is valid for a wider range of parameters than the GL approach and also allows for spatial variation of the complex gap function, reveals striking physics beyond that in Ref. [4]. Most remarkably, as the magnetic field is increased, a transition takes place into a phase that consists of stripes of the BCS phase separated by inter-layer vortices (namely 2 π phase winding of the order parameter between two layers), as shown schematically in the inset of Fig. 1(c). This 'BCS stripe phase' is an example of a highly nonuniform FFLO state. An immediately testable prediction of our work is that the phase occurring close to the transition line in the experiment of Ref. [15] is a BCS stripe phase, the striped character of which can be revealed by a measurement of the local order parameter, whose amplitude should show periodic oscillations in the direction perpendicular to the magnetic field, as well as by the distribution of local supercurrent. This phase continuously evolves into the layer helical phase of Ref. [4] in the limit of large B where the layers become effectively decoupled. An observation of the stripe phase through substantial variation in the gap amplitude will, we believe, provide a more direct confirmation of the underlying FFLO physics.

Indirect evidence for the phase transition predicted in Ref. [4] has been observed in bilayer or multi-layer Ising

Model Hamiltonian and Method We model each layer of the double-layer SC as a square lattice, with intralayer hopping t 1 = 1 and an inter-layer hopping t (inset of Fig. 1(a)). Superconductivity is incorporated through an on-site attractive interaction term, i.e. a negativeU Hubbard model. The Hamiltonian is given by H =

FIG. 1. (a) Transition from BCS phase (green) to stripe phase (black) as a function of Φ / Φ 0 for U = 4, t = 0 6, . µ = -3, N = 60, and T = 0. In the vicinity of the transition we find two self-consistent solutions. The inset shows the side view of the lattice model. (b) Determination of the thermodynamic limit of for the critical flux per plaquette, Φ p,c , for several values of t . (c) B -t phase diagram at T = 0, U = 4, and µ = -3. The error in the thermodynamic value of the critical flux is determined by performing linear fits, as in (b), with different numbers of points. The insets show schematics of the current flow patterns for the two phases. (d) The B -T phase diagram for U = 4, µ = -3, N = 40, and t = 0 5. . The green, red and blue dots show BCS, stripe and normal phases. (e) The average gap amplitude as a function of temperature for three different values of Φ p / Φ . 0 (f) The relative phase at 4 different points in the phase diagram 1(d) as we first move along the transition line from BCS to stripe phase and then along the BCS to normal phase. The plots for T = 0, Φ p = 0 062Φ . 0 and T = 0 12, Φ . p = 0 08Φ . 0 coincide.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00689v1-figure-324df63d-2be8-4011-9dbb-6b20b2cdf69a.jpg)

H

gauge choice does not break the periodicity of the lattice, implying that the in-plane magnetic field does not break the translational symmetry; in other words, we do not need to define a magnetic unit cell (MUC) and the flux Φ can take arbitrary values. As shown in the Supplementary Material, this model describes the low energy physics of bilayer Ising superconductors for which the spin is pinned to the out-of-plane directions and therefore the effect of Zeeman coupling is negligible.

$$\mathcal { H } _ { \mathcal { I } } = - U \sum _ { j, l } c _ { l, j, \uparrow } ^ { \dagger } c _ { l, j, \downarrow } ^ { \dagger } c _ { l, j, \downarrow } c _ { l, j, \uparrow }.$$

$$\mathcal { H } _ { 0 } + \mathcal { H } _ { I } \text{ with } \quad \text{gauge cl} \\ \mathcal { H } _ { 0 } = & - \sum _ { \tilde { \delta }, l, j, \sigma } \left ( e ^ { i A _ { \tilde { \delta }, j } } c _ { l, j + \tilde { \delta }, \sigma } ^ { \dagger } c _ { l, j, \sigma } + e ^ { - i A _ { \tilde { \delta }, j } ^ { \dagger } } c _ { l, j, \sigma } ^ { \dagger } c _ { l, j + \tilde { \delta }, \sigma } \right ) \text{impl/in} } \text{ the trar} \text{ need to} \\ - t \sum _ { j, \sigma } \left ( e ^ { i A _ { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \that \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \ t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t | \t |                                                                                                                                                                                                         }                                                                                                                                                                                                        \end{                                                                                                                                                                                                       
                                                                                                                                                                                                        - t \sum _ { j, l, \sigma } } \left ( e ^ { i A _ { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \t \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde                                                                                                                                                                                                       .                                                                                                                                                                                                        + t \sum _ { j, \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde } \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde [ \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde { \tilde	 \t                                                                                                                                                                   {$$

The layers are taken to be parallel to the x-y plane, l = ± and σ = ↑ ↓ , are the layer and spin indices, j ≡ ( j x , j y ) labels the lattice sites in a plane, and ˆ = ˆ δ x, y ˆ denotes in-plane unit vectors. The operator c l, j ,σ annihilates an electron with spin σ at site j in layer l . The magnetic field is taken as B = By ˆ. It is incorporated through A l δ, j and A z, j , which are the Peierls phases for intraand inter-layer hoppings, determined by the condition that the phases around each plaquette add to 2 π Φ p / Φ , 0 where Φ p is the flux per plaquette and Φ 0 = h/e is the flux quantum. We will work with the gauge choice, A l ˆ x, j = lπ × Φ p / Φ , and 0 A l ˆ y, j = 0 = A z, j . Note that this

The mean field BdG Hamiltonian H MF , bilinear in electron operators, is obtained by replacing the H I term by

$$\stackrel { \sim } { x, j _ { y } } } _ { \text{notes} } \quad \mathcal { H } _ { \mathcal { I } } \rightarrow - \sum _ { l, j } \left ( \Delta _ { l, j } c _ { l, j, \uparrow } ^ { \dagger } c _ { l, j, \downarrow } ^ { \dagger } + \Delta _ { l, j } ^ { * } c _ { l, j, \downarrow } c _ { l, j, \uparrow } - \frac { | \Delta _ { l, j } | ^ { 2 } } { U } \right )$$

$$\Delta _ { l, j } = U \langle c _ { l, j, \downarrow } c _ { l, j, \uparrow } \rangle.$$

where ⟨· · · ⟩ denotes thermal average. We solve for this Hamiltonian for a square system of size N × N in the x-y plane. We take a unit cell of size 1 in the ˆ direction, and y let it span the entire system ( N sites) in the ˆ direction. x So, counting two layers, we have N × 1 × 2 sites in the

FIG. 2. BCS stripe phase for a system of length N = 60 for different values of the magnetic field. (For clarity, we show only half of the system; this behavior repeats in the other half.) For small Φ / Φ 0 the state (not shown) has zero relative phase, spatially uniform gap, equal and opposite currents in two layers, and no interlayer current. The first column corresponds to a magnetic flux Φ just above the critical value, the second to an intermediate Φ, and the third to a high Φ. Spatial variations are shown for: (a-c) the relative phase (difference between the phases in two layers); (d-f) the gap amplitude; (g-i) the phase of ∆ [determined up to a global U(1) rotation]; and (j-l) the current density J x in the top layer. The magnitudes of the gap are same in both layers. The current in the bottom layer is -J x , from which the interlayer current may be deduced. The black dots in (g-i) show the positions of interlayer vortices in the gap parameter (determined up to a global translation). The parameters chosen are U = 4 , µ = -3 , T = 0 K,t = 0 6. . Note that for a large Φ the relative phase varies linearly with the position and the gap amplitude becomes nearly constant.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00689v1-figure-7cd9a328-0e9b-422d-9bd0-ad5f030ec547.jpg)

unit cell. We assume anti-periodic boundary conditions in both ˆ and ˆ direction. x y We perform our calculation for various N values and extrapolate to 1 N → 0 to obtain the thermodynamic limits for various quantities (Fig. 1 (b)).

To go to the momentum space, we write j = R j + αx ˆ, where R j = j y ˆ is y the position vector of the unit cell containing the site j , and α = 1 · · · N is the in-plane site index within the unit cell. Wethen define the transformation: c l, j ,σ = 1 N ∑ k e i k R · j c l,α,σ ( k ), where k = π NL n y y ˆ; n y = -( N -1) , -( N -3) . . . ( N -1) ( n y only takes odd integer values). Defining an 8 N dimensional vector Ψ( k ) = [ { c l,α, ↑ ( k ) } { , c l,α, ↓ ( k ) } { , c † l,α, ↑ ( -k ) } { , c † l,α, ↓ ( -k ) } ] T , we obtain (closely following Ref. [17])

$$H _ { \text{MF} } & \ = \ \frac { 1 } { 2 } \sum _ { k } \Psi ^ { \dagger } ( k ) \hat { \mathcal { H } } _ { \text{BdG} } ( k ) \Psi ( k ) + \frac { 1 } { 2 } \sum _ { k } \text{Tr} [ H _ { 0 } ( k ) ] \quad \begin{array} { c c } \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { u } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \tilde { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu  & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar{ \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } & \bar { \mu } \\ & \bar { \mu } & \bar { \mu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \tilde { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar{ \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } \\ & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar \nu & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } &\bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nolimits & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \Bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu } & \bar { \nu$$

To obtain self-consistent solutions we (i) begin with an initial guess ∆ guess l,α = ∆ | | e ilnπα N on the site α = 1 , · · · N , which corresponds to n vortices in the entire system; (ii) construct the BdG Hamiltonian to calculate ∆ l,α using Eq. (4); (iii) use this as our new guess; and (iv) repeat this process until we get a self consistent result, defined by the condition that the relative difference between the absolute values of the gaps in two successive steps is less than some tolerance (typically taken to be 10 -5 ) at each site. We perform this for all values of 0 ≤ n ≤ Int[2Φ / Φ ] + 1 0 and take the self-consistent solution with the lowest free energy F = ⟨ H MF -TS ⟩ , where the thermal average of the first term on the right hand side of Eq. (5) is given by 1 2 ∑ k ,β E β ( k ) n F ( E β ( k )) with n F ( E ) = ( e E/ k ( B T ) + 1) -1 , and ⟨ S ⟩ = -k B ∑ k,β [( n F ( E β ( k )) log n F ( E β ( k )) +(1 -n F ( E β ( k ))) log(1 -n F ( E β ( k )))]. We have tested the stability of the lowest energy solution by adding a random variation to the initial guess. Our solutions have zero net current, as required by Bloch's theorem [18].

where ˆ H BdG is an 8 N × 8 N BdG matrix and H 0 ( k ) is the single particle part of the Hamiltonian in Fourier space. The eigen-energy spectrum of H MF can be obtained by the diagonalization of ˆ H BdG .

The in-plane current and that between two layers are given by J δ,l ( r j ) = -∂ H ∂A δ,l and J z ( r j ) = -∂ H ∂A z , respectively. The average current densities in the ˆ and ˆ dix z

$$\left \langle j ^ { l } _ { x, j } \right \rangle = i \sum _ { \sigma } \left ( e ^ { i A ^ { l } _ { \hat { x }, j } } \left \langle c ^ { \dagger } _ { l, j + \hat { x }, \sigma } c _ { l, j, \sigma } \right \rangle - e ^ { - i A ^ { l } _ { \hat { x }, j } } \left \langle c ^ { \dagger } _ { l, j, \sigma } c _ { l, j + \hat { x }, \sigma } \right \rangle \right ) \, \text{$4\pi\Phi/\Phi_{0}$} \, i$} \, \text{per flux} \, \text{$2$} \, \text{start}.$$

$$\langle J _ { z, j } \rangle = i t \sum _ { \sigma } \left ( e ^ { i A _ { z, j } } \left \langle c _ { 1, j, \sigma } ^ { \dagger } c _ { 2, j, \sigma } \right \rangle - e ^ { - i A _ { z, j } } \left \langle c _ { 2, j, \sigma } ^ { \dagger } c _ { 1, j, \sigma } \right \rangle \right ) \, \text{ using from } i \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \, \text{$i$} \right \}$$

and the average current density in the ˆ direction vany ishes. The current density satisfies the continuity equation ⟨∇ · J ⟩ = ⟨ ∂ ρ t ⟩ = i ⟨ H [ , ρ ] ⟩ = 0.

Phase diagram We first consider the T = 0 phase diagram and obtain for a given N the lowest energy state to determine the critical values of Φ / Φ 0 below which the BCS phase survives. We obtain the thermodynamic limit by extrapolation. (The number of vortices in the ground state varies as a function of N .) The resulting phase diagram is shown in Fig. 1(c) as a function of B and t .

To bring out the nature of the ground state, an important quantity is the relative phase, namely the gaugeinvariant difference between phases of the order parameter at two sites that lie directly across one another (i.e. have the same x -y coordinates) in the two layers. The Hamiltonian is invariant under the gauge transformation: c l,j,σ → ˜ c l,j,σ = e iϕ j,l c l,j,σ , ∆ j,l → ˜ ∆ j,l = e 2 iϕ j,l ∆ j,l , A l ˆ δ, j → ˜ A l ˆ δ, j = A l ˆ δ, j -ϕ j +ˆ δ,l + ϕ j ,l , and A z,j → ˜ A z, j = A z,j -ϕ j , + + ϕ j , -(the last two are needed to ensure the invariance of the hopping term). The gauge invariant relative phase is defined as arg(∆ + ,α / ∆ -,α ) + 2 A z,α with A z,α = 0 for our chosen gauge.

For sufficiently small Φ, the ground state is in the 'BCS phase.' Here the relative phase is uniformly zero (the phases in the two layers are locked by interlayer tunneling and the gap can be taken to be real); the gap amplitude is constant; there are no interlayer supercurrents; and there is no spontaneous net currents, although the top and the bottom layers have equal and opposite diamagnetic currents[19].

The BCS phase survives up to a critical magnetic flux. Fig. 2(a) depicts the relative phase for a magnetic flux just above the critical value. Evidently, there are regions where the relative phase remains nearly constant separated by regions where it rapidly changes by 2 π . Thus, rather than a relative phase linearly varying with x , the system finds it energetically favorable to create BCS stripes running parallel to the direction of the magnetic field. Associated with the formation of stripes is a periodic variation in the amplitude of the gap Fig. 2(d), and also the phase of the gap Fig. 2(g). The stripe phase also has a complex current pattern. Fig. 2(j) displays the current j x in the top layer; the current in the bottom layer is -j x and the interlayer current can be deduced from current conservation. The second column in Fig. 2 depicts these quantities for an intermediate flux and the third for a large flux. In the limit of large flux, we find the relative phase changes linearly with x , recovering the layer helical phase in Ref.[4], with a net change of about 4 π Φ Φ / 0 across the system corresponding to two vortices per flux quantum, and the gap amplitude is nearly constant.

Intuitively, the BCS stripe phase can be viewed as arising from a competition between the inter-layer tunneling and the in-plane magnetic flux. While the inter-layer tunneling tends to lock the SC order parameter phases between two layers, thus favoring the BCS phase, the in-plane magnetic flux tends to drive a variation of this relative phase by introducing inter-layer vortices to form the layer helical phase. The highly non-uniform BCS stripe phase emerges as a compromise between these two competing tendencies.

We have also evaluated the finite temperature phase diagram, shown in Fig. 1(d) for t = 0 5. . Fig. 1(e) shows how the gap amplitude varies as a function of T, where the gap | ∆ in the stripe phase denotes the spatial av-| erage of the gap amplitude. The transitions from either the BCS or the stripe phase to the normal phase follow the standard behavior. The average gap amplitude varies smoothly across the transition from the BCS to the stripe phase because the number of vortices per flux quantum rises continuously from zero. Fig. 1(f) shows the behavior of the relative phase near the BCS to stripe phase boundary as a function of T ; the phase becomes more linear as T approaches T , consistent with Ref.[4]. c

Finally, a remark is in order on the gauge invariant pair momentum given by p = ∇ ϕ -2 A , where ∆ = | ∆ | e iϕ . In our lattice model this can be written as p l x, j = ϕ l, j +ˆ x -ϕ l, j -2 A l ˆ x, j . Using this definition, we see that in the BCS phase, though there is no phase modulation of the gap parameter, there is an equal and opposite constant momentum in the top and bottom layers, which results in equal and opposite diamagnetic currents in the two layers. At high magnetic fields, there is equal and opposite contribution from both the phase and the gauge term in both layers. Since for high magnetic fields, ϕ = lq · x with constant q ( l = ± ), the gauge invariant momentum is spatially uniform and equal to zero. This is not the case, however, in the stripe phase (Fig. 2 (g) and (h)).

Discussion -The interplay between superconductivity and magnetic field in 2D systems has been of interest in numerous contexts, including the possibility of topological superconductivity supporting Majorana zero modes [20-33]. Here we predict a stripe phase in a bilayer superconductor exposed to a parallel magnetic field.

An earlier Landau-Ginzburg treatment of this problem by Liu [4] assumes the layer helical phase ansatz ∆ ( l r ) = ∆ | | e ilqx where l = ± is the layer index and q is spatially uniform. He finds q c = 2 eBz 0 2 ℏ in the limit of large B , where z 0 is the distance between two layers; the relative phase difference 2 q x c changes by 2 π in a distance ∆ x = 2 π/ q 2 c , which contains Bz 0 ∆ x = Φ 0 / 2 flux quanta, thus producing twice as many vortices as

the number of flux quanta. In contrast, the highly nonuniform stripe phase found in this work cannot be represented by layer helical phase with a single q in each layer. Qui and Zhou [19] have also studied this problem within the Landau-Ginzburg formalism, while assuming that the gap function has the periodicity of an MUC containing one flux quantum (the ratio of the number of vortices to the number of flux quanta can take only discrete values in this model); this periodicity results in a Bloch wave solution for the superconducting order parameter and a commensurate-incommensurate phase transition. The BCS stripe phase in our work does not necessarily belong to the Bloch wave solution as it in general does not respect the periodicity of the MUC.

By obtaining the self-consistent solutions of the BdG equations, we have shown that the application of a parallel B causes a transition into a stripe phase. We thus predict that the phase transition observed in 2H-stacked NbSe 2 in Ref. [34] is into a highly nonuniform stripe phase, which adiabatically evolves into the inter the helical phase of Ref.[4] for large B . This stripe phase may be identified most directly by measuring the oscillations in the amplitude of the gap, which are predicted to be strongest at low temperatures near the phase boundary.

Acknowledgments UN, JS (at Penn State) and JKJ were supported in part by the U. S. Department of Energy, Office of Basic Energy Sciences, under Grant no. DE-SC-0005042. Work by JS and ER was supported by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES) under Award DESC0022245. C.-X.L acknowledges the support from the NSF through The Pennsylvania State University Materials Research Science and Engineering Center [DMR2011839]. We acknowledge Advanced CyberInfrastructure computational resources provided by The Institute for CyberScience at The Pennsylvania State University.

## Appendix: Mapping bilayer TMD system with Ising spin-orbit coupling to our model

In this section, we will demonstrate that the model we used can serve as the low-energy theory of the bilayer transition-metal-dichalcogenide (TMD) system with Ising spin-orbit coupling. We start from the continuous model for a TMD monolayer. The single particle Hamiltonian for the layer l ( l = 1 2) at zero magnetic , field reads

$$\mathcal { H } _ { 0, l } ( p = \epsilon K + k ) = \xi _ { k } + ( - 1 ) ^ { l } \epsilon \beta _ { S O C } s _ { z }, \quad ( 6 )$$

where ξ k is the kinetic energy term, ϵ = ± is the valley index, β SOC is the Ising spin-orbit coupling strength, and s z is the Pauli matrix for the spin degree of freedom. In the second quantization notation with the electron annihilation operator denoted as c k ,ϵ,l,s ( s being the spin index), the Hamiltonian takes the form

$$\hat { H } _ { 0 } = \sum _ { k, \epsilon, s, s ^ { \prime } } c _ { k, \epsilon, l, s } ^ { \dagger } [ \mathcal { H } _ { 0, l } ( \epsilon K + k ) ] _ { s, s ^ { \prime } } c _ { k, \epsilon, l, s ^ { \prime } }. \quad ( 7 )$$

One can see that this Hamiltonian has a diagonal form with the spin-valley-layer locking, namely ⟨ + ↑| ˆ H 0 ,l | + ↑⟩ = ξ k + ( -1) l β SOC , ⟨ + ↓| ˆ H 0 ,l | + ↓⟩ = ξ k - -( 1) l β SOC , ⟨- ↑| ˆ H 0 ,l |- ↑⟩ = ξ k - -( 1) l β SOC , and ⟨ + ↑| ˆ H 0 ,l | + ↑⟩ = ξ k +( -1) l β SOC for non-zero Hamiltonian matrix elements. Since we only wish to consider the lowest energy bands and assume β SOC &gt; 0, we project the full Hamiltonian onto the basis A = {| l = 1 , ϵ = + , ↑⟩ , | l = 1 , ϵ = - ↓⟩ , , | l = 2 , ϵ = + , ↓⟩ , | l = 2 , ϵ = - ↑⟩} , . The low energy Hamiltonian then becomes:

$$\hat { H } _ { 0 } ^ { l o w } = \sum _ { k, \alpha \in A } ( \xi _ { k } - \beta _ { S O C } ) c _ { k, \alpha } ^ { \dagger } c _ { k, \alpha }.$$

Note that the spin-valley indices are opposite for the two layers in the above basis wavefunctions, and this will have substantial influence on the inter-layer tunneling.

The inter-layer tunneling term in the original basis (before projecting onto the lowest energy bands) preserves both spin and valley:

$$\hat { H } _ { t } = - t \sum _ { k, \epsilon, s } \left ( c ^ { \dagger } _ { k, \epsilon, 1, s } c _ { k, \epsilon, 2, s } + c ^ { \dagger } _ { k, \epsilon, 2, s } c _ { k, \epsilon, 1, s } \right ).$$

But, for l = 1 and l = 2, the spin and valley indices are opposite for low energy bands in the basis set A . Therefore the direct tunneling from this term in the low energy theory is zero. Inter-layer tunneling between low-energy bands can, however, be mediated by the combination of ˆ H t and the Zeeman coupling from the in-plane magnetic field.

We consider the in-plane magnetic field along the ˆ x direction. The Zeeman coupling is then given by

$$\hat { H } _ { Z } = g B _ { x } \sum _ { k, l, \epsilon, s } c _ { k, \epsilon, l, s } ^ { \dagger } ( \sigma _ { x } ) _ { s, s ^ { \prime } } c _ { k, \epsilon, l, s ^ { \prime } },$$

where g is the g-factor. From the Hamiltonian forms of ˆ H t and ˆ H Z , we see that while the inter-layer tunneling term preserves the spin-valley index while changing the layer index, the Zeeman term for an in-plane magnetic field changes the spin index while preserving the other two. Thus, we get an effective tunneling between the lowest energy bands from the second order term in perturbation theory:

$$\text{valley} _ { \text{h. and} } \quad \hat { H } _ { t } ^ { l o w } = \frac { t g B _ { x } } { \beta _ { S O C } } \sum _ { k } \left [ c _ { k, +, 1, \uparrow } ^ { \dagger } c _ { k, +, 2, \downarrow } + c _ { k, -, 1, \downarrow } ^ { \dagger } c _ { k, -, 2, \uparrow } + h. c. \right ].$$

The effective low energy theory is given by:

$$\hat { H } ^ { l o w } = \hat { H } _ { 0 } ^ { l o w } + \hat { H } _ { t } ^ { l o w }.$$

Due to the redundancy in the labelling of the basis wavefunctions, we can remove the valley index in the notation, so that | l = 1 , ϵ = + , ↑⟩ → | l = 1 , ↑⟩ , | l = 1 , ϵ = - ↓⟩ , → | l = 1 , ↓⟩ , | l = 2 , ϵ = + , ↓⟩ → | l = 2 , ↑⟩ , | l = 2 , ϵ = - ↑⟩ → | , l = 2 , ↓⟩ , and tgB x β SOC → t . The resulting Hamiltonian is a model for the spinful electrons in a bilayer system. We may include the orbital effect of magnetic field by performing the Peierls substitution and then implement the lattice regularization for ˆ H low to transform it to the lattice model, and with that, we can recover the noninteracting Hamiltonian H 0 in the main text. Thus, the in-plane Zeeman effect only provides a correction to the inter-layer tunneling parameter t .

- [1] M. P. Lilly, K. B. Cooper, J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Physical Review Letters 82 , 394 (1999).
- [2] M. M. Fogler, A. A. Koulakov, and B. I. Shklovskii, Physical Review B 54 , 1853 (1996).
- [3] V. J. Emery, S. A. Kivelson, and J. M. Tranquada, Proceedings of the National Academy of Sciences 96 , 8814 (1999).
- [4] C.-X. Liu, Physical Review Letters 118 , 087001 (2017).
- [5] P. Fulde and R. A. Ferrell, Physical Review 135 , A550 (1964).
- [6] A. Larkin and Y. N. Ovchinnikov, Soviet Physics-JETP 20 , 762 (1965).
- [7] D. Agterberg, Non-Centrosymmetric Superconductors: Introduction and Overview , 155 (2012).
- [8] J. J. Kinnunen, J. E. Baarsma, J.-P. Martikainen, and P. T¨ orm¨, a Reports on Progress in Physics 81 , 046401 (2018).
- [9] A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and J. L. Sarrao, Physical Review Letters 91 , 187004 (2003).
- [10] Y. Matsuda and H. Shimahara, Journal of the Physical Society of Japan 76 , 051005 (2007).
- [11] G. B. Partridge, W. Li, R. I. Kamar, Y. an Liao, and R. G. Hulet, Science 311 , 503 (2006).
- [12] M. M. Parish, S. K. Baur, E. J. Mueller, and D. A. Huse, Physical Review Letters 99 , 250403 (2007).
- [13] J. J. Kinnunen, J. E. Baarsma, J.-P. Martikainen, and
- P. T¨ orm¨, a Reports on Progress in Physics 81 , 046401 (2018).
- [14] D. Zhao, L. Debbeler, M. K¨hne, S. Fecher, N. Gross, u and J. Smet, Nature Physics 19 , 1599 (2023).
- [15] P. Wan, O. Zheliuk, N. F. Q. Yuan, X. Peng, L. Zhang, M. Liang, U. Zeitler, S. Wiedmann, N. E. Hussey, T. T. M. Palstra, and J. Ye, Nature 619 , 46 (2023).
- [16] M. Croitoru and A. Buzdin, Condensed Matter 2 , 30 (2017).
- [17] J. Schirmer, J. K. Jain, and C. X. Liu, Physical Review B 109 , 134518 (2024).
- [18] H. Watanabe, Journal of Statistical Physics 177 , 717 (2019).
- [19] G.-W. Qiu and Y. Zhou, Physical Review B 105 , L100506 (2022).
- [20] M. Stone and S.-B. Chung, Physical Review B 73 , 014505 (2006).
- [21] M. Oshikawa, Y. B. Kim, K. Shtengel, C. Nayak, and S. Tewari, Annals of Physics 322 , 1477 (2007).
- [22] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Reviews of Modern Physics 80 , 1083 (2008).
- [23] J. D. Sau, R. M. Lutchyn, S. Tewari, and S. D. Sarma, Physical Review Letters 104 , 040502 (2010).
- [24] T. D. Stanescu, J. D. Sau, R. M. Lutchyn, and S. D. Sarma, Physical Review B 81 , 241310 (2010).
- [25] J. Alicea, Reports on Progress in Physics 75 , 076501 (2012).
- [26] H.-H. Hung, P. Ghaemi, T. L. Hughes, and M. J. Gilbert, Physical Review B 87 , 035401 (2013).
- [27] C. Beenakker, Annual Review of Condensed Matter Physics 4 , 113 (2013).
- [28] J. Zhou, Y.-J. Wu, R.-W. Li, J. He, and S.-P. Kou, Europhysics Letters 102 , 47005 (2013).
- [29] R. R. Biswas, Physical Review Letters 111 , 136401 (2013).
- [30] S. D. Sarma, M. Freedman, and C. Nayak, npj Quantum Information 1 , 1 (2015).
- [31] T. Liu and M. Franz, Physical Review B 92 , 134519 (2015).
- [32] J. M. Murray and O. Vafek, Physical Review B 92 , 134520 (2015).
- [33] E. D. Smith, K. Tanaka, and Y. Nagai, Physical Review B 94 , 064515 (2016).
- [34] P. Wan, O. Zheliuk, N. F. Yuan, X. Peng, L. Zhang, M. Liang, U. Zeitler, S. Wiedmann, N. E. Hussey, T. T. Palstra, et al. , Nature , 1 (2023).
Uddalok Nag Jonathan Schirmer Enrico Rossi C. -X. Liu J. K. Jain
BCS Stripe Phase in Coupled Bilayer Superconductors
As a signature of competing correlations, stripes occur in a variety of
strongly correlated systems, such as high temperature superconductors (SCs) and
quantum Hall effect. We study a double layer SC in the presence of a parallel
magnetic field $B$ within the Bogoliubov-de Gennes framework. We find that for
low $B$ the system remains in the ``Bardeen-Cooper-Schrieffer (BCS) phase" with
a spatially uniform gap, but with increasing $B$, a transition occurs into a
phase which contains stripes of the BCS phase separated by regions where the
interlayer phase difference rotates by $2\pi$ due to the presence of
inter-layer vortices. This stripe phase is predicted to manifest through
oscillations in the amplitude of the SC gap and an alternating pattern of
supercurrents. We will comment on the relation to previous works based on the
Landau-Ginzburg formalism as well as on the possible experimental realization
and signature of this phase.