Orbitally limited pair-density-wave phase of multilayer superconductors

We investigate the magnetic field dependence of an ideal superconducting vortex lattice in the parity-mixed pair-density-wave phase of multilayer superconductors within a circular cell Ginzburg-Landau approach. In multilayer systems, due to local inversion symmetry breaking, a Rashba spin-orbit coupling is induced at the outer layers. This combined with a perpendicular paramagnetic (Pauli) limiting magnetic field stabilizes a staggered layer dependent pair-density-wave phase in the superconducting singlet channel. The high-field pair-density-wave phase is separated from the low-field BCS phase by a first-order phase transition. The motivating guiding question in this paper is: What is the minimal necessary Maki parameter αM for the appearance of the pair-density-wave phase of a superconducting trilayer system? To address this problem we generalize the circular cell method for the regular flux-line lattice of a type-II superconductor to include paramagnetic depairing effects. Then, we apply the model to the trilayer system, where each of the layers are characterized by Ginzburg-Landau parameter κ₀ and a Maki parameter αM. We find that when the spin-orbit Rashba interaction compares to the superconducting condensation energy, the orbitally limited pair-density-wave phase stabilizes for Maki parameters αM>10.