Superconductivity mediated by Cooper quartets of charge 4e is a phenomenon of key importance for both the understanding of certain exotic superconductors and the engineering of quantum memories protected at the hardware level by topological order. Josephson junction arrays in the shape of the dice lattice constitute one of the main candidates for the realization of this phase of matter. Here, we analyze numerical signatures of the emergence of this exotic phase when superconducting dice arrays are frustrated by inserting one third of a flux quantum per rhombic plaquette. We adopt simulations of relaxation dynamics to study the critical current of such devices and analyse the Fourier decomposition of their bulk supercurrents. A further characterization of these systems at finite temperature is obtained through Monte Carlo techniques and two-dimensional infinite tensor networks. Our results indicate that the observed peaks of the critical current and temperature at frustration 1/3 correspond to a superconductor-insulator phase transition compatible with the deconfinement of half-vortices, while the low-temperature correlation functions of the model confirm the onset of a 4e superconducting phase mediated by Cooper quartets. We finally address the effects of Josephson energy and flux disorder typical of experimental arrays, and comment on the role of charging energies in the corresponding two-dimensional quantum model.
Superconductivity mediated by Cooper quartets of charge 4e is a phenomenon of key importance for both the understanding of certain exotic superconductors and the engineering of quantum memories protected at the hardware level by topological order. Josephson junction arrays in the shape of the dice lattice constitute one of the main candidates for the realization of this phase of matter. Here, we analyze numerical signatures of the emergence of this exotic phase when superconducting dice arrays are frustrated by inserting one third of a flux quantum per rhombic plaquette. We adopt simulations of relaxation dynamics to study the critical current of such devices and analyse the Fourier decomposition of their bulk supercurrents. A further characterization of these systems at finite temperature is obtained through Monte Carlo techniques and two-dimensional infinite tensor networks. Our results indicate that the observed peaks of the critical current and temperature at frustration 1/3 correspond to a superconductor-insulator phase transition compatible with the deconfinement of half-vortices, while the low-temperature correlation functions of the model confirm the onset of a 4e superconducting phase mediated by Cooper quartets. We finally address the effects of Josephson energy and flux disorder typical of experimental arrays, and comment on the role of charging energies in the corresponding two-dimensional quantum model.
We study the quantum phase transition between the insulating and the globally coherent superfluid phases in the Bose-Hubbard model with T3 structure, the ``dice lattice.'' Even in the absence of any frustration the superfluid phase is characterized by modulation of the order parameter on the different sublattices of the T3 structure. The zero-temperature critical point as a function of magnetic field shows the characteristic ``butterfly'' form. At full frustration the superfluid region is strongly suppressed. In addition, due to the existence of the Aharonov-Bohm cages at f=1∕2, we find some evidence for the existence of an intermediate insulating phase characterized by a zero superfluid stiffness but finite compressibility. In this intermediate phase bosons are localized due to the external frustration and the topology of the T3 lattice. We name this new phase the Aharonov-Bohm insulator. In the presence of charge frustration the phase diagram acquires the typical lobe structure. The form and hierarchy of the Mott insulating states with fractional fillings are dictated by the particular topology of the T3 lattice. The results presented were obtained by a variety of analytical methods: mean-field and variational techniques to approach the phase boundary from the superconducting side and a strongly coupled expansion appropriate for the Mott insulating region. In addition we performed quantum Monte Carlo simulations of the corresponding (2+1)-dimensional XY model to corroborate the analytical calculations with a more accurate quantitative analysis. We finally discuss experimental realization of the T3 lattice both with optical lattices and with Josephson junction arrays.
We study the differential resistivity transition of two-dimensional superconducting arrays induced by an external driving current, in the presence of thermal fluctuations and a magnetic field corresponding to f flux quantum per plaquette. Recent experiments have identified this transition as a dynamic vortex Mott insulator transition at vortex densities near rational values of f. The critical behavior is determined from a scaling analysis of the current-voltage relation near the transition, obtained by Monte Carlo simulations of a Josephson-junction array model in the vortex representation. For a square-lattice array, the critical exponents obtained near f=1/2 are consistent with the experimental observations. The same scaling behavior is observed near f=1/3. For a honeycomb array, although similar results are obtained for f=1/3, the transition is absent for f=1/2, consistent with an incommensurate vortex phase.
Quantum design of Cooper quartets in a double quantum dot system coupled to ordinary superconducting leads is presented as a novel platform for the study of an elusive many-body state of matter, that is at the basis of the phenomenon of charge-4e superconductivity. A fundamentally novel, maximally correlated ground state, in the form of a superposition of vacuum |0⟩ and four-electron state |4e⟩, emerges as a narrow resonance and it is promoted by an attractive interdot interaction. A novel phenomenology in the dissipationless transport regime is elucidated, that yields typical flux quantization in units of h/4e and manifests in non-local multi-terminal coherence and in two-Cooper pair transport properties mediated by the quartet ground state. The results open the way to the exploration of correlation effects and non-local coherence in hybrid superconducting devices, parity-protected quantum computing schemes and more generally, the work poses the basis for the design and simulation of novel correlated states of matter starting from ordinary ingredients available in a quantum solid state lab.
Unlike conventional charge-2e superconductors, a charge-4e superconductor exhibits long-range coherence of electron quartets rather than Cooper pairs. Clear zero-temperature realizations of charge-4e superconductivity remain rare. Here, we investigate the zero-temperature phase diagram of the attractive SU(4) Hubbard model with numerically exact, large-scale quantum Monte Carlo (QMC) simulations overcoming major technical hurdles. We identify both charge-2e and charge-4e superconducting phases. Upon increasing interaction, charge-2e correlations are suppressed and eventually vanish, while the charge-4e correlations remain robust and converge with system size, signaling the onset of a quartet-condensed phase. Interestingly, across the charge-2e--charge-4e transition, single electrons remain gapped, while charge-2e correlations exhibit a scaling behavior inconsistent with a conventional Landau description. These features are naturally captured by a fractionalized framework in which the physical charge-2e order parameter is a composite field coupled to an emergent non-Abelian gauge structure. We formulate an Sp(4) gauge-Higgs theory that realizes deconfined quantum pseudocriticality between the Higgs (charge-2e) phase and the confined (charge-4e) phase. The Sp(4) gauge-Higgs theory yields pseudocriticality through a fixed-point collision, and its one-loop collision-point exponents quantitatively track the QMC results. Our results establish charge-4e superconductivity as a bona fide zero-temperature phase, provide a simple model for future studies in a numerically exact framework, and reveal an unconventional route to superconducting criticality.
Superconductivity provides a canonical example of a quantum phase of matter. When superconducting islands are connected by Josephson junctions in a lattice, the low-temperature state of the system can map to the celebrated XY model and its associated universality classes. This has been used to experimentally implement realizations of Mott insulator and Berezinskii-Kosterlitz-Thouless transitions to vortex dynamics analogous to those in type-II superconductors. When an external magnetic field is added, the effective spins of the XY model become frustrated, leading to the formation of topological defects (vortices). Here, we observe the many-body dynamics of such an array, including frustration, via its coupling to a superconducting microwave cavity. We take the design of the transmon qubit, but replace the single junction between two antenna pads with the complete array. This allows us to probe the system at 10 mK with minimal self-heating by using weak coherent states at the single (microwave) photon level to probe the resonance frequency of the cavity. We observe signatures of an ordered vortex lattice at rational flux fillings of the array.
A one-dimensional quantum version of the frustrated XY (planar rotor) model is considered which can be physically realized as a ladder of Josephson junctions at half a flux quantum per plaquette. This system undergoes a superconductor to insulator transition at zero temperature as a function of charging energy. The critical behavior is studied using a Monte Carlo transfer matrix applied to the path-integral representation of the model and a finite-size-scaling analysis of data on small system sizes. Depending on the ratio between the interchain and intrachain couplings the system can have single or double transitions which is consistent with the prediction that its critical behavior should be described by the two-dimensional classical XY-Ising model.