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.
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.
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
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).
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Ran, S. et al. Extreme magnetic field-boosted superconductivity. Nat. Phys. 15, 1250–1254 (2019).
Article CAS Google Scholar
Braithwaite, D. et al. Multiple superconducting phases in a nearly ferromagnetic system. Commun. Phys. 2, 147 (2019).
Article Google Scholar
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
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
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
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
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
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
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
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
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
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
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
Kamat, S., et al. Thermodynamic discovery of tetracriticality and emergent multicomponent superconductivity in UTe2. Preprint at https://arxiv.org/abs/2603.17905 (2026).
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
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
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).
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
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
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
Tinkham, M. Introduction to Superconductivity (Courier Corporation, 2004).
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
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
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
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
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
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).
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
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
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).
Jermain, C., Rowlands, G. et al. PyMeasure package. https://pymeasure.readthedocs.io/en/latest/index.html (2023).
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).
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
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and permissions
About this article
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