Cascade of vestigial orders in two-component superconductors: Nematic, ferromagnetic, s -wave charge- 4e, and d -wave charge- 4e states

Electronically ordered states that break multiple symmetries can melt in multiple stages, similarly to liquid crystals. In a partially melted phase, known as vestigial phase, a bilinear made out of combinations of the multiple components of the primary order parameter condenses. Multicomponent superconductors are thus natural candidates for vestigial order since they break both the U(1)-gauge and also time-reversal or lattice symmetries. Here, we use group theory to classify all possible real-valued and complex-valued bilinears of a generic two-component superconductor on a tetragonal or hexagonal lattice. While the more widely investigated real-valued bilinears correspond to vestigial nematic or ferromagnetic order, the little explored complex-valued bilinears correspond to a vestigial charge-4e condensate, which itself can have an underlying s-wave, dx2-y2-wave, or dxy-wave symmetry. To properly describe the fluctuating regime of the superconducting Ginzburg-Landau action and thus access these competing vestigial phases, we employ both a large-N and a variational method. We show that while vestigial order can be understood as a weak-coupling effect in the large-N approach, it is akin to a moderate-coupling effect in the variational method. Despite these distinctions, both methods yield similar results in wide regions of the parameter space spanned by the quartic Landau coefficients. Specifically, we find that the nematic and ferromagnetic phases are the leading vestigial instabilities, whereas the various types of charge-4e order are attractive albeit subleading vestigial channels. The only exception is for the hexagonal case, in which the nematic and s-wave charge-4e vestigial states are degenerate. We discuss the limitations of our approach, as well as the implications of our results for the realization of exotic charge-4e states in material candidates.