Moving flux quanta cool superconductors by a microwave breath

Almost any use of a superconductor implies a nonequilibrium state. Remarkably, the non-equilibrium states induced by a microwave stimulus and the dynamics of magnetic flux quanta (Abrikosov vortices) can give rise to strikingly contrary effects: A sufficiently high-power electromagnetic field of GHz frequency can stimulate superconductivity, whereas fast vortex motion can trigger an instability abruptly quenching the superconducting state. Here, we advance or delay such dynamical quenching of the vortex state in Nb thin films by tuning the power and frequency of the microwave ac stimulus added to a dc bias current. The experimental findings are supported by time-dependent Ginzburg-Landau simulations and they can be explained qualitatively based on a model of"breathing mobile hot spots", implying a competition of heating and cooling of quasiparticles along the trajectories of moving fluxons whose core sizes vary in time. In addition, we demonstrate universality of the stimulation effect on the thermodynamic and transport properties of type II superconductors.

Remarkably, the non-equilibrium states in superconductors generated by a microwave stimulus and the Abrikosov vortex dynamics can imply conceptually opposite effects. Thus, while superconductivity is known to be destroyed by high temperatures, currents and magnetic fields, a sufficiently high-power electromagnetic field of subgap frequency may stimulate superconductivity itself [30]. This counterintuitive effect was explained as a consequence of an irradiation-induced redistribution of quasiparticles away from the superconducting gap edge * oleksandr.dobrovolskiy@univie.ac.at [26]. While microwave-stimulated superconductivity was since then observed in various systems [31][32][33], no theory of this effect in the presence of vortices is available so far. At the same time, investigations of stimulated superconductivity in the vortex state could have important implications both inside and outside the condensed matter physics community, as mapping solutions of astrophysics problems to scalar condensates allows for modeling the physics of black holes and gravity upon holographic superconductors [34]. While the latter can exhibit vortices [35] and vortex-flow states [36], the possibility of stimulation of holographic superconductors is under current debate [37,38].
The presence and motion of vortices is expected to modify the signatures of microwave-stimulated superconductivity due to the variation of the order parameter in space and time, the friction-induced heating of quasiparticles in the vortex cores [27,39,40], and the energy leaking outside of the core at large vortex velocities [20][21][22]. The latter results in a vortex-core shrinkage, a further acceleration of the flux flow, and the associated quench of the low-resistive state due to the flux-flow instability [20][21][22]. While superconductivity in the presence of a high-frequency excitation and non-equilibrium effects related to vortex dynamics have recently received much attention theoretically [29,[41][42][43][44][45] and experimentally [46][47][48][49][50][51][52][53][54][55][56], the microwave-stimulated and the vortex-dynamicsgenerated nonequilibrium states have never been studied simultaneously so far.
Here, we investigate the competition between quenching and stimulation of superconductivity in the GHzfrequency (ac+dc)-driven nonlinear resistive regime in Nb thin films. Under optimized excitation conditions we demonstrate a pronounced extension of the lowdissipative state towards higher vortex velocities as compared to the unexcited regime. Our experimental findings are largely reproduced by time-dependent Ginzburg-Landau simulations and can be explained by a competition of heating and cooling of quasiparticles escaping from vortices in conjunction with a periodic variation of their core size. Additionally, we demonstrate universality of the stimulation effect on the thermodynamic and transport properties of type II superconductors.

I. EXPERIMENTAL RESULTS
Microwave control of the vortex dynamics. We study the vortex dynamics under superimposed direct and microwave current drives in a coplanar waveguide (CPW) made of a Nb film with thickness d = 50 nm and exhibiting a superconducting transition at T c = 8.544 K. The experimental geometry is shown in Fig. 1(a). The perpendicular-to-film-plane magnetic field with induction B = µ 0 H populates the CPW with a lattice of Abrikosov vortices. The sum of applied direct and microwave currents exerts a Lorentz force on the vortices that causes their motion with velocity v across the central conductor of the CPW. The associated voltage drop along the central conductor allows one to distinguish between the different resistive states of the sample.
We demonstrate the control of the nonlinear resistive regime via advancing or delaying the dc-bias-induced breakdown of the non-equilibrium superconducting state. This control can be clearly seen in Fig. 1(c-e) where the current-voltage (I-V ) curves at T = 0.988T c are shown for a series of magnetic fields at three ac frequencies for the ac power levels −60 dBm (corresponding to 1 nW, solid lines) and −6 dBm (correspodning to 0.25 mW, symbols). In what follows, the excitation level −60 dBm will be referred to as the unexcited state. Thus, in all I-V curves one can recognize the nearly linear regime of flux flow followed by an upward bending at the foots of abrupt jumps to the normally conducting state. These jumps are the hallmark of the Larkin-Ovchinnikov (LO) instability [20][21][22] occurring at the instability current I * relating to the instability voltage V * . The definition of I * and V * , as well as the respective quantities in the presence of a microwave current, I * mw and V * mw , is shown in Fig. 1(e).
The most striking observation in Fig. 1(c-e) is that depending on the ac frequency, the microwave ac stimulus affects the onset of the flux-flow instability differently. At 64.1 MHz in Fig. 1(c), which is exemplary for relatively low frequencies, the instability jumps occur at I * mw (H) < I * (H). The occurrence of the instability at smaller I values in the excited I-V curves can be understood as a consequence of the replacement of I by the sum of direct and microwave currents in the excited regime. By contrast, at 13.9 GHz, which is representative for the highest frequencies available in our experiment, the onset of the instability is shifted towards higher current values. This shift is accompanied by the development of a pronounced upturn bending in the excited I-V curves as compared to the unexcited ones. The expansion of the nonlinear regime in the excited I-V curves in Fig.  1(e) can be clearly seen as a deviation from the reference dashed straight line. These extended nonlinear I-V sections at the foots of the instability jumps are the fingerprint of the stabilization effect of the microwave current on the flux flow in the investigated system. For completeness, in Fig. 1(d) we present the data at 1.02 GHz which has been chosen as an intermediate frequency at which the excited and unexcited I-V curves almost perfectly overlap. This coincidence can be explained by the assumption that the contribution of the microwave current into the triggering of the instability is nearly completely compensated by its stabilization effect on the flux flow.
Microwave control of instability parameters. To illustrate further the stabilization effect of the microwave current on the nonequilibrium state generated by moving vortices, from the last data point before the instability jump we deduce the instability velocity v * and the instability current density j * by the standard relations v * = V * /(HL) and j * = I * /(wd). Here, L = 1 mm is distance between the voltage contacts while w = 50 µm and d = 50 nm is the conductor width and thickness, respectively. The magnetic field dependences of the instability parameters deduced from the I-V curves are presented in Fig. 2(a) and (b). One clearly sees a decrease of v * (H) and j * (H) for each dataset, as well as v * (f ) and j * (f ) at a fixed H value. To describe the evolution of the instability parameters in the presence of the microwave current quantitatively, we introduce new parameters ∆v * = (v * mw − v * )/v * and ∆j * = (j * mw − j * )/j * for their relative changes. Their field dependence at various microwave frequencies is shown in Fig. 2(d) and (e). From the field dependences of ∆v * (H) and ∆j * (H) it follows that the enhancement of the critical velocity in the presence of a microwave current can reach up to 40 % while for the instability current the maximum enhancement can be about 10 %. Remarkably, the enhancement of the instability parameters attains a maximum at about 10 mT, and the stimulation effect becomes weaker with decrease of the ac power and vanishes upon reaching an ac power level of about −20 dBm. This is clearly different from the excitation of the sample with a frequency of 64.1 MHz at which the instability parameters are reduced by up to about 7% with respect to the unexcited state.
The evolution of the critical velocity in the broad range of ac power levels and frequencies is summarized in Fig.  2(f). Remarkably, one can distinguish two frequency regimes, f 1 GHz and f 1 GHz, in which the behavior of the v * mw /v * differs qualitatively. In the regime f 1 GHz, which we will call the low-frequency regime, the microwave current either has no effect on v * mw /v * at lower ac power levels, or it suppresses the flux-flow stability at higher power levels, P > P sup (f ). By contrast, at f 1 GHz, which will be referred to as a high-frequency regime, the v * mw /v * ratio first remains constant, then grows in the power range P st (f ) < P < P sup (f ), and finally decreases at yet higher microwave power levels P > P sup (f ). Here, P st (f ) and P sup (f ) are the frequencydependent ac power levels above which a stimulation and a suppression effect are observed, respectively.
A similar stimulation effect for the maximum current I * mw , up to which the low-resistive state is maintained in the presence of a microwave ac stimulus, is demonstrated in Fig. 3. In Fig. 3(a) the frequency dependence of the normalized current I * mw /I * exhibits a systematic evolution from a weak frequency dependence at low ac power levels to a notably growing tendency at higher power levels. While the dashed lines in each panel in Fig. 3(a) depict the reference level I * mw /I * = 1 in the absence of a microwave excitation, an enhancement of the instability current is clearly seen at higher frequencies and power levels exceeding about −12 dBm. At the same time, I * mw /I * begins to decrease with a further increase of the microwave power. The dependence I * mw (P, f )/I * in Fig. 3(b) and (c) is qualitatively very similar to that for v * mw (P, f )/v * in Fig. 2(f). Furthermore, both depen-  Interestingly, the maximal increase of the transition temperature in the presence of an ac stimulus is only of the order of 1%, while that for H c2 reaches about 6%. At the same time, the relative enhancement of the instability velocity and the instability current is significantly larger and reaches up to 40% and 10%, respectively.
The effect of a microwave ac stimulus on the critical transport and thermodynamic parameters can be summarized in the f -P phase diagram of the instability current in Fig. 3(c), where four different regimes can be identified. Namely, at low ac frequencies and power levels there is virtually no effect of the microwave stimulus, as expected (region 1 in Fig. 3(c)). With increasing ac power, at 1 GHz an enhancement of the parameters is observed (region 2 in Fig. 3(c)). At a further increase of the ac power up to P h the parameters at higher frequencies stop to increase even further, while at low frequencies the parameters begin to decrease (region 3 in Fig. 3(c)). Finally, at high ac power levels a suppression of the superconducting critical parameters is observed regardless of the ac frequency (region 4 in Fig. 3(c)). In all panels T = 0.988Tc. fst and Pst are the ac frequency and the microwave power, respectively, above which H mw c2 becomes larger than Hc2. Psup designates the microwave power above which H mw c2 becomes smaller than Hc2.

SIMULATION RESULTS
Since analytical theories of microwave-stimulated mixed state and the flux-flow instability in the presence of a high-frequency ac current are unavailable so far, further insights into the evolution of the dynamic state generated by vortex motion can be gained on the basis of computer simulations relying upon the solution of the time-dependent Ginzburg-Landau (TDGL) equation. While the microscopic derivation of the TDGL was originally done for a gapless superconductor with paramagnetic impurities [57], it is also widely used for studying various aspects of current-driven vortex matter in superconductors with gap and non-magnetic impurities [58,59], including the vortex dynamics at high vortex velocities and at GHz ac frequencies [16,60,61]. In the TDGL, not only the vortex-vortex interaction is taken into account, but also vortex-core structures and interaction with pinning centers can be described. At the same time, far from T c , the TDGL equationdoes not reproduces the physics in the vortex core quantitatively, but it still describes the spatiotemporal evolution of vortex matter qualitatively [59]. For the simulations we adopt the link variable method [51,62] for 2D and 3D systems of adjustable size and shape and use a solver that approximates numerically the solution of the TDGL equation. The gauge for the TDGL solver is such that the scalar potential is zero, so no electric currents can be introduced directly in the simulations and the effect of the currents is introduced by the stray field they would generate. The simulations take place in a two-dimensional 360 × 100 cell superconducting strip, with a constant and uniform external field H = 0.2H c2 applied perpendicular to the plane of the strip. The magnetic field H dc of opposite directions at each boundary is applied to simulate the field generated by the dc bias current. The magnitude of H dc is varied in steps of 0.005H c2 every 10 4 steps of the simulation, giving enough time to reach a stationary state before increasing H dc further. Finally, an additional field H mw is applied with opposite and alternating directions at each boundary, simulating the high-frequency ac current of constant amplitude and frequency f , that affects the dynamic state of the strip. The geometry of the applied fields is illustrated in Fig. 6(a). More details on the simulation procedure are given in the Methods section.
The simulations are done for Nb with the Ginzburg-Landau parameter κ = 6 deduced from the experiment and the critical temperature T c = 8.5 K. With these values, simulations were performed at T = 4 K with ac frequencies f = 0.1f 0 , 0.5f 0 , f 0 , 5f 0 and 10f 0 and with an ac field amplitude of 0.2H c2 . Here, f 0 = 1/t 0 , where t 0 = π /96k B T c is the characteristic time of the relaxation of the order parameter [63]. The low temperature was chosen for the only reason of a better contrast in the spatial dependence of the order parameter in the simulations of the vortex patterns. After setting the temperature, the external magnetic field and the main parameters defining the material, the order parameter is initialized in the superconducting state |ψ| = 1 and its evolution is simulated as a function of the dc current in the presence of a high-frequency ac current stimulus, Fig. 6(b). Soon after start of the simulations the Nb strip is quickly filled with vortices, due to the applied external field. As soon as the dc current begins to increase, there appears a net flux of vortices in the direction of the y axis, from top to bottom. A further increase of the dc current causes a shift of vortices that turns the upper part of the strip into the normal state, and eventually breaks the superconducting channel. The different stages of this process can be seen in Fig. 6(c) illustrating the vortex patterns and the evolution of the superconducting channel with increasing dc current value. In the simulations, the critical current I c is defined as a current at which the super- conducting channel breaks down. The relative variation of this current is measured as a function of the normalized drive frequency f , exhibiting a 60% increase as the frequency is varied from 0.1f 0 to 10f 0 . The increase of I c (f ) is nonlinear, with a faster increase between 0.1f 0 and f 0 , see Fig. 6(f). We have also checked that the acstimulated increase of I c varies with κ, as can be seen in Fig. 6(g). Thus, for κ = 1 the variation is smaller than about 3.5 %, while for κ = 6 it can reach 60 %. The dependences of the critical current on the microstrip aspect ratio and temperature are illustrated in Fig. 7.

DISCUSSION
The reported experimental findings relate to a fundamental problem of the interaction of a transport current containing dc and high-frequency ac components with a superconductor in the vortex state. As is well known, the vortex state is characterized by a spatially modulated superconducting order parameter which vanishes in the vortex cores and attains a maximal value between them. Accordingly, a type II superconductor can be regarded as a continuum medium consisting of bunches of quasiparticles in the vortex cores surrounded by a bath of the superconducting condensate formed by the superfluid of Cooper pairs [64]. To the authors' best knowledge, so far there is no available theory addressing the complex interaction of the superimposed dc and microwave currents with the quasiparticles and the condensate at the same time. Yet, some ingredients of this interaction, which were studied separately so far, should be mentioned prior to discussing the experimental findings.
The effect of a dc current on the superconducting condensate is well known [65]. With an increase of a dc current, the absolute value of the order parameter is decreasing and the peak in the density of states at the edge of the superconducting gap is smeared. This is because of gaining a finite momentum by the Cooper pairs that form a coherent excited state that plays a central role in the explanation of the gauge invariance of the Meiss-ner effect [66]. The equivalence of depairing due to an electric current and due to a magnetic field is also well known, both theoretically [67] and experimentally [65], and we have in fact used this equivalence in the TDGL simulations.
A general theory of depairing by a microwave field was formulated quite recently [41]. It was shown that the ground state of a superconductor is altered qualitatively in analogy to the depairing due to a dc current. However, in contrast to dc depairing, the density of states acquires steps at multiples of the microwave photon energy and shows an exponential-like tail in the subgap regime [41]. Additionally, depending on temperature, one can consider two regimes in which the response of a superconductor is dominated either by the response of the superfluid (at low temperatures, T /T c ≪ 1) or by the quasiparticles (close to the critical temperature, (T − T c )/T c ≪ 1, as in our experiment). It is also known that at T T c , microwave radiation can be absorbed by quasiparticles, leading to a nonequilibrium distribution over the energies [26].
In general, the response of the condensate to an external microwave field becomes apparent via a change of the kinetic impedance (imaginary part of the complex resistivity) while the quasiparticles give rise to the microwave loss (real part of the complex resistivity). In the presence of an external magnetic field inducing vortices in the superconductor, the vortex-induced resistive loss dominates the response of the superconductor. We note that the ac frequencies in our experiment are much smaller than the superconducting gap frequency f ≪ f gap (0.998T c ) ≃ 75 GHz and the experiment is done in the vicinity of T c where the superconducting gap is small and the GL, LO, and Eliashberg theories are justified [21,26,57,59,68,69].
It is important to stress that the known effects of stimulation of superconductivity by a microwave stimulus in the absence of an external magnetic field include an enhancement of the Ginzburg-Landau depairing current (in narrow channels) implying a transition to a resistive state due to the formation of phase-slip centers [31] or the Aslamazov-Lempitskii maximum current (in wide films) at which the vortex structure induced by the self-field evolves into the first phase-slip line [70]. By contrast, in the presence of an external magnetic field when Abrikosov vortices move under the action of the transport current, there is an additional phenomenon leading to an abrupt quenching of the superconductor to the normally conducting state earlier than the Ginzburg-Landau or Aslamazov-Lempitskii critical current is reached. Within the LO theoretical framework [20,21], this quenching is because of the flux-flow instability caused by the nonlinear dependence of the film conductivity on the electric field. However, the LO theory was developed in the dirty limit near T c and for weak magnetic fields when the heating of the superconductor can be neglected. To account for a finite density of vortices, the LO theory was extended by Bezuglyj and Shklovskij in Ref. [22]. In that work [22], the complete set of instability points in I-V curves for a series magnetic field values is described by the system of equations where t = [1+b+(b 2 +8b+4) 1/2 ]/3(1+2b) and b = B/B T is the dimensionless magnetic field with the parameter Here, k B is the Boltzmann constant, c the speed of light in vacuum, R = (σ n d) −1 the film sheet resistance, e 0 the electron charge, h the heat removal coefficient, and τ ε the quasiparticle energy relaxation time. In the system of equations (1), the parameters E 0 and j 0 are defined as The curve calculated by Eqs. (1) and (3) is shown in Fig. 2(c) which also contains the experimentally measured instability points in the normalized voltage V * /V 0 ≡ E * /E 0 versus normalized current density j * /j 0 representation. The normalization parameters used for fitting the experimental data to the theoretical curve are B T = 12 mT with V 0 = 0.155; 0.14; 0.17 V/cm and j 0 = 184; 176; 192 kA/cm 2 for the unexcited state and the excitation with 64.1 MHz and 13.9 GHz, respectively. From the figure it follows that the instability points E * (j * /j 0 )/E 0 in the absence of ac current and in the presence of an ac current with f = 64.1 MHz nicely fit to the expressions (1). By contrast, in the presence of an ac current with f = 13.9 GHz a noticeable deviation from the expressions (1) is observed. While it was recently revealed that possible inhomogeneities in the distribution and strength of pinning sites do not alter the dependence given by Eq. (1) qualitatively [55], but rather require a renormalization of the parameters V 0 and E 0 , our findings suggest that the superimposed microwave ac current at high enough power levels becomes a crucial ingredient which qualitatively modifies the physical picture and causes a deviation from the flux-flow instability theory developed for the sole case of a dc bias current. Namely, the microscopic scenario of the flux-flow instability under a dc current drive implies a decrease in the number of quasiparticles in the vortex cores under the action of an electric field. In return, the decrease in the number of quasiparticles leads to a shrinkage of the vortex cores and a decrease in the vortex viscosity with increasing vortex velocity. As a consequence, the viscous force has a maximum as a function of the vortex velocity, and as soon as the Lorentz force exceeds this maximum, the viscous flow of the vortices becomes unstable.
From the specific power at the instability point, P 0 = j 0 E 0 = (h/d)(T c − T ) [22,48,55,71], following from Eqs. (3) with σ n = 1/(R d), one can deduce the heat removal coefficient h ≈ 1.2 WK −1 cm −2 . Substitution of h and B T = 12 mT into Eq. (2) yields the energy relaxation time τ ε ≈ 0.23 ns. For the relaxation time associated with the electron-phonon scattering in the LO model, we find this estimate to be in reasonable agreement with the τ ε estimates of 0.15 ns [72], 0.2-0.4 ns [48], and 0.3-0.7 ns [46] for Nb thin films. We note that if one uses the order parameter relaxation time (1.8 · 10 −11 s for Nb [68,72]) for the estimate of the minimal frequency f min = 1.73/(2πτ ε ) [68], above which an enhancement of superconductivity via the Eliashberg mechanism is possible, than no microwave stimulation effect is expected at f < f min ≈ 15 GHz [68]. However, for gapped superconductors, such as Nb, it was argued that the gap change is much slower, with a relaxation rate dominated by the electron-phonon scattering [73,74]. Accordingly, if one takes τ ε = 0.23 ns for the estimate of the threshold frequency in Nb, one obtains f min ≈ 1.2 GHz. While the microscopic description of the studied system should clearly go beyond the Eliashberg theory [26,68], we underline that the enhancement of the critical parameters in our experiment is observed at f 1 GHz. In this way, the electron-phonon interaction time appears to be an important ingredient governing microwave-stimulated superconductivity in the vortex state in Nb films.
To elucidate the extension of the low-resistive flux-flow regime in the presence of (dc+ac) current drives, we suggest the following qualitative explanation. Namely, when an ac current is added to the dc current, at rather high microwave power levels the vortices start oscillating near their equilibrium positions, which in return are displaced due to vortex motion under the action of the dc driving current. In other words, the dc current leads to a translational motion of vortices while the ac stimulus induces a breathing mode in addition to their translation. This breathing mode appears due to the variation of the vortex core sizes in time due to the periodically-modulated quasiparticle escape from the cores. Accordingly, the strong oscillations of vortices are expected to lead to the formation of "clouds" of quasiparticles around the vortex cores, whose relaxation should now take place in a larger volume as compared to the non-excited case. At the same time, the dissipation is decreasing due to the combined effects of the shrinkage of the vortex cores and the redistribution of quasiparticles away from the gap edge via the Eliashberg mechanism. While a detailed theory of microwave-stimulated superconductivity in the vortex state at high vortex velocities is yet to be elaborated, we believe that the model of "moving breathing hot spots" can be applied for a qualitative explanation of the observed effects.
To summarize, we have investigated the effect of an ac current with frequencies in the MHz and the lower GHz range on the resistive state of superconducting Nb films subjected to a dc bias current in perpendicular magnetic fields. While without ac excitation and in the presence of an ac excitation at frequencies in the MHz range the quenches of the superconductor to a highly-resistive state are nicely described by the expressions derived within the framework of the Larkin-Ovchinnikov theory of flux-flow instability, the addition of a GHz-frequency ac current at moderately high power levels leads to a notable deviation from the Larkin-Ovchinnikov theory and to an extension of the low-resistive dynamical state to larger current values. This key experimental observation is supported by simulations based on the time-dependent Ginzburg-Landau equation and can be explained qualitatively by a model of "breathing mobile hot spots" implying a competition of heating and cooling of quasiparticles along the trajectories of moving fluxons whose core sizes vary in time. While a complete theory of microwave-stimulated superconductivity in the vortex state remains unavailable so far, our results are relevant for superconducting microwave circuits exploited in quantum computing and remote sensing, as well as they furthermore might have implications for diverse physical problems allowing mapping of their solutions to scalar condensates.

A. Film growth and characterization.
The CPWs were fabricated by photolithography and Ar etching from epitaxial (110) Nb films on a-cut sapphire substrates. The films were grown by dc magnetron sputtering in a setup with a base pressure in the 10 −8 mbar range. The films were sputtered at the substrate temperature T = 800 • C. The Ar pressure was 5 × 10 −3 mbar and the film growth rate was 1 nm/s. The (110) orientation of the films was inferred from X-ray diffraction measurements [75]. The epitaxy of the films was confirmed by reflection high-energy electron diffraction. The as-grown films have a smooth surface with an rms surface roughness of less than 0.5 nm, as inferred from atomic force microscopy inspection in the range 1 µm×1 µm. The films are characterized by a superconducting transition temperature T c in zero field of 8.544 K, as determined by the 90 % resistance criterion. Their normal-state resistivity just above the superconducting transition amounts to 2.48 µΩcm and the upper critical field at zero temperature is estimated as H c2 (0) = H c2 (T )/[1 − (T /T c ) 2 ] ≈ 1.2 T corresponding to a superconducting coherence length ξ(0) = (Φ 0 /2πH c2 ) 1/2 ≈ 17 nm. Here, Φ 0 is the magnetic flux quantum. The magnetic field penetration depth λ(0) at zero temperature in the films can be estimated as 100 nm [76] yielding the Ginzburg-Landau parameter κ ≈ 6.
B. Microwave and dc electrical measurements.
Combined broadband microwave and dc electrical measurements were done in a 4 He cryostat with magnetic field H directed perpendicular to the film surface. A custom-made cryogenic sample probe with coaxial cables was employed. The microwave and dc currents were superimposed and uncoupled by using two bias-tees mounted at the VNA ports. The dc voltage measurements were performed employing a Keithley Sourcemeter 2635B and an Agilent 34420A nanovoltmeter. The microwave signal, with a frequency between 30 MHz and 14 GHz, was generated and analyzed by a Keysight-Agilent E5071C vector network analyzer (VNA). In the transmission line, the microwave power level −6 dBm (0.25 mW) corresponds to the nominal microwave current ≃ 2.2 mA. The sample under study is a 50 Ohmmatched Nb coplanar waveguide (CPW), with the width of the central conductor w = 50 µm, the center-to-ground distance a = 20 µm and the length of the active area L = 1 mm. The coaxial cables were connected to the CPW via SMPB connectors spring-loaded to the goldplated contact pads sputtered on top of the Nb film using a shadow mask. With the estimate for the superconducting gap frequency f gap ≈ 700 GHz [77,78]  In the simulations, we solve numerically the Time Dependent Ginzburg-Landau (TDGL) equations [58,59] using the link variable method. To accomplish this, 2D and 3D TDGL solvers were developed following Ref. [79]. When testing the first simulations and comparing the results, both solvers gave very similar results on vortex motion and critical currents. Because of this, the 2D solver was used since it requires considerably less computational effort.
The simulations take place in a 360x100 cells domain of superconducting material. The size of the cells is scaled by the coherence length ξ 0 = ξ(0), which was experimentally estimated to be around ξ 0 = 17 nm. The lateral size of each cell is 0.5ξ 0 , so the simulated strip is about 3 × 0.8 µm 2 . In order to simulate a strip such as the one used in the experiments, we could either make the cells represent a much bigger area or use a domain with ∼ 10 5 cells on each side. The former would result in a simulation unable to resolve individual vortices, and the latter would take huge amounts of time and computational resources for a single simulation. Our solution to this was choosing a size such that we can both have a big number of vortices (more than 100 vortices fit easily in each simulation, as can be seen in Fig. 6) and resolve their individual shape and motion, while being able to perform several simulations in a reasonable time. Therefore, these simulations must be interpreted in qualitative terms to understand the effects of different dc and ac currents on superconductivity, and not as a quantitative prediction.
The simulation procedure is as following: An external dc magnetic field, uniform everywhere, is applied and populates the sample with vortices. An ac magnetic field (generated by the ac current) of different sign is applied at each boundary, of constant amplitude and frequency throughout each simulation. Finally, a dc magnetic field, of different signs at each boundary, is applied. Its magnitude is changed in steps of 0.005H c2 every 10 4 steps, with the time step defined as a x,y,z = 0.5ξ 0 being the lateral size of each cell, and time in units of the relaxation time for the order parameter, t 0 = π /96k B T c . For the smallest time step used (using κ = 6, thus dt = 0.00104t 0 ) this equals to 10 4 ×0.00104 = 10.4t 0 , giving enough time to reach a stationary state before this dc stray field was increased once more. This is similar to changing the applied dc current in the current-voltage curve. More details on the simulation procedure are provided in Ref. [80]. 644348 (MagIC). Research leading to these results was conducted within the framework of the COST Action CA16218 (NANOCOHYBRI) of the European Cooperation in Science and Technology. Support through the Frankfurt Center of Electron Microscopy (FCEM) is gratefully acknowledged.