Electrically controllable superconducting memory effect in UTe<sub>2</sub>

Nature正文已收录本站

Main

If a computer could be assembled from superconducting components, the energy efficiency would far surpass that of conventional electronics. Historic research efforts towards this goal yielded pivotal breakthroughs in the development and discovery of scanning tunnelling microscopy15 and high-temperature superconductivity16. Although recent strides have been taken in advancing superconducting rectification17,18 and switching19 technologies, realizing read/writable memory functionality in superconducting platforms has remained challenging.

Encodable superconducting memory functionality has been demonstrated in ferromagnet–superconductor heterostructures12,13, as well as through history-dependent trapping of magnetic flux20,21. In some type II superconductors such as NbSe2, a hysteretic modulation of the critical current density Jc has been observed near a so-called ‘peak-effect’22 region in which Jc suddenly increases on approaching the upper critical field. Such behaviour is understood to be governed by a dynamic competition between the injection of a highly disordered, strongly pinned vortex phase through surface edge barriers and the subsequent annealing of this disorder by a bulk transport current23. By tuning the amplitude, frequency and direction of the driving current, the spatial footprint of this disordered vortex matter may be programmably manipulated24. The material effectively archives its electrical history within this spatial distribution, giving the vortex landscape the ability to effect non-volatile information storage25.

Here we study vortex dynamics within the multiphase spin-triplet superconductor candidate UTe2 (refs. 6,7,8). In a specific portion of the complex phase landscape, we observe hysteretically tunable Jc properties analogous to those of NbSe2 (ref. 23). However, rather than becoming pronounced near to the critical temperature, as per the peak effect, in UTe2, this ‘memory effect’ manifests at very low temperatures over a magnetic field interval in which this material is known to transition between two distinct superconducting phases9,10. We therefore posit that, while in NbSe2 thermally induced competition between elastic and pinning energies yields metastable memory phenomenology23, by contrast, in UTe2, it seems that competing interactions between two distinct vortex structures—each native to their separate superconducting phases—drives the observed superconducting memory effect.

Multiphase superconductivity

UTe2 is an orthorhombic heavy-fermion superconductor, which, under various conditions of applied pressure and magnetic field, has been reported to have up to six different superconducting phases8,9,26,27,28,29,30,31,32,33,34,35,36,37. The most easily accessible transition between two superconducting states is found at ambient pressure for a magnetic field B aligned along the hard \(\hat{b}\text{-axis}\). Bulk-sensitive specific heat measurements9 have discerned a thermodynamic phase boundary between two superconducting states at around 15–20 T. These superconductive phases are typically referred to as SC1 (the low-B state) and SC2 (at higher B; Fig. 1a; ref. 31). There is a substantial, growing body of evidence that SC1 and SC2 have distinct order parameters8,9,10,34,37,38,39.

Fig. 1: Electrical switching of an intrinsic superconducting memory effect.

a, The low-temperature phase diagram of UTe2 for magnetic field B applied along the hard crystallographic \(\hat{b}\) direction31. The region of phase space in which memory effects are manifested is coloured purple, lying between two distinct superconducting phases (SC1 and SC2). b, The measured voltage V across a UTe2 single crystal as a function of dc current density J. Panels are arranged chronologically from left to right, with the insets (and colour coding) showing how J was modulated as a function of time. When J is ramped smoothly, the profile of V(J) is in its equilibrium state and no hysteresis is observed. By contrast, after a discontinuous change in J from a large magnitude abruptly to zero (labelled ‘Pump’), hysteresis in V(J) is observed, plotted in red and orange. Smoothly sweeping the current from a large negative value to zero (‘Erase’) then resets the system to the original V(J) profile with lower Jc. All data were acquired at 50 mK with \(B\Vert \hat{b}\). c, Data collected on sample S2 at 7 T and 50 mK. Here we plot just two traces, to highlight the hysteretic out-of-equilibrium form of V(J) in the memory state. d, The lower part shows V measured as a function of time for modulations of J depicted in the upper part, cycling through a sequence of perturbation, measurement and erasure (resetting) protocols, as described in e. By sitting at this point in the J–V curve (sourcing J = 2 Acm−2), successively switching in and out of the memory state yields a voltage versus time profile similar to supercurrent rectification by the superconducting diode effect17.

Source data

Recent measurements of the surface properties of the SC1 state at low B, made by scanning tunnelling microscopy and scanning SQUID, have revealed several anomalous vortex features. These include a mirror-asymmetric profile of the vortex cores40, which form in doublets41,42, with a tendency to organize into extended stripes of vortices43. These vortex structures form for B > Bc1 ≈ 4 mT (ref. 44) and, so far, have been found to persist up to at least 8 T (ref. 41).

The vortex behaviour at higher B is also unusual. Current-dependent magnetotransport measurements have discerned a region of anomalously low Jc between the SC1 and SC2 states at around 15–20 T (ref. 45). Alongside signatures in the magnetic susceptibility, these observations have been taken to indicate a region of coexistence between the SC1 and SC2 states10 that, hereafter, we refer to as the SC1.5 regime.

We investigated the current–voltage (J–V) characteristics of UTe2 in this region of phase space in which SC1 and SC2 show signs of coexistence. When the material is entirely in either the SC1 or the SC2 state, we observe forms of V(J) that are typical for a type II superconductor in a magnetic field. By contrast, in the intermediate SC1.5 region, we find that abrupt perturbative modulations of J lead to unusual hysteretic features in the form of V(J), increasing the value of Jc. The system then stays in this new hysteretic high-Jc state on subsequent measurements, only returning to the original V(J) profile after the application of a sufficiently large perturbation that acts to reset the system. This hysteretic tuning seems robustly long-lived (persisting over a timescale of at least several hours), constituting a previously unknown electrically controllable superconducting memory effect.

Superconducting memory effect in UTe2

We plot a chronological sequence of J–V curves in Fig. 1b. Throughout this study, all excitations were applied as direct currents (Methods), unless otherwise specified. With UTe2 in its intermediate SC1.5 state at low temperatures, when the current is swept smoothly, V(J) is a single-valued function, with no hysteresis observed for increasing or decreasing current ramps. By contrast, after suddenly decreasing J from a high value down to zero, a larger value of Jc is then observed on smoothly sweeping J back up again (red curve, top-right panel of Fig. 1b). This hysteretic new form of V(J) is also manifested for opposite polarity current ramps (orange curve, subsequent panel). Then, on smoothly sweeping down from a large magnitude of current back to zero, the initial state is recovered, with the original form of V(J) retraced for subsequent rising and falling current sweeps (blue curves). We identify this hysteretic tuning of Jc as a current-controlled memory effect and refer to the hysteretically higher Jc value as the characterization of the system being in the ‘memory state’.

In Fig. 1c, we show data measured on a second sample (labelled S2). For clarity, here we show just the non-hysteretic V(J) curve (in blue, labelled ‘Equilibrium’) and the post-perturbation memory state curve (in red, labelled ‘Out-of-equilibrium’). On successive measurements, we find that the hysteretic form of V(J) in the memory state persists over a timescale of at least several hours. Therefore, although the memory state is away from the equilibrium state of the system, it does not seem to be transiently metastable in character, giving promise for possible future development of non-volatile memory functionalities.

The tunability of the memory state enables us to perform pump/probe operations to switch between different states (Fig. 1d,e). For example, if we choose a J value that lies at the bottom of the hysteresis loop, then perturbing the system to switch from the equilibrium (non-memory) state to the out-of-equilibrium (memory) state changes V between non-zero and zero values by raising Jc (blue and red data points in Fig. 1d). These high and low states could, for instance, be labelled 1 and 0 for computational memory purposes.

We find that the memory effect can be controlled by tuning the system in three different ways: by modulating the amplitude, duration and relaxation rate of the perturbative excitation current Jex (Fig. 2). For each of these tuning parameters, we find that larger hysteresis loops are induced for stronger perturbations, up to some saturated maximal values (Extended Data Fig. 1). Furthermore, each of these tuning axes yields a similar evolution in temperature of the memory effect. Figure 3 depicts the temperature T dependence over the interval 50 mK ≤ T ≤ 1.0 K for amplitude tuning. At low temperatures below 0.4 K, bright yellow regions in Fig. 3b identify a strong memory effect characterized by the observation of a large ΔJ for high Jex. However, at higher temperatures, this quickly diminishes, with hysteresis loops closing for T > 0.6 K. Similar behaviour is also observed under duration and relaxation tuning, as shown in Extended Data Fig. 2. Such a low temperature scale of the memory effect in UTe2 is surprising, given that the superconducting Tc = 2.1 K in zero field and is still >1 K in this field range. This points to some intrinsic low-energy phenomenon, distinct from the superconductivity itself, being responsible for the observed hysteretic V(J) phenomena.

Fig. 2: Three tuning parameters to control the memory effect.

a–c, Hysteretic response profiles for modulations of an applied dc excitation’s amplitude (a), duration (b) and relaxation rate (c). In each case, the greater the perturbation to the system—be it a larger amplitude Jex, longer duration Δt or more rapid relaxation rate α—the greater the resulting hysteresis loop in V(J), up to some saturation value (see Extended Data Fig. 1 for further analysis). For the heatmaps, ΔJ is calculated as the change in J at fixed values of V. The regions of highest ΔJ, coloured yellow, correspond to the strongest memory effect.

Source data

Fig. 3: Temperature dependence of the memory effect.

a, V(J) at incremental temperatures as indicated. The inset shows the sequence of current modulations. Although large hysteresis loops are recorded at low temperature, these have closed by 650 mK. At 1.0 K, the curvature of V(J) is markedly different compared with the equilibrium curve at 50 mK (see also Supplementary Fig. 7). b, Heatmaps of ΔJ at incremental temperatures for modulations of the excitation amplitude Jex up to 20 A cm−2 measured on sample S2 at 14 T. Again, yellow colouring indicates the strongest memory effect, which sharply diminishes at higher temperatures.

Source data

Mapping the memory effect

To reveal the full domain of phase space in which memory effects are present, we performed a detailed mapping of the hysteretic V(J) characteristics of UTe2 as a function of magnetic field strength and rotation angle up to 30 T at 50 mK, which is summarized in Fig. 4. At low B, the form of V(J) is single-valued and does not exhibit hysteresis under perturbative amplitude, timescale or relaxation rate tuning. As B is increased, hysteretic signatures of the memory effect start to appear at around 7 T, which persist up to higher B until the phase boundary out of SC1 into SC2 is crossed. The exact value of B at which this occurs varies slightly from sample to sample, depending on crystalline quality9,31—for sample S1, this is just above 20 T, whereas for S2, it is just below 20 T (Fig. 4e and Extended Data Figs. 3 and 4). At high B around 30 T, the material is then fully in the SC2 phase and the memory effect is no longer exhibited (Fig. 4c).

Fig. 4: Mapping the memory region between SC1 and SC2.

a, Schematic phase landscape of UTe2 for rotations of B by an angle θ from \(\hat{b}\) towards \(\hat{c}\). The memory region (coloured purple)—identified by non-zero ΔJ at low T—is located at the intersection of the SC1 and SC2 domains. We refer to this as SC1.5. b, The equilibrium profile of V(B) at 50 mK for successive J values, as indicated by the colour scale. An anomalous non-monotonic kink in V(B) at high J is observed at the boundary of the memory region. c, Sequential modulations of J at different fields for \(B\Vert \hat{b}\) (0°). d, Modulations of J at different magnetic field tilt angles θ at B = 15 T. The inset defines the measurement protocol. e,f, Evolution of ΔJ as a function of B at 0° and 15° (e) and as a function of θ at 15 T and 29 T (f). g, Equilibrium V(θ) for J = 2.5 A cm−2, B = 15 T, T = 50 mK measured on sample S2. The transition from the memory region into SC1 is marked by a sudden jump in V at around θ = 20°, indicating a change in the vortex properties. This angle is where, in panel f, ΔJ reaches zero within the resolution of the measurement. h, Magnetic field sweeps of the effective resistivity ρ for incremental J values as indicated. The data in this panel were acquired by low-frequency ac measurements, whereas all other data in this article are from dc measurements (Methods). Although zero resistance is observed over all B in the low J limit, at higher J, a complex profile of flux-flow behaviour is exhibited.

Source data

Tilting B provides insight into the character of the memory effect. The SC2 phase of UTe2 only resides in a narrow range of θ close to \(B\parallel \hat{b}\), which, for high-quality salt-flux-grown crystals such as those studied here, extends out to θ ≈ 20° (ref. 31). We find that the memory region closely tracks this. In Fig. 4e,f, we plot the value of ΔJ observed for a perturbative Jex for various B and θ (see Extended Data Fig. 5 for raw V(J) curves). At 0°, no change in ΔJ is observed at low B, which quickly rises for B > 7 T and peaks just above 0.2 A cm−2 at around 10 T. A very similar profile is exhibited at θ = 15°, indicating that the strength of the memory effect is not strongly peaked at the \(\hat{b}\text{-axis}\) but is similar in magnitude throughout the B–θ parameter space directly below SC2. Then, on rotating beyond the range of SC2, ΔJ abruptly diminishes back to zero (Fig. 4f). This truncation of the memory region is also reflected in the variation of V measured in equilibrium for a fixed J (Fig. 4g), which suddenly jumps on departing the memory region into exclusively SC1 close to θ = 20°. The close connection between the angular range in which SC2 sits above SC1, and the observed domain of phase space that hosts the memory effect, underscores that the interplay of these two distinct superconducting states is imperative in generating the observed hysteretic V(J) phenomena.

Discussion

Now we turn to considering the likely microscopic origin of the memory effect in UTe2. First, given that we are applying substantial electrical current pulses at dilution fridge temperatures, we should consider any extrinsic heating effects that might be present. In the Supplementary Information, we plot the measured temperature as a function of time during the current pulse measurement protocols and present further discussion as to why we do not believe that Joule heating effects greatly affect our measurements. One strong argument against Joule heating artefacts comes from the sensitivity of the memory effect to the relaxation rate of the dc stimulus. In Fig. 2, we showed that stronger memory effect correlates with larger amplitude and duration tuning. For both of these, the larger stimulus correlates with more energy deposited into the system (through lead/contact resistances). However, the memory effect is also strongly correlated with faster relaxation rates α, which anticorrelates with the total energy transferred into the system. If the memory effect were because of some relatively straightforward Bean–Livingston barrier or Bean model pinning effect process46, this should be insensitive to how the current is turned off, depending only on the total energy imparted. We note, however, that it is therefore hard to rule out the possibility of thermal quenching playing a role. We posit that the acute sensitivity of the memory effect to the turn-off speed implies some dynamic reorganization or annealing of vortex domains, probably because of some interesting intrinsic physical property (or properties) of UTe2.

Further arguments against extrinsic heating effects originate from the fact that this phenomenon is only observed in the narrow regime of B–θ space directly underneath SC2, with no hysteretic V(J) phenomena observed in purely SC1 or SC2. The memory effect is most pronounced at the lowest temperatures and quickly diminishes on warming (Fig. 3), providing perhaps the strongest argument against any heating-induced artefacts. This is in sharp contrast to the case of NbSe2, in which hysteretic V(J) profiles are enabled close to Tc (ref. 23). Instead, our observations strongly indicate that the hysteretic V(J) behaviour of UTe2 is because of some intrinsic low-energy phenomenon, with an energy scale much lower than Tc.

We propose that the memory effect of UTe2 may be the result of competition between two different vortex structures—one native to SC1 and the other to SC2. In the equilibrium state, the vortices arrange themselves in a highly ordered configuration to minimize entropy and energetic costs. On perturbing the system, current densities above Jc are applied, inducing a flux-flow regime in which vortices are dissipatively travelling through the sample. If the perturbation is sufficiently sharp—for example, if the relaxation rate α is very fast—then when J drops back below Jc, vortex flow suddenly halts and the vortices are abruptly locked into their new positions. We can imagine a glassy vortex state has now been quenched in the material, which has much higher disorder than the original equilibrium vortex lattice. Higher disorder implies stronger pinning forces between vortices and hence the higher Jc values that characterize the memory effect. The system can then be reset by smoothly raising J back above Jc before gradually ramping back down to zero current. This could be thought of as annealing the system, thereby returning to the equilibrium low-Jc state. In the Supplementary Information, we provide phenomenological modelling of the memory effect and show how our experimental observations can be qualitatively understood within this interpretation.

Outlook

An attractive technological use case for the superconducting memory effect is as an ultralow-power alternative to solid-state or hard disk drive memory functionality. Moreover, the inherent plasticity of the effect presents a promising route towards cryogenic neuromorphic processing. For the bulk millimetre-sized single-crystal specimens investigated here, to probe (read) whether the material is or isn’t in the memory state—which could correspond to a binary 1 or 0—required a power consumption of about 1 nW. Scaling down even just mesoscopically to the order of one square micron in size would reduce this to ≲1 pW—orders of magnitude less than for nanometre-sized silicon elements47. Furthermore, whereas cryotron-based48 superconducting devices require transitioning to the normal state, UTe2 has the distinct advantage that both read and write operations can be performed while superconducting. Recent studies have demonstrated the ability of this material to be carved at the micron-scale by focused ion beams32,49, making the prospect of lithographically processing UTe2 wafers50 a realistic, albeit ambitious, goal.

Further work is required to assess whether similar memory phenomena could be exhibited by other multiphase superconductors such as CeRh2As2 (ref. 1) or UPt3 (ref. 2). We hope that our discovery here may help catalyse efforts to discover new multiphase superconductors composed of more abundant elements, which might exhibit the memory effect to higher temperatures. Although the list of known multiphase superconductors is short, notable recent additions from materials as diverse as moiré graphene51 and high-Tc nickelates52 suggest that further promising candidates await discovery.

In summary, we studied the J–V characteristics of the multiphase superconductor UTe2. In an intermediate phase space regime bordering two superconducting phases, we observed anomalous hysteresis in the value of the critical current, dependent on the history of applied direct current ramps and pulses. We showed that modulating the amplitudes and timescales of dc electrical stimuli can selectively tune between low and high critical current states, with the material then ‘remembering’ which state it is in. We interpret this to be because of different vortex configurations being either abruptly quenched or smoothly annealed. Our work demonstrates that, in principle, the vortex matter of an unconventional superconductor can be electrically manipulated and used to encode and process information, thereby opening up new possibilities for low-energy electronic devices, cryogenic neuromorphic memory and quantum computational hardware.

Methods

Sample preparation

Single-crystal UTe2 specimens were grown in a flux of molten salts53 by the recipe detailed in ref. 54. Samples were contacted by spot-welding four gold wires, each of 50 μm in diameter, for respective current and voltage leads. These were subsequently connected to copper twisted pairs, which were soldered to the wiring of the probe. Each of the three samples investigated in this study had J sourced along the [100] direction, with leads affixed to the (001) surface. Residual resistivity ratio measurements were performed using a Quantum Design Physical Properties Measurement System in Cambridge with low-frequency ac excitations. The residual resistivity ratio values for the three samples are S1:105, S2:44 and S3:395.

Electrical transport measurements

The hysteretic J–V properties of UTe2 were investigated by performing direct current measurements using a Keithley 6221/2182A combined current source–multimeter system. The in-built pulse-delta sweep mode and arbitrary waveform generator features were used. Measurements were performed in a 30 T superconducting magnet with Oxford Instruments dilution fridge system at the Synergetic Extreme Condition User Facility (SECUF) in Beijing55. All measurements presented in this study were performed in that system, except those in Fig. 4h and Supplementary Fig. 2 that were performed at the National High Magnetic Field Laboratory in Florida, in a resistive magnet using a 3He system, in which low-frequency ac measurements were obtained using a Keithley 6221 current source and a Stanford Research 860-series lock-in amplifier. For the electrical switching study at SECUF, a sequence of pulse-delta measurements with different applied currents was defined by entering the starting current value, ending current value and the number of steps for the sweep. To control the amplitude, duration and relaxation rate, a varying shape excitation pulse was defined by 100 points and generated with the arbitrary waveform feature of the source–multimeter system. Control of the instrument was realized by custom-developed Python code with the pymeasure package56, which is included as a supplement to this article (Supplementary Information). The reset protocol before performing a perturbative excitation was achieved by sweeping the current from 40 mA to 0 mA smoothly with a 2-mA step over a 20-s time interval with the pulse-delta sweep mode. The cycling sequence in Fig. 1d was performed by the pulse-delta mode. The time sequence is defined by entering the current amplitude and number of steps. Further electrical transport measurements are presented in Extended Data Figs. 6, 7, 8 and 9.

Data availability

The datasets supporting the findings of this study are available from the University of Cambridge Apollo Repository57. Source data are provided with this paper.

References

  1. Khim, S. et al. Field-induced transition within the superconducting state of CeRh2As2. Science 373, 1012–1016 https://doi.org/10.1126/science.abe7518 (2021).

    Article  ADS  CAS  PubMed  Google Scholar 

  2. Joynt, R. & Taillefer, L. The superconducting phases of UPt3. Rev. Mod. Phys. 74, 235–294 https://doi.org/10.1103/RevModPhys.74.235 (2002).

    Article  ADS  CAS  Google Scholar 

  3. Lévy, F., Sheikin, I., Grenier, B. & Huxley, A. D. Magnetic field-induced superconductivity in the ferromagnet URhGe. Science 309, 1343–1346 https://doi.org/10.1126/science.1115498 (2005).

    Article  ADS  CAS  PubMed  Google Scholar 

  4. Lounasmaa, O. V. & Thuneberg, E. Vortices in rotating superfluid 3He. Proc. Natl Acad. Sci. USA 96, 7760–7767 https://doi.org/10.1073/pnas.96.14.7760 (1999).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  5. Finne, A. et al. Dynamics of vortices and interfaces in superfluid 3He. Rep. Prog. Phys. 69, 3157–3230 https://doi.org/10.1088/0034-4885/69/12/R03 (2006).

    Article  ADS  CAS  Google Scholar 

  6. Ran, S. et al. Nearly ferromagnetic spin-triplet superconductivity. Science 365, 684–687 https://doi.org/10.1126/science.aav8645 (2019).

    Article  ADS  CAS  PubMed  Google Scholar 

  7. Aoki, D. et al. Unconventional superconductivity in UTe2. J. Phys. Condens. Matter 34, 243002 https://doi.org/10.1088/1361-648X/ac5863 (2022).

    Article  ADS  CAS  Google Scholar 

  8. Lewin, S. K., Frank, C. E., Ran, S., Paglione, J. & Butch, N. P. A review of UTe2 at high magnetic fields. Rep. Prog. Phys. 86, 114501 https://doi.org/10.1088/1361-6633/acfb93 (2023).

    Article  ADS  CAS  Google Scholar 

  9. Rosuel, A. et al. Field-induced tuning of the pairing state in a superconductor. Phys. Rev. X 13, 011022 https://doi.org/10.1103/PhysRevX.13.011022 (2023).

    Article  CAS  Google Scholar 

  10. Sakai, H. et al. Field induced multiple superconducting phases in UTe2 along hard magnetic axis. Phys. Rev. Lett. 130, 196002 https://doi.org/10.1103/PhysRevLett.130.196002 (2023).

    Article  ADS  CAS  PubMed  Google Scholar 

  11. Baek, B., Rippard, W. H., Benz, S. P., Russek, S. E. & Dresselhaus, P. D. Hybrid superconducting-magnetic memory device using competing order parameters. Nat. Commun. 5, 3888 https://doi.org/10.1038/ncomms4888 (2014).

    Article  ADS  CAS  PubMed  Google Scholar 

  12. Fermin, R., Scheinowitz, N. M. A., Aarts, J. & Lahabi, K. Mesoscopic superconducting memory based on bistable magnetic textures. Phys. Rev. Res. 4, 033136 https://doi.org/10.1103/PhysRevResearch.4.033136 (2022).

    Article  CAS  Google Scholar 

  13. Günkel, T. et al. Field-induced phase transitions in cuprate superconductors for cryogenic in-memory computing. Small 21, 2411908 https://doi.org/10.1002/smll.202411908 (2025).

    Article  CAS  Google Scholar 

  14. Cheng, Y., Shu, Q., He, H., Dai, B. & Wang, K. L. Current-driven magnetization switching for superconducting diode memory. Adv. Mater. 37, 2415480 https://doi.org/10.1002/adma.202415480 (2025).

    Article  CAS  Google Scholar 

  15. Binnig, G., Rohrer, H., Gerber, C. & Weibel, E. Surface studies by scanning tunneling microscopy. Phys. Rev. Lett. 49, 57–61 https://doi.org/10.1103/PhysRevLett.49.57 (1982).

    Article  ADS  Google Scholar 

  16. Bednorz, J. & Müller, K. Possible high Tc superconductivity in the Ba–La–Cu–O system. Z. Phys. B 64, 189–193 https://doi.org/10.1007/BF01303701 (1986).

    Article  ADS  CAS  Google Scholar 

  17. Ando, F. et al. Observation of superconducting diode effect. Nature 584, 373–376 https://doi.org/10.1038/s41586-020-2590-4 (2020).

    Article  ADS  CAS  PubMed  Google Scholar 

  18. Nadeem, M., Fuhrer, M. S. & Wang, X. The superconducting diode effect. Nat. Rev. Phys. 5, 558–577 https://doi.org/10.1038/s42254-023-00632-w (2023).

    Article  Google Scholar 

  19. Matsuki, H. et al. Realisation of de Gennes’ absolute superconducting switch with a heavy metal interface. Nat. Commun. 16, 5674 https://doi.org/10.1038/s41467-025-61267-2 (2025).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. Golod, T., Morlet-Decarnin, L. & Krasnov, V. M. Word and bit line operation of a 1 × 1 μm2 superconducting vortex-based memory. Nat. Commun. 14, 4926 https://doi.org/10.1038/s41467-023-40654-7 (2023).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  21. Semenok, D. V. et al. Superconducting memory and trapped magnetic flux in ternary lanthanum polyhydrides. Mater. Today Phys. 49, 101595 https://doi.org/10.1016/j.mtphys.2024.101595 (2024).

    Article  CAS  Google Scholar 

  22. Pippard, A. A possible mechanism for the peak effect in type II superconductors. Philos. Mag. 19, 217–220 https://doi.org/10.1080/14786436908217779 (1969).

    Article  ADS  Google Scholar 

  23. Paltiel, Y. et al. Dynamic instabilities and memory effects in vortex matter. Nature 403, 398–401 https://doi.org/10.1038/35000145 (2000).

    Article  ADS  CAS  PubMed  Google Scholar 

  24. Marchevsky, M., Higgins, M. & Bhattacharya, S. Two coexisting vortex phases in the peak effect regime in a superconductor. Nature 409, 591–594 https://doi.org/10.1038/35054512 (2001).

    Article  ADS  CAS  PubMed  Google Scholar 

  25. Xiao, Z. L., Andrei, E. Y. & Higgins, M. J. Flow induced organization and memory of a vortex lattice. Phys. Rev. Lett. 83, 1664–1667 https://doi.org/10.1103/PhysRevLett.83.1664 (1999).

    Article  ADS  CAS  Google Scholar 

  26. Ran, S. et al. Extreme magnetic field-boosted superconductivity. Nat. Phys. 15, 1250–1254 (2019).

    Article  CAS  Google Scholar 

  27. Braithwaite, D. et al. Multiple superconducting phases in a nearly ferromagnetic system. Commun. Phys. 2, 147 (2019).

    Article  Google Scholar 

  28. Aoki, D. et al. Multiple superconducting phases and unusual enhancement of the upper critical field in UTe2. J. Phys. Soc. Jpn. 89, 053705 https://doi.org/10.7566/JPSJ.89.053705 (2020).

    Article  ADS  Google Scholar 

  29. Thomas, S. M. et al. Evidence for a pressure-induced antiferromagnetic quantum critical point in intermediate-valence UTe2. Sci. Adv. 6, 8709–8723 https://doi.org/10.1126/sciadv.abc8709 (2020).

    Article  ADS  CAS  Google Scholar 

  30. Kinjo, K. et al. Superconducting spin reorientation in spin-triplet multiple superconducting phases of UTe2. Sci. Adv. 9, eadg2736 https://doi.org/10.1126/sciadv.adg2736 (2023).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wu, Z. et al. Enhanced triplet superconductivity in next-generation ultraclean UTe2. Proc. Natl Acad. Sci. USA 121, e2403067121 https://doi.org/10.1073/pnas.2403067121 (2024).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Helm, T. et al. Field-induced compensation of magnetic exchange as the possible origin of reentrant superconductivity in UTe2. Nat. Commun. 15, 37 https://doi.org/10.1038/s41467-023-44183-1 (2024).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  33. Wu, Z. et al. Superconducting critical temperature elevated by intense magnetic fields. Proc. Natl Acad. Sci. USA 122, e2422156122 https://doi.org/10.1073/pnas.2422156122 (2025).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Vasina, T. et al. Connecting high-field and high-pressure superconductivity in UTe2. Phys. Rev. Lett. 134, 096501 https://doi.org/10.1103/PhysRevLett.134.096501 (2025).

    Article  ADS  CAS  PubMed  Google Scholar 

  35. Wu, Z. et al. A quantum critical line bounds the high field metamagnetic transition surface in UTe2. Phys. Rev. X 15, 021019 https://doi.org/10.1103/PhysRevX.15.021019 (2025).

    Article  CAS  Google Scholar 

  36. Lewin, S. K. et al. High-field superconducting halo in UTe2. Science 389, 512–515 https://doi.org/10.1126/science.adn7673 (2025).

    Article  ADS  CAS  PubMed  Google Scholar 

  37. Wu, Z. et al. Magnetic signatures of pressure-induced multicomponent superconductivity in UTe2. Phys. Rev. Lett. 134, 236501 https://doi.org/10.1103/PhysRevLett.134.236501 (2025).

    Article  ADS  CAS  PubMed  Google Scholar 

  38. Tokunaga, Y. et al. Longitudinal spin fluctuations driving field-reinforced superconductivity in UTe2. Phys. Rev. Lett. 131, 226503 https://doi.org/10.1103/PhysRevLett.131.226503 (2023).

    Article  ADS  CAS  PubMed  Google Scholar 

  39. Kamat, S., et al. Thermodynamic discovery of tetracriticality and emergent multicomponent superconductivity in UTe2. Preprint at https://arxiv.org/abs/2603.17905 (2026).

  40. Yang, Z. Spectroscopic evidence of symmetry breaking in the superconducting vortices of UTe2. Natl Sci. Rev. 12, nwaf267 https://doi.org/10.1093/nsr/nwaf267 (2025).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Sharma, N. et al. Observation of persistent zero modes and superconducting vortex doublets in UTe2. ACS Nano 19, 31539–31550 https://doi.org/10.1021/acsnano.5c08406 (2025).

    Article  CAS  PubMed  Google Scholar 

  42. Yin, R., et al. Yin-Yang vortex on UTe2 (011) surface. Nat. Commun. 17, 5394 https://doi.org/10.1038/s41467-026-72162-9 (2026).

  43. Wang, Y. et al. Observation of vortex stripes in UTe2. Nano Lett. 25, 12824–12831 https://doi.org/10.1021/acs.nanolett.5c02265 (2025).

    Article  ADS  CAS  PubMed  Google Scholar 

  44. Ishihara, K. et al. Anisotropic enhancement of lower critical field in ultraclean crystals of spin-triplet superconductor candidate UTe2. Phys. Rev. Res. 5, L022002 https://doi.org/10.1103/PhysRevResearch.5.L022002 (2023).

    Article  CAS  Google Scholar 

  45. Tokiwa, Y. et al. Anomalous vortex dynamics in the spin-triplet superconductor UTe2. Phys. Rev. B 108, 144502 https://doi.org/10.1103/PhysRevB.108.144502 (2023).

    Article  ADS  CAS  Google Scholar 

  46. Tinkham, M. Introduction to Superconductivity (Courier Corporation, 2004).

  47. Green, M. A. Intrinsic concentration, effective densities of states, and effective mass in silicon. J. App. Phys. 67, 2944–2954 https://doi.org/10.1063/1.345414 (1990).

    Article  ADS  CAS  Google Scholar 

  48. Buck, D. A. The cryotron–a superconductive computer component. Proc. IRE 44, 482–493 https://doi.org/10.1109/JRPROC.1956.274927 (1956).

    Article  ADS  Google Scholar 

  49. Zhang, L. et al. Dimensionality of the reinforced superconductivity in UTe2. Nat. Commun. 16, 10308 https://doi.org/10.1038/s41467-025-66288-5 (2025).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  50. Yoon, H., Baek, S., Saha, S. R., Butera, R. E. & Paglione, J. Submicrometer-thick UTe2 flake achieved by mechanical exfoliation. Supercond. Sci. Technol. 39, 045005 https://doi.org/10.1088/1361-6668/ae5742 (2026).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  51. Cao, Y., Park, J. M., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Pauli-limit violation and re-entrant superconductivity in moiré graphene. Nature 595, 526–531 https://doi.org/10.1038/s41586-021-03685-y (2021).

    Article  ADS  CAS  PubMed  Google Scholar 

  52. Rubi, K. et al. High-field-stabilized reentrant superconductivity in infinite-layer nickelate thin films. Nat. Commun. 17, 9267 https://doi.org/10.1038/s41467-026-75922-9 (2026).

  53. Sakai, H. et al. Single crystal growth of superconducting UTe2 by molten salt flux method. Phys. Rev. Mater. 6, 073401 https://doi.org/10.1103/PhysRevMaterials.6.073401 (2022).

    Article  CAS  Google Scholar 

  54. Eaton, A. G. et al. Quasi-2D Fermi surface in the anomalous superconductor UTe2. Nat. Commun. 15, 223 https://doi.org/10.1038/s41467-023-44110-4 (2024).

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  55. Li, G. et al. An experimental station with high-field all-superconducting magnet focusing on quantum oscillation studies. Chin. Phys. B https://doi.org/10.1088/1674-1056/ae7e75 (2026).

  56. Jermain, C., Rowlands, G. et al. PyMeasure package. https://pymeasure.readthedocs.io/en/latest/index.html (2023).

  57. Wu, Z. & Eaton, A. G. Research data supporting: electrically controllable superconducting memory effect in UTe2. https://doi.org/10.17863/CAM.131466 (2026).

Download references

Acknowledgements

We gratefully acknowledge stimulating discussions with D. Agterberg, E. Babaev, P. Coleman, J. Durrell, A. Greer, A. Huxley, M. Ireland, P. Littlewood, G. Lonzarich, C. Reichhardt, T. Winyard and especially A. Levchenko.

Funding

This project was supported by the UK Engineering and Physical Sciences Research Council (EPSRC) through grants EP/X011992/1 and EP/R513180/1. Part of this work was carried out at the Synergetic Extreme Condition User Facility (SECUF, https://cstr.cn/31123.02.SECUF). Part of this work was performed at the National High Magnetic Field Laboratory, which is supported by National Science Foundation (NSF) cooperative agreement nos. DMR-1644779 and DMR-2128556 and the State of Florida. We acknowledge financial support by the Czech Science Foundation GAČR under the JUNIOR STAR grant no. 26-21795M (STiUS). The work of D.S. was financially supported by the NSF Quantum Leap Challenge Institute for Hybrid Quantum Architectures and Networks grant no. OMA-2016136. D.V.C. acknowledges financial support from the National High Magnetic Field Laboratory through a Dirac Fellowship, which is financed by the NSF (grant no. DMR-2128556) and the State of Florida, and from Washington University in St. Louis through the Edwin Thompson Jaynes Postdoctoral Fellowship. T.I.W. acknowledges support from Murray Edwards College (University of Cambridge) and the Cambridge Philosophical Society through a Henslow Fellowship. T.I.W. and A.G.E. acknowledge support from the Institute for Complex Adaptive Matter (ICAM) through US NSF grant number 2201516 under the AccelNet programme of the Office of International Science and Engineering and from QuantEmX grants from ICAM and the Gordon and Betty Moore Foundation through grants GBMF5305 and GBMF9616. G.L. acknowledges support from the National Key Research and Development Projects of China (grant nos. 2022YFA1602803 and 2023YFA1607402). R.Z. acknowledges support from the National Key Research and Development Projects of China (grant nos. 2024YFA1611302 and 2025YFE0202100) and CAS PIFI programme (2024PG0003). A.G.E. acknowledges support from the Henry Royce Institute for Advanced Materials through the Equipment Access Scheme enabling access to the Advanced Materials Characterisation Suite at Cambridge, grant numbers EP/P024947/1, EP/M000524/1 and EP/R00661X/1, and from Sidney Sussex College (University of Cambridge).

Author information

Authors and Affiliations

  1. Cavendish Laboratory, University of Cambridge, Cambridge, UK

    Zheyu Wu, Hanyi Chen, Mengmeng Long, Theodore I. Weinberger, Alexander J. Hickey, F. Malte Grosche & Alexander G. Eaton

  2. National High Magnetic Field Laboratory, Tallahassee, FL, USA

    Zheyu Wu, Dmitry V. Chichinadze & Dave Graf

  3. Department of Physics, University of Wisconsin–Madison, Madison, WI, USA

    Daniel Shaffer

  4. Department of Physics, Washington University in St. Louis, St. Louis, MO, USA

    Dmitry V. Chichinadze

  5. Department of Condensed Matter Physics, Faculty of Mathematics and Physics, Charles University, Prague 2, Czech Republic

    Andrej Cabala, Vladimir Sechovský & Michal Vališka

  6. School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai, China

    Jinxu Pu

  7. Clarendon Laboratory, Department of Physics, University of Oxford, Oxford, UK

    Jinxu Pu

  8. Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing, China

    Gang Li & Rui Zhou

  9. School of Physical Sciences, University of Chinese Academy of Sciences, Beijing, China

    Gang Li & Rui Zhou

Authors

  1. Zheyu Wu
  2. Hanyi Chen
  3. Mengmeng Long
  4. Daniel Shaffer
  5. Dmitry V. Chichinadze
  6. Andrej Cabala
  7. Theodore I. Weinberger
  8. Alexander J. Hickey
  9. Jinxu Pu
  10. Dave Graf
  11. Vladimir Sechovský
  12. Michal Vališka
  13. Gang Li
  14. Rui Zhou
  15. F. Malte Grosche
  16. Alexander G. Eaton

Contributions

A.C., V.S. and M.V. grew single-crystal UTe2 specimens. Z.W. discovered the memory effect. Z.W., H.C., M.L., T.I.W., A.J.H., J.P., D.G., G.L., R.Z. and A.G.E. performed measurements. Z.W., H.C., M.L. and A.G.E. analysed data. Z.W. and A.G.E. interpreted data. D.S. and D.V.C. performed theoretical modelling. F.M.G. and A.G.E. supervised the project. A.G.E. wrote the paper, with input from all co-authors.

Corresponding author

Correspondence to Alexander G. Eaton.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature thanks Lei Chen, Wade DeGottardi, Yun Suk Eo and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Saturation of the memory effect.

a–c, Amplitude (a), duration (b) and relaxation rate tuning (c) of the memory effect of UTe2. For sufficiently large perturbative amplitudes Jex or durations Δt, the memory effect approaches a saturation value. ΔJ is calculated here at a nominal fixed value of V = 1 μV, as indicated by the dashed lines. All data were acquired on sample S2 at 14 T and 50 mK.

Extended Data Fig. 2 Further temperature dependence of the memory effect.

a,b, Results for duration tuning (a) and relaxation rate tuning (b), to complement the amplitude tuning depicted in Fig. 3. A similar evolution, of diminishing memory effect at higher temperatures, is seen for all three types of perturbative stimuli.

Extended Data Fig. 3 Evolution of the memory effect as a function of magnetic field strength.

a–c, Amplitude (a), duration (b) and relaxation rate tuning (c) of the memory effect. Measurements were performed on sample S2 for B aligned along the \(\hat{b}\text{-axis}\).

Extended Data Fig. 4 Magnetic field dependence of the slope of V(J) in the equilibrium state.

\(\frac{\partial V}{\partial J}\) at incremental field strengths as indicated. Jc is minimal around 20 T and is actually higher at 28 T than it is at 9 T. Anomalous non-monotonic behaviour is observed in the memory region around 20 T. Note that, for this sample (S1), the memory region is located over a higher range of B than for sample S2, probably because of differences in crystalline quality (and hence vortex pinning forces). All measurements were performed at T = 50 mK.

Extended Data Fig. 5 Rotation study of the superconducting memory effect.

a–c, Amplitude (a), duration (b) and relaxation rate tuning (c) of the memory effect under rotations of B from the \(\hat{b}\text{-axis}\) (θ = 0°) towards the \(\hat{c}\text{-axis}\) (θ = 90°). ΔJ abruptly diminishes for θ ⪆ 19°.

Extended Data Fig. 6 Tuning of memory hysteresis.

Controllability of the superconducting memory effect by modulation of the perturbation amplitude. Insets show the chronological sequence of how J was swept. In Fig. 1, we always perturbed the system by discontinuously changing J from the same magnitude—by contrast, here we perform the perturbation from two different values of J. A larger hysteresis in V(J) is observed for the larger perturbation. This shows that the magnitude of the memory effect is tunable by the size of the excitation, similar to Fig. 2.

Extended Data Fig. 7 Presence of the memory effect under both polarities of current flow.

a,b, Results for positive (a) and negative (b) polarity of J. The same (positive) polarity of erasure and perturbation were performed in both instances, with the response then measured upon sweeping J to large positive and to large negative values (as depicted in the schematics on the left). Consistent behaviour is observed for both polarities.

Extended Data Fig. 8 Raw temperature-dependent data.

Hysteretic J–V measurements performed on sample S2 at 14 T at 50 mK, 100 mK, 300 mK, 430 mK and 600 mK for amplitude, duration and relaxation rate tuning. These data were used to generate the heatplots presented in Fig. 3 and Extended Data Fig. 2.

Extended Data Fig. 9 Rotation study showing no memory effect in SC2.

V(J) at B = 29 T for rotations of B from the \(\hat{b}\text{-axis}\) (θ = 0°) towards the \(\hat{c}\text{-axis}\) (θ = 90°). The measurement protocol is specified in the top-left inset. No memory effect is discerned at any θ.

Supplementary information

Source data

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, Z., Chen, H., Long, M. et al. Electrically controllable superconducting memory effect in UTe2. Nature 657, 632–637 (2026). https://doi.org/10.1038/s41586-026-11015-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • Version of record:

  • Issue date:

  • DOI: https://doi.org/10.1038/s41586-026-11015-3