Scientific Reports

Sci Rep. 7: 13740

Mobile fluxons as coherent probes of periodic pinning in superconductors

http://orcid.org/0000-0002-7895-8265Oleksandr V. Dobrovolskiy

Dobrovolskiy@Physik.uni-frankfurt.de

¹², http://orcid.org/0000-0001-7415-465XMichael Huth¹, Valerij A. Shklovskij², Ruslan V. Vovk²

1. 0000 0004 1936 9721grid.7839.5Physikalisches Institut, Goethe University, Frankfurt am Main, 60438 Germany

2. 0000 0004 0517 6080grid.18999.30Physics Department, V. Karazin Kharkiv National University, Kharkiv, 61022 Ukraine

© The Author(s) 2017

DOI: 10.1038/s41598-017-14232-z

Published online: 23 October 2017

Abstract

The interaction of (quasi)particles with a periodic potential arises in various domains of science and engineering, such as solid-state physics, chemical physics, and communication theory. An attractive test ground to investigate this interaction is represented by superconductors with artificial pinning sites, where magnetic flux quanta (Abrikosov vortices) interact with the pinning potential U(r) = U(r + R) induced by a nanostructure. At a combination of microwave and dc currents, fluxons act as mobile probes of U(r): The ac component shakes the fluxons in the vicinity of their equilibrium points which are unequivocally determined by the local pinning force counterbalanced by the Lorentz force induced by the dc current, linked to the curvature of U(r) which can then be used for a successful fitting of the voltage responses. A good correlation of the deduced dependences U(r) with the cross sections of the nanostructures points to that pinning is primarily caused by vortex length reduction. Our findings pave a new route to a non-destructive evaluation of periodic pinning in superconductor thin films. The approach should also apply to a broad class of systems whose evolution in time can be described by the coherent motion of (quasi)particles in a periodic potential.

Introduction

The interaction of (quasi)particles with a periodic potential arises in various domains of science and engineering, such as solid-state physics, chemical physics, and communication theory. Exemplary systems include Josephson junctions¹, superionic conductors², ring laser gyroscopes³, the dynamics of spin⁴ and charge-density⁵ waves, phase-locking loops⁶ in radioengineering, the motion of domain walls⁷, the magnetization dynamics of interacting spins⁸ and the diffusion of colloidal particles in periodic structures⁹. An interesting analogy can be found in superconductivity: Type II superconductors with periodic pinning sites substantiate an attractive test ground to investigate this interaction, thus allowing one to study dynamic properties of the vortex system and scrutinize its force-velocity (current-voltage) response.

The symmetry, shape and intensity of the pinning potential strongly affect the dynamics of Abrikosov vortices in superconductors, as reflected in both, the dc resistance and the microwave (mw) loss¹⁰. In this regard, tailoring periodic pinning sites has been revealed to be most efficient for controlling the resistive response of superconductors¹¹–¹³, dynamic mode-locking¹⁴, stimulation of superconductivity at microwave frequencies¹⁵, enhancing speed limits for the dynamics of Abrikosov vortices¹⁶–¹⁹ as well as enabling functionality of such fluxonic devices as rectifiers²⁰, transistors²¹, and high-frequency filters²². The availability of high-resolution nanofabrication tools has advanced the use of more and more sophisticated, asymmetric pinning geometries (vortex ratchets)¹¹,²⁰. For instance, in a vortex ratchet device (diode), the pinning site array may consist of elements of different sizes and shapes, such as triangles²⁰, grading circles²³, and arrow-shaped wedged cages²⁴. The coordinate dependence of the resulting pinning potential U(x, y) in these is rather complex, which is why for a quantitative interpretation of data an assumption is usually made for U(x, y) a priori. In the case of a simpler washboard pinning potential (WPP)²⁵, U(x, y) ≡ U(x) = U(x + a), where a is the period, an asymmetry can be induced by pre-defining the steepness of the left and right slopes of linear-extended pinning “sites” differently. The respective 3D washboard ratchet nanostructures can be realized, e.g., by milling of uniaxial nanogrooves by a focused ion beam (FIB) into the film²⁶, or via ferromagnetic decoration of the films by focused electron beam induced deposition (FEBID)²⁷. Since the geometrical appearance of an asymmetric washboard nanostructure does not reveal the asymmetry and the shape of the resulting WPP, non-destructive approaches for the determination of U(x, y) are needed.

The ideas for the determination of U(x, y) from experiment have been attracting attention of researchers for half a century²⁸–³². Early approaches for the determination of the distribution, density and strength of pinning sites included proposals to quantify the pinning force from measurements implying a small ripple magnetic field superimposed on a larger dc magnetic field²⁸ and from measurements with a small ac current²⁹. Recently, elastic properties of individual vortices were investigated by magnetic force microscopy³⁰. Further, scanning SQUID microscopy has been used to probe the dynamics and pinning of single vortices under combined dc and small ac drives, and the dependence of the elementary pinning force of multiple defects on the vortex displacement has been measured³¹. It has also been shown³² that vortices respond to local mechanical stress applied in the vicinity of a vortex thus allowing one to manipulate individual vortices without magnetic field or current. In all these experiments the samples contained randomly distributed pinning sites and either the vortex response in a non-coherent regime²⁸,²⁹ or that of an individual vortex³⁰–³² was probed.

An intriguing situation ensues in a washboard pinning landscape at a particular value of the magnetic field H m when each row of vortices is pinned at the bottom of the linearly-extended pinning sites (nanogrooves) and there are neither vacant nanogrooves nor vortices pinned between them. This field H m is the fundamental matching field at which the vortex lattice is commensurate with the pinning landscape. Computer simulations revealed³³ that in this case the vortex lattice has a crystalline structure, the effective intervortex interaction is cancelled, and each vortex experiences the same pinning force. An important consequence emerges from this: In the coherent regime, the moving ensemble of vortices in a WPP can be regarded as a moving vortex crystal, thus allowing one to interpret data employing single-vortex mechanistic models³⁴–³⁸.

Here, in contradistinction to previous works²⁸–³², we use the coherent vortex dynamics at the fundamental matching field to examine a theoretical mechanistic approach³⁹ for the determination of the coordinate dependence of a periodic pinning potential in superconductors under combined dc and mw current stimuli. Specifically, we evaluate the dc current-induced reduction of the depinning frequency in nanopatterned Nb films with different groove slope asymmetries and determine the dependences U(x) from ensemble-integrated microwave power absorption data. Further, using the pinning asymmetry parameters thus deduced, we augment the validity of the presented approach by a good fitting of mode-locking steps in the electrical voltage response in the presence of an ac drive to expressions derived within the framework of a stochastic model⁴⁰ of anisotropic pinning.

Results

We studied the vortex-state resistive response of Nb microstrips by combined microwave transmission spectroscopy and electrical voltage measurements, Fig. 1. The 500 μm long and 150 μm wide microstrips were patterned with arrays of 500 nm-spaced nanogrooves milled by FIB. The nanogrooves were parallel to the direction of the transport current density j which, in a perpendicular magnetic field B = μ ₀ H, exerted a transverse Lorentz force F L = j × B on vortices leading to their dissipative motion in the direction perpendicular to the grooves. The nanogroove arrays induced a periodic pinning potential of the washboard type U(x, y) ≡ U(x) = U(x + a). Sample S was patterned with nanogrooves having a symmetric cross-section, sample A1 with nanogrooves having a weak asymmetry of the groove slopes, and the control sample A2 with nanogrooves having a strong asymmetry of the slopes. Atomic force microscope images of the nanogrooves are presented further in the Discussion section. The asymmetric grooves are oriented in such a way that a positive dc bias makes the vortices to probe the gentle slope of the WPP: The vortices are shaken at the bottom of the pinning potential, which is gradually shifted in the negative x direction as the dc bias magnitude is increased. Samples S, A1, and A2 are 40 nm, 56 nm and 70 nm thick and have a superconducting transition temperature T c of 8.66 K, 8.85 K, and 8.94 K, respectively. The upper critical field for all samples at zero temperature is estimated as H c2(0) = H c2(T)/[1 − (T/T c)²] ≈ 1 T corresponding to a coherence length of ξ(0) = [Φ₀/2πH c2(0)]1/2 ≈ 17 nm. The diameter of the vortex core  ≃  2ξ at T = 0.3T c is by a factor of 6 smaller than the full width at half depth of the grooves in samples A1 and A2, and by a factor of 3 for sample S.

Figure 1View larger versionFigure 1. Sketch of the experimental setup: A combination of dc and mw current densities j + j mw is applied to a superconducting Nb microstrip. The microstrip contains an array of 500 nm-spaced nanogrooves milled by focused ion beam. The Lorentz force causes the vortices to move parallel to the x-axis. The measurable quantities are the dc voltage and the absolute value of the forward transmission coefficient S ₂₁ of the mw power at port 2 with respect to that at port 1.

“Ratchet window” in microwave power absorption

We investigated the frequency dependence of the mw power absorbed by vortices in Nb microstrips with symmetric and asymmetric nanogrooves at the fundamental matching field H m = 7.2 mT and T/T c = 0.3 for a series of dc bias values of both polarities. The data for the absorbed mw power were acquired in terms of the relative change of the absolute value of the forward transmission coefficient ΔS ₂₁(f) ≡ S ₂₁ − S 21ref, where S 21ref is the reference mw loss in the transmission line (all cables, connectors etc.) and ΔS ₂₁ is a measure for the mw loss due to vortex motion in the sample under study. Here all quantities are expressed in dB since it is the ratio of the transmission coefficients with and without the microstrip device which is a quantity of our interest.

In what follows we focus on the data acquired at a mw excitation power of P = −20 dBm, leaving an analysis of nonlinear effects observed at higher power levels beyond the scope of this work. In the order of magnitude, the amplitude of the vortex displacement at P = −20 dBm at a frequency of 1 GHz can be estimated as 10 nm on the basis of complementary fluctuation spectroscopy measurements on a series of narrow Nb microstrips fabricated from as-grown flat and grooved Nb films. The data at P = −30 dBm lead to essentially the same results but the data are more noisy. Contrarily, further effects emerge at high excitation power levels of P = −6 dBm and −3 dBm in the nonlinear regime which should be analyzed with an in-depth analysis of overheating effects and the associated with them quasiparticle escape from the vortex cores¹⁵. In this way, in the absence of the dc current, the ac-driven vortex dynamics is substantially in the linear regime (vortex displacements are small) and it transits into the dc-induced nonlinear regime as the dc bias is increasing so that the vortices are shaken while overcoming the WPP barriers.

Figure 2(a,b) presents the raw data ΔS ₂₁(f) for sample A1, while the respective data for samples S and A2 are reported in Supplementary. The arrangement of vortices at 7.2 mT in the pinning nanolandscape is shown in Fig. 2(c) for the assumed triangular vortex lattice with lattice parameter aΔ=(2Φ0/B3)1/2 and the matching condition aΔ=2a/3. In Fig. 2 ΔS ₂₁(f) has a smooth crossover from the weakly-dissipative regime at low frequencies (ΔS ₂₁(f) = 0 corresponds to the minimal mw loss) to the strongly-dissipative regime at high frequencies (ΔS ₂₁(f) ≈ −6.3 dB is the maximum mw excess loss due to vortices). The crossover takes place at the depinning frequency f d determined at the point where the phase difference between the viscous and pinning forces amounts to π/2³⁹ and corresponds to about −3 dB excess loss level in our measurements, which is indicated as a determination criterion for f d in Fig. 2(a). At low frequencies the vortices visit many potential wells (delocalized regime), whereas at high frequencies the vortices are shaken at the groove bottoms (localized regime). We note that the pinning force Fp = − ∇U(x) attains a maximum at the groove slopes, whereas F p ≈ 0 in the vicinity of the groove bottoms where U(x) is nearly constant. Accordingly, overcoming of the pinning barriers by vortices at low frequencies causes the pinning force to dominate the viscous force and the response is weakly dissipative, whereas in the high-frequency regime the viscous force dominates the pinning force, thus leading to a strongly dissipative response³⁴. The zero-bias depinning frequency f ₀ ≡ f d(j = 0) of samples S, A1 and A2 at T/T c = 0.3 and H m = 7.2 mT amounts to 5.72, 3.95, and 3.12 GHz, respectively.

Figure 2View larger versionFigure 2. Frequency dependence of the relative change of the absolute value of the forward transmission coefficient ΔS ₂₁(f) of sample A1 at positive (a) and negative (b) dc densities, as indicated. The arrow in (a) indicates the definition of the depinning frequency f d by the −3 dB criterion. (c) Vortex lattice configuration with lattice parameter aΔ=(2Φ0/H3)1/2 and the matching condition aΔ=2a/3 in a washboard nanolandscape with period a = 500 nm at the fundamental matching field H m = 7.2 mT.

For both dc bias polarities f d shifts towards low frequencies with increasing dc bias value, Fig. 2(a,b). Under dc polarity reversal, the magnitudes of the f d(j) shifts noticeably differ for sample A1, even more so for sample A2, but nearly coincide for sample S, see Supplementary. The reduction of the depinning frequency upon increasing the dc bias can be understood as a consequence of the effective lowering of one of the pinning potential wells due to its tilt caused by the dc bias⁴¹. The mechanistic consideration of a vortex as a particle leads to the conclusion that during an ac semiperiod, while the dc-tilted pinning potential well is broadening⁴¹, with increasing frequency f the vortex has no longer time ∝1/f to reach the groove slopes where the pinning forces dominate. Accordingly, the response becomes stronger dissipative already at lower frequencies, as compared with the zero-bias curve.

The difference in the magnitudes of the current-induced shifts of f d of samples A1 and A2 is best seen in Fig. 3 displaying contour plots of the difference signal ΔS21±=ΔS21(j+)−ΔS21(j−) which serves as a measure of the ratchet properties of the samples. The bright areas in the contour plots in Fig. 3(b,c) represent the “ratchet windows” in the microwave power absorption, i.e. the frequency - current range where the system exhibits rectifying properties, i.e. works as a diode. As expected, sample S with symmetric grooves does not exhibit ratchet properties, as the absolute value of the depinning current density does not change under dc bias polarity reversal in this sample. The “ratchet window” of samples A1 and A2 is bound by the gentle-slope jd+ and the steep-slope jd− depinning current densities (extrapolated to f → 0 corresponding to the dc regime), as well as by the value of the zero-bias depinning frequency f d(j = 0), as indicated in Fig. 3. The depinning current density j d amounts to 0.75 MA/cm² for both CVC branches of sample S, 0.91 MA/cm² and 0.7 MA/cm² for the steep-slope and the gentle-slope direction of the vortex motion in sample A1, and 1.25 MA/cm² and 0.52 MA/cm² for the steep-slope and the gentle-slope direction of the vortex motion in sample A2, respectively. The depinning current densities jd± are determined by the 10 μV/cm electric field strength criterion from the current-voltage curves (CVCs) discussed next.

Figure 3View larger versionFigure 3. Difference signals ΔS21±=ΔS21(j+)−ΔS21(j−) for sample S (a), A1 (b), and A2 (c) at T = 0.3T c and the fundamental matching field H m = 7.2 mT. The arrows indicate the zero-bias depinning frequencies f d(j = 0) and the depinning current densities jd+ and jd− for the motion of vortices against the gentle and steep slope of the grooves, respectively. The dashed lines are fits of the general form f d/f d(j = 0) = [1 − (j/j d)k/l]m/n, with exponents k, l, m, n indicated in Fig. 6.

Mode-locking steps in the current-voltage curves

Figure 4(a–c) displays the dc electrical voltage as a function of the normalized dc density and the ac amplitude at T = 0.98T c and the fundamental matching field H m = 7.2 mT. In all panels f = 0.3f d, with f d amounting to 444, 306, and 242 MHz, for samples S, A1, and A2, respectively. In the absence of microwave excitation, the CVC of Sample S in Fig. 4(a) is symmetric, while the CVCs of Sample A1 and A2 demonstrate two different absolute values of the depinning current density jd+>|jd−| for the positive and the negative branch. The addition of the microwave stimulus leads to the appearance of Shapiro steps in the CVC. The steps occur at voltages⁴² V = nV ₀ ≡ nNΦ₀ f, where n is an integer, N is the number of vortex rows between the voltage leads, f is the microwave frequency, and Φ₀ = 2.07 × 10⁻¹⁵ Vs is the magnetic flux quantum. The steps in the CVCs, which are best seen in the conventional current-voltage representation in Fig. 5(a), arise when the hopping distance of Abrikosov vortices during one ac halfwave coincides with one or a multiple of the nanostructure period. These interference steps are a fingerprint of the coherent vortex dynamics. They remain visible in the field range from 7 to 7.5 mT and disappear as H is tuned further away from the fundamental matching field H m = 7.2 mT.

Figure 4View larger versionFigure 4. Electrical dc voltage as a function of the normalized dc current density and the microwave amplitude for all samples at T = 0.98T c and H m = 7.2 mT. (a–c) Experimental data. (d–f) Theoretical modeling results, as detailed in Discussion. The color scale is the same in all panels with the electrical voltage normalized to the first Shapiro step voltage V ₀. Vertical dashed lines mark mode-locking fringes. Arrows depict the dc loading capability of the vortex ratchets, as deduced from crossing the straight lines j d(j mw) at the point where the internal asymmetry of the pinning potential is effectively compensated by the extrinsic asymmetry of the WPP caused by the dc bias.

Figure 5View larger versionFigure 5. (a) Shapiro steps in the CVC of microstrips S and A2 at jmw/jdmw=1.83, T = 0.98T c, and the fundamental matching field H m = 7.2 mT. Symbols: Experiment. Solid lines: Fits to Eq. (6). (b) Mode-locking “fringes” in the dc voltage with increasing mw amplitude at T = 0.98T c, H m = 7.2 mT, and a series of dc densities as indicated.

We note that temperatures close to T c are fortunate for the observation of Shapiro steps which are best seen in the flux flow regime at j > j d. This is because of the depinning current decreasing with increasing temperature and the flux flow setting on at smaller dc densities j d < j < j * thus allowing one to operate in an extended current range where overheating effects are negligible. Here j * is the current density corresponding to an abrupt transition of the sample into the normal state due to the Larkin-Ovchinnikov instability¹⁶,¹⁸,¹⁹,⁴³–⁴⁸.

In Fig. 4(a–c) one sees that j d nearly linearly decreases with increasing microwave amplitude and at some point the lines j d(j mw) for different polarities cross, thus allowing one to deduce a dc bias value for which an effective “symmetrization” of the WPP is predicted⁴⁰: At this dc bias value the internal asymmetry of the pinning potential is effectively compensated by the extrinsic asymmetry of the WPP caused by the dc bias. Here we refer to Fig. 2 in Supplementary for a sketch. Quantitatively, these dc values characterize the loading capability of the ratchet⁴⁹ and amount to j/j d ≈ 0.1 for sample A1 and j/j d ≈ 0.25 for sample A2. Obviously, the lines jd±(jmw) cross at j/j d ≈ 0 for sample S as this sample has symmetric grooves. As the ac amplitude increases, the dc voltage for all samples exhibits mode-locking peculiarities (“fringes”). Up to four lowest-order fringes for sample A2 are exemplified in Fig. 5(b) for a series of sub-depinning dc biases of both polarities. Both, Shapiro steps and mode-locking fringes were theoretically predicted earlier for symmetric³⁸ and asymmetric⁴⁰ WPPs. To examine whether the voltage responses can be described theoretically, we determine the coordinate dependence U(x) of the pinning potential from the microwave power absorption and use the deduced expressions for U(x) for modeling of the voltage responses in Fig. 4(d–f), as detailed next.

Discussion

Determination of the coordinate dependence of the pinning potential

Figure 6(a–c) displays the dependences f d/f ₀ versus j/j d deduced from the raw experimental data ΔS ₂₁(f), as shown in Fig. 2 for a series of dc bias values of both polarities for all samples. The data acquired for positive dc biases are shown by solid red symbols, while those for negative dc biases by open blue ones. The relative uncertainty in the determination of the depinning frequency does not exceed 7% at j → j d and it is smaller than 3% at dc current densities j < 0.8j d. This is indicated by the error bars in Fig. 6(a–c), which become smaller than the symbol size at j < 0.8j d. The larger error in the determination of f d at j → j d is associated with a smeared functional shape of the curves ΔS ₂₁(f) in Fig. 2 upon vanish of the WPP barriers due to their tilt by the dc current.

Figure 6View larger versionFigure 6. (a–c) Reduction of the depinning frequency upon increasing the dc density as deduced from the mw power absorption data. The experimental data for the positive (●) and the negative (□) dc polarity are approximated by fits (solid lines of respective color) of the general form f d/f ₀ = [1 − (j/j d)k/l]m/n, with the exponents k, l, m, n as indicated. (d–f) Atomic force microscope images of the WPP landscapes milled by FIB on the surface of samples S, A1, and A2. When subjected to a combination of microwave and dc currents, fluxons act as movable coherent “sensors” probing the curvature of U(x): The ac component shakes the fluxons in the vicinity of their equilibrium points which are unequivocally determined by the local pinning force counterbalanced by the driving (Lorentz) force induced by the dc bias current. A positive dc bias makes the vortices to “probe” the gentle groove slope, whereas the steep groove slope is probed by vortices when they are subject to a negative dc bias. (g–i) The coordinate dependences of the pinning potentials U(x) deduced numerically using the approach outlined in Supplementary (solid lines) are plotted together with fits (dashed lines) to Eq. (2) and the cross-sectional AFM line scans across the grooves (symbols). The light red areas in panels (g–i) correspond to the left-hand WPP groove slopes in panels (a–c) probed at positive dc bias values (red solid symbols in panels (a–c)). The light blue areas in panels (g–i) correspond to the right-hand WPP groove slopes in panels (a–c) probed at negative dc bias values (blue open symbols in panels (a–c)).

To interpret the dc-induced reduction of the depinning frequency f d/f ₀ in Fig. 6(a–c) as a function of the dc current density j/j d at 0.3T c for the coherent vortex at the fundamental matching field H m = 7.2 mT, we employ a single-vortex theoretical approach outlined in Supplementary. This approach is based on the Gittleman-Rosenblum model³⁴ generalized for the presence of a dc current bias³⁹,⁵⁰,⁵¹. The procedure of determination of the coordinate dependence of the pinning potential from the microwave power absorbed by vortices requires to approximate the reduction of the depinning frequency for both dc polarities by some expression. The experimental data for all our samples fit to expressions of the universal “mean-field” type

1. fd/f0=[1−(j/jd)k/l]m/n[Image: 41598_2017_14232_Article_Equ1.gif]

with the exponents k/l and m/n labeled at the respective fit curves in Fig. 6(a–c). For consistency, we also plot the respective fit curves in the contour plots in Fig. 3 and note that they nicely encage the “ratchet windows” quantitatively.

The employment of the procedure outlined in Supplementary to the fit curves in Fig. 6(a–c) results in the U(x) curves shown in Fig. 6(g–i) by solid lines. The curves U(x) in Fig. 6(g–i) are plotted for one period of the pinning potential scaled to its depth U ₀ and period a. The line thickness reflects the uncertainty in the determination of U(x), which does not exceed 10% and 5% on the dome and at the bottom of the potential, respectively. Analytically, the dependences U(x) deduced from the mw data can be approximated by expressions of the general form

2. U(x)/U0=[(1−cos2πx/a)+e(1−sin4πx/a)/2]/2[Image: 41598_2017_14232_Article_Equ2.gif]

with e being the asymmetry parameter and amounting to 0, 0.13 and 0.56 for samples S, A1, and A2, respectively.

Figure 6(g–i) displays the cross-sections of the nanogrooves measured by atomic force microscopy (AFM) in comparison with the dependences U(x) deduced from the mw data and the respective approximations by Eq. (2). Obviously, the deduced coordinate dependence of the pinning potential U(x) in Sample A1 corresponds quite well to the coordinate dependence of the film thickness d(x) modulated by the nanogrooves, as measured by AFM. The AFM profile for Sample A2 can also be fitted satisfactorily to the same expression as U(x). At the same time, U(x) for sample S reproduces the AFM profile qualitatively, but does not follow it quantitatively. The modulated thickness d(x) in Sample S can be fitted to

3. d(x)=[(1−cos2πx)+0.56(1−cos4πx)/2]/2.[Image: 41598_2017_14232_Article_Equ3.gif]

The stronger discrepancy between the dependences U(x) and d(x) in sample S may be caused by a factor of about two smaller groove width in Sample S as compared to Samples A1 and A2. In all, we draw the conclusion that the employed procedure³⁹ is not only sensitive to the asymmetry of the pinning potential, but it also reveals a correlation between U(x) and d(x) thus indicating that pinning in the two thicker samples is primarily caused by vortex length reduction. This can be understood as the vortex energy ε(x) is proportional to its length L(x) determined by the film thickness ε(x) ∝ L(x) = d(x).

Modeling of electrical resistance responses

To augment the validity of the employed approach for the determinations of the coordinate dependences U(x) of the pinning potential from the microwave power absorption at T = 0.3T c we employ a stochastic model of uniaxial periodic pinning³⁷,³⁸,⁴⁰ to fit the peculiarities observed in electrical voltage responses at noticeably higher temperatures T = 0.98T c. The theoretical treatment relies upon the Langevin equation for a vortex moving with velocity v in a magnetic field B = z B, where B ≡ |B|, and z is the unit vector in the z direction

4. ηv = FL + Fp + Fth, [Image: 41598_2017_14232_Article_Equ4.gif]

where F L = (Φ₀/c)j × z is the Lorentz force, Φ₀ is the magnetic flux quantum, and c is the speed of light. In Eq. (4) j ≡ j(t) = j + j mw cos ωt, where j and j mw are the dc and ac current density amplitudes and ω is the angular frequency. Fp = − ∇Up(x) is the anisotropic pinning force, where U p(x) is a ratchet WPP. F th is the thermal fluctuation force represented by Gaussian white noise and η is the vortex viscosity.

The ratchet WPP is modeled by

5. Up(x)=(Up/2)[1−coskx+e(1−sin2kx)/2],[Image: 41598_2017_14232_Article_Equ5.gif]

where k = 2π/a with a being the period and U p the depth of the WPP. In Eq. (5) e is the asymmetry parameter allowing for tuning the asymmetry strength. It is this parameter which must be determined from experiment.

The asymmetry parameter deduced from the experimental data amounts to e = 0 for Sample S, e = 0.13 for Sample A1, and e = 0.56 for Sample A2. These values of e are used in the simulations of the electric field response given by⁴⁰

6. E = ν(j, jmw, f, t)j.[Image: 41598_2017_14232_Article_Equ6.gif]

Here ν is the (j, j mw, f, t)-dependent effective nonlinear mobility of the vortex under the influence of the Lorentz force and it is expressed in terms of matrix continued fractions⁴⁰. In Eq. (6), the electric field is scaled to the electric field at the first Shapiro step, while the dc density j and the mw density amplitude j mw to the depinning current density, frequency f to the depinning frequency, coordinate x to the WPP period, and temperature t to the WPP depth.

We use the mean-square parameter jd=jd+jd− for the presentation of the data in dimensionless form allowing for a direct comparison of experiment with theory in Fig. 4. The voltage responses are calculated by Eq. (6) for a series of values of the asymmetry parameter e at the reduced temperature t = T/U p = 0.002 which corresponds to the pinning activation energy U p of about 5000 K⁵² and the temperature T = 0.98T c at which the experiment is conducted.

We now proceed to a deeper analysis of the mode-locking fringes reported in Fig. 5(b). The appearance of these voltage peaks (for the positive dc polarity) and dips (for the negative dc polarity) can be explained by higher-order ratchet effects (labeled with “1” to “4” in the plot) as follows: We first consider the curve V/V0(jmw/jdmw) in Fig. 5(b) at j/j d = 0.05, that is in the limit of very small dc biases. In the conventional ratchet effect resulting in the appearance of the “1” voltage peak, at a given ac amplitude jmw>˜jdgentlemw, a vortex overcomes the gentle slope of the WPP during one half of the ac period. However, the amplitude j mw is yet smaller than the strong-slope depinning current density and the vortex can not return into its original WPP well. This results in a net motion of the vortex by one WPP period and the associated rectified voltage. As the ac amplitude increases, the vortex can overcome the steep barrier of the WPP and the net motion disappears thus resulting in an almost full suppression of the rectified voltage. With a further increase of the ac amplitude the vortex can consequentially overcome two gentle barriers of WPP, and two options appear for its backward motion depending on the strength of the WPP asymmetry and the ac frequency, please refer to Fig. 2 of Supplementary: (i) The vortex can only overcome one barrier in the backward direction, and after one ac period it appears in the WPP well in one period away from the original WPP well. (ii) The vortex can not overcome the barrier in the backward direction and it remains in the WPP well in two periods away from the original WPP well. Obviously, with a further increase of the ac amplitude the “splitting” of ways increases for the vortex in what WPP well it will appear after one ac period and this leads to a smearing of the higher-order ratchet peaks in the dc voltage.

The curve at j = 0.25 in Fig. 5(b) corresponds to the particular interesting case of effective compensation of the intrinsic WPP asymmetry due to the difference in the steepness of the groove slopes by the external asymmetry of the potential induced by the dc bias. In consequence of this effective “compensation of the anisotropies” the “floor” of the rectified voltage can be adjusted to almost zero which is why the loading capability of the ratchet amounts to about j/j d = 0.25 for sample A2. In all, the simulation results allow us not only to explain the main features observed in the experiment, but we find a very good quantitative agreement between the experimental data and theoretical modeling.

In conclusion, we have presented an approach allowing for the determination of a periodic pinning potential from the microwave power absorbed by vortices under dc bias reversal in superconductors with periodic pinning at low temperatures. For thicker films, the deduced coordinate dependences of the washboard pinning potentials U(x) largely mimic the coordinate dependences of the film thickness, as visible from cross-sectional AFM line scans. The presented procedure allowed us to directly determine the asymmetry of the pinning potential and to further use it for modeling of voltage responses at T = 0.98T c. Our findings pave a new route to the non-destructive evaluation of periodic pinning in superconductor thin films. Moreover, the reported approach should also apply to a broad class of systems whose evolution in time can be described by the coherent motion of (quasi)particles in a periodic potential.

Methods

Film growth and characterization

The 150 × 500 μm² microstrips 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⁻⁸ mbar range. In the sputtering process the substrate temperature was 850 °C, the Ar pressure 4 × 10⁻³ mbar, and the growth rate was about 1 nm/s. X-ray diffraction measurements revealed the (110) orientation of the films⁵³. The epitaxy of the films has been 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 deduced from AFM scans in the range 1 μm × 1 μm. The background pinning in the epitaxial (110) Nb films is very weak as compared to the anisotropic pinning induced by nanogrooves milled by FIB. The parameters of the samples are compiled in Table 1.

See full tableTable 1. Main parameters of the investigated Nb microstrips.

Fabrication of nanogrooves

Patterning of the samples was done in a high-resolution dual-beam scanning electron microscope (FEI, Nova Nanolab 600). In the patterning process, the asymmetry of the groove slopes was achieved by defining the grooves in the FIB bitmap file for sample S as a single line for the beam to pass, whereas a step-wise increasing number of FIB beam passes was assigned to each groove defined as a 3-step and a 5-step “stair” for samples A1 and A2, respectively. Due to blurring effects, the symmetric grooves in Sample S have rounded corners while smoothed straight slopes resulted instead of the “stairs” in samples A1 and A2. For all samples the beam parameters were 30 kV/50 pA, 1 μs dwell time and 50 nm pitch. The grooves are parallel to the microstrip edges with a misalignment of less than 0.2°. The microstrip width is an integer multiple number (N = 300) of the nanopattern period to prevent possible ratchet effects due to the edge barrier asymmetry¹³.

Atomic force microscopy

A Nanosurf easyScan 2 atomic force microscope (AFM) under ambient conditions in non-contact, dynamic force mode was used. The cantilever tip was shaped like a polygon-based pyramid, with a tip radius of less than 7 nm (Nanosensors PPP-NCLR). Convolution effects due to the finite tip radius can be neglected, as is corroborated by the invariance of the AFM images taken with the cantilever scanning at different angles with respect to the grooves. The cross-sectional AFM profiles are the results of averaging of 250 line scans acquired for a scanning field of 500 × 500 nm².

Microwave spectroscopy

Combined broadband mw and dc electrical measurements were done in a ⁴He cryostat with magnetic field H directed perpendicular to the film surface. A custom-made cryogenic sample probe with coaxial cables was employed⁵⁴. The mw signal was generated and analyzed by an Agilent E5071C vector network analyzer (VNA). The mw and dc signals were superimposed and uncoupled by using two bias-tees mounted at the VNA ports. The VNA operated in the frequency sweep mode, with 1548 frequency points scanned with an exponentially growing increment between 300 KHz and 14 GHz. For all frequencies f, the mw excitation power at the sample was P = −20 dBm (10 μW) kept by the VNA in accordance with the pre-saved calibration data for S 21ref(f, T). At small fields, as investigated by us, S 21ref does not depend on the magnetic field value H.

Modeling of voltage responses

The infinite matrix continued fractions for the calculation of the voltage responses by Eq. (29) of ref.⁴⁰ were approximated by matrix continued fractions of finite order. This has been done by putting Q m = 0 at some m = M, whereas the dimension of the submatrices Q m and the vectors C m was confined to some finite number K. Both M and K depend on temperature and mw density and on the number of harmonics to be taken into account. These numbers were chosen as K = 51 and M = 500 for the reliable calculation of the components Fk1(ω) for up to 10 harmonics, for j mw up to 5, and for t ≡ T/U ₀ = 0.002, respectively. This ensures a calculation accuracy of at least three digits. The Hall effect was neglected in the calculations.

Supplementary Material

Open In Web Browser

Footnotes

Electronic supplementary material

Supplementary information accompanies this paper at 10.1038/s41598-017-14232-z.

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

Acknowledgements

Roland Sachser is acknowledged for automating the data acquisition and his help with nanopatterning. This work was supported through DFG project DO1511/3-1 and conducted within the framework of the COST Action CA16218 (NANOCOHYBRI) of the European Cooperation in Science and Technology. Further, funding from the European Commission in the framework of the program Marie Sklodowska-Curie Actions - Research and Innovation Staff Exchange (MSCA-RISE) under Grant Agreement No. 644348 (MagIC) is acknowledged.

Author Contributions

O.V.D. conducted the experiment and wrote the manuscript. O.V.D. and M.H. designated the samples. V.A.S. suggested the concept of the experiment and analyzed the data. R.V.V. took part in the discussion of the results and analyzed the electrical voltage responses. M.H., V.A.S., and R.V.V. participated in the manuscript writing.

Notes

Competing Interests

The authors declare that they have no competing interests.

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Barone, A. & Patterno, G. Physics and Applications of the Josephson Effect (John Wiley & Sons, New York, 1982).

Fulde P, Pietronero L, Schneider WR, Strässler S, authors. Problem of brownian motion in a periodic potential. Phys. Rev. Lett. 1975;35:1776–1779. [Cross Ref]

Chow WW, et al., authors. The ring laser gyro. Rev. Mod. Phys. 1985;57:61–104. [Cross Ref]

Barthel E, et al., authors. Conduction noise and motional narrowing of the nuclear magnetic resonance line in sliding spin-density waves. Phys. Rev. Lett. 1993;71:2825–2828. [Cross Ref] [PubMed]

Zybtsev SG, Pokrovskii VY, authors. Coherence of a sliding charge-density wave stimulated by asynchronous rf irradiation. Phys. Rev. B. 2013;88:125144 [Cross Ref]

Risken, H. The Fokker-Planck Equation (Springer, Berlin, 1989).

Pérez-Junquera A, et al., authors. Crossed-ratchet effects for magnetic domain wall motion. Phys. Rev. Lett. 2008;100:037203 [Cross Ref] [PubMed]

Titov SV, Kachkachi H, Kalmykov YP, Coffey WT, authors. Magnetization dynamics of two interacting spins in an external magnetic field. Phys. Rev. B. 2005;72:134425 [Cross Ref]

Evstigneev M, et al., authors. Diffusion of colloidal particles in a tilted periodic potential: Theory versus experiment. Phys. Rev. E. 2008;77:041107 [Cross Ref]

Brandt, E. H. The flux-line lattice in superconductors. Rep. Progr. Phys. 58, 1465–1594, http://stacks.iop.org/0034-4885/58/i=11/a=003 (1995).

Plourde, B. L. T. Nanostructured superconductors with asymmetric pinning potentials: Vortex ratchets. IEEE Trans. Appl. Supercond. 19, 3698–3714, http://ieeexplore.ieee.org/xpl/articleDetails.jsp?arnumber=5232856, 10.1109/TASC.2009.2028873 (2009).

Cuadra-Solis P-d-J, et al., authors. Observation of commensurability effects in a patterned thin superconducting pb film using microwave reflection spectrometry. Phys. Rev. B. 2014;89:054517 [Cross Ref]

Pryadun, V. V., Sierra, J., Aliev, F. G., Golubovic, D. S. & Moshchalkov, V. V. Plain superconducting films as magnetic field tunable two-dimensional rectifiers. Appl. Phys. Lett. 88, 10.1063/1.2171788 (2006).

Jelic, Z. L., Milosevic, M. V., Van de Vondel, J. & Silhanek, A. V. Stroboscopic phenomena in superconductors with dynamic pinning landscape. Sci. Rep. 5, 14604 EP–, 10.1038/srep14604, Article (2015).

Lara A, Aliev FG, Silhanek AV, Moshchalkov VV, authors. Microwave-stimulated superconductivity due to presence of vortices. Sci. Rep. 2015;5:9187 [Cross Ref] [PubMed]

Silhanek, A. V. et al. Influence of artificial pinning on vortex lattice instability in superconducting films. New J. Phys. 14, 053006, http://stacks.iop.org/1367-2630/14/i=5/a=053006 (2012).

Grimaldi G, et al., authors. Speed limit to the abrikosov lattice in mesoscopic superconductors. Phys. Rev. B. 2015;92:024513 [Cross Ref]

Shklovskij VA, Nazipova AP, Dobrovolskiy OV, authors. Pinning effects on self-heating and flux-flow instability in superconducting films near Tc. Phys. Rev. B. 2017;95:184517 [Cross Ref]

Dobrovolskiy, O. V. et al. Pinning effects on flux flow instability in epitaxial nb thin films. Supercond. Sci. Technol. 30, 085002, http://stacks.iop.org/0953-2048/30/i=8/a=085002 (2017).

Villegas, J. E. et al. A superconducting reversible rectifier that controls the motion of magnetic flux quanta. Science302, 1188–1191, http://www.sciencemag.org/content/302/5648/1188.abstract, 10.1126/science.1090390 (2003).

Vlasko-Vlasov, V. K., Colauto, F., Benseman, T., Rosenmann, D. & Kwok, W.-K. Triode for magnetic flux quanta. Sci. Rep. 6, 36847 EP–, 10.1038/srep36847, Article (2016).

Dobrovolskiy, O. V. & Huth, M. Dual cut-off direct current-tunable microwave low-pass filter on superconducting nb microstrips with asymmetric nanogrooves. Appl. Phys. Lett. 106, 142601–1–5, 10.1063/1.4917229 (2015).

Gillijns W, Silhanek AV, Moshchalkov VV, Reichhardt CJO, Reichhardt C, authors. Origin of reversed vortex ratchet motion. Phys. Rev. Lett. 2007;99:247002 [Cross Ref] [PubMed]

Togawa Y, et al., authors. Direct observation of rectified motion of vortices in a niobium superconductor. Phys. Rev. Lett. 2005;95:087002 [Cross Ref] [PubMed]

Dobrovolskiy, O. V. Abrikosov fluxonics in washboard nanolandscapes. Physica C533, 80–90, http://www.sciencedirect.com/science/article/pii/S092145341630096X, 10.1016/j.physc.2016.07.008 (2017).

Dobrovolskiy, O. V., Huth, M. & Shklovskij, V. A. Alternating current-driven microwave loss modulation in a fluxonic metamaterial. Appl. Phys. Lett. 107, 162603–1–5, 10.1063/1.4934487 (2015).

Dobrovolskiy, O. V., Huth, M. & Shklovskij, V. A. Anisotropic magnetoresistive response in thin Nb films decorated by an array of Co stripes. Supercond. Sci. Technol. 23, 125014–1–5, 10.1088/0953-2048/23/12/125014 (2010).

Campbell AM, Evetts JE, authors. Flux vortices and transport currents in type ii superconductors. Adv. Phys. 1972;21:199–428. [Cross Ref]

Lowell, J. Investigation of pinning sites in superconductors by means of alternating currents. ii. experiments. J. Phys. F: Metal Phys. 2, 559, http://stacks.iop.org/0305-4608/2/i=3/a=023 (1972).

Auslaender OM, et al., authors. Mechanics of individual isolated vortices in a cuprate superconductor. Nat Phys. 2009;5:35–39. [Cross Ref]

Embon L, et al., authors. Probing dynamics and pinning of single vortices in superconductors at nanometer scales. Sci Rep. 2015;5:7598 [Cross Ref] [PubMed]

Kremen, A. et al. Mechanical control of individual superconducting vortices. Nano Letters0, null, 10.1021/acs.nanolett.5b04444, PMID: 26836018 (2016).

Lu Q, Reichhardt CJO, Reichhardt C, authors. Reversible vortex ratchet effects and ordering in superconductors with simple asymmetric potential arrays. Phys. Rev. B. 2007;75:054502 [Cross Ref]

Gittleman JI, Rosenblum B, authors. Radio-frequency resistance in the mixed state for subcritical currents. Phys. Rev. Lett. 1966;16:734–736. [Cross Ref]

Coffey MW, Clem JR, authors. Unified theory of effects of vortex pinning and flux creep upon the rf surface impedance of type-II superconductors. Phys. Rev. Lett. 1991;67:386–389. [Cross Ref] [PubMed]

Pompeo, N. & Silva, E. Reliable determination of vortex parameters from measurements of the microwave complex resistivity. Phys. Rev. B78, 094503–1–10, 10.1103/PhysRevB.78.094503 (2008).

Shklovskij, V. A. & Dobrovolskiy, O. V. AC-driven vortices and the Hall effect in a superconductor with a tilted washboard pinning potential. Phys. Rev. B78, 104526–1–12, 10.1103/PhysRevB.78.104526 (2008).

Shklovskij, V. A. & Dobrovolskiy, O. V. Frequency-dependent ratchet effect in superconducting films with a tilted washboard pinning potential. Phys. Rev. B84, 054515–1–12, 10.1103/PhysRevB.84.054515 (2011).

Shklovskij, V. A. & Dobrovolskiy, O. V. Microwave Absorption by Vortices in Superconductors with a Washboard Pinning Potential, chap. 11, 263–288 (InTech, Rijeka, 2012).

Shklovskij, V. A., Sosedkin, V. V. & Dobrovolskiy, O. V. Vortex ratchet reversal in an asymmetric washboard pinning potential subject to combined dc and ac stimuli. J. Phys.: Cond. Matt. 26, 025703, http://stacks.iop.org/0953-8984/26/i=2/a=025703 (2014).

Shklovskij VA, Hop DTB, authors. The Hall effect and microwave absorption by vortices in an anisotropic superconductor with a periodic pinning potential. Low Temp. Phys. 2010;36:71–80. [Cross Ref]

Fiory AT, author. Quantum interference effects of a moving vortex lattice in al films. Phys. Rev. Lett. 1971;27:501–503. [Cross Ref]

Larkin, A. I. & Ovchinnikov, Y. N. Nonlinear conductivity of superconductors in the mixed state. J. Exp. Theor. Phys. 41, 960–965, http://www.jetp.ac.ru/cgi-bin/index/e/41/5/p960?a=list (1975).

Musienko, L. E., Dmitrenko, I. M. & Volotskaya, V. G. JETP Lett. 31, 567 (1980).

Klein W, Huebener RP, Gauss S, Parisi J, authors. Nonlinearity in the flux-flow behavior of thin-film superconductors. J. Low Temp. Phys. 1985;61:413–432. [Cross Ref]

Volotskaya, V. G., Dmitrenko, I. M., Koretskaya, O. A. & Musienko, L. E. Fiz. Nizk. Temp. 18, 973 (1992).

Peroz, C. & Villard, C. Flux flow properties of niobium thin films in clean and dirty superconducting limits. Phys. Rev. B72, 014515–1–6, 10.1103/PhysRevB.72.014515 (2005).

Leo A, et al., authors. Competition between intrinsic and extrinsic effects in the quenching of the superconducting state in fe(se,te) thin films. Phys. Rev. B. 2016;93:054503 [Cross Ref]

Knufinke, M. et al. Deterministic josephson vortex ratchet with a load. Phys. Rev. E85, 011122–1–9, 10.1103/PhysRevE.85.011122 (2012).

Shklovskij, V. A. Determination of coordinate dependence of the washboard pinning potential from the dynamic experiment with vortices. In Procedings of the Fifth International Conference on Mathematical Modeling and Computer Simulation of Materials Technologies MMT-2008, Ariel, Israel (2008).

Shklovskij VA, Dobrovolskiy OV, authors. Determination of coordinate dependence of a pinning potential from a microwave experiment with vortices. Low Temp. Phys. 2013;39:120–124. [Cross Ref]

Dobrovolskiy, O. V., Begun, E., Huth, M. & Shklovskij, V. A. Electrical transport and pinning properties of Nb thin films patterned with focused ion beam-milled washboard nanostructures. New J. Phys. 14, 113027–1–27, http://stacks.iop.org/1367-2630/14/i=11/a=113027 (2012).

Dobrovolskiy, O. V. & Huth, M. Crossover from dirty to clean superconducting limit in dc magnetron-sputtered thin Nb films. Thin Solid Films520, 5985–5990, http://www.sciencedirect.com/science/article/pii/S0040609012005718, 10.1016/j.tsf.2012.04.083 (2012).

Dobrovolskiy, O. V., Franke, J. & Huth, M. ⁴He sample probe for combined microwave and dc transport measurements. Meas. Sci. Technol. 26, 035502, http://stacks.iop.org/0957-0233/26/i=3/a=035502 (2015).

[Back]

Figure 1.

Sketch of the experimental setup: A combination of dc and mw current densities j + j mw is applied to a superconducting Nb microstrip. The microstrip contains an array of 500 nm-spaced nanogrooves milled by focused ion beam. The Lorentz force causes the vortices to move parallel to the x-axis. The measurable quantities are the dc voltage and the absolute value of the forward transmission coefficient S ₂₁ of the mw power at port 2 with respect to that at port 1.

[Back]

[Back]

Figure 2.

Frequency dependence of the relative change of the absolute value of the forward transmission coefficient ΔS ₂₁(f) of sample A1 at positive (a) and negative (b) dc densities, as indicated. The arrow in (a) indicates the definition of the depinning frequency f d by the −3 dB criterion. (c) Vortex lattice configuration with lattice parameter aΔ=(2Φ0/H3)1/2 and the matching condition aΔ=2a/3 in a washboard nanolandscape with period a = 500 nm at the fundamental matching field H m = 7.2 mT.

[Back]

[Back]

Figure 3.

Difference signals ΔS21±=ΔS21(j+)−ΔS21(j−) for sample S (a), A1 (b), and A2 (c) at T = 0.3T c and the fundamental matching field H m = 7.2 mT. The arrows indicate the zero-bias depinning frequencies f d(j = 0) and the depinning current densities jd+ and jd− for the motion of vortices against the gentle and steep slope of the grooves, respectively. The dashed lines are fits of the general form f d/f d(j = 0) = [1 − (j/j d)k/l]m/n, with exponents k, l, m, n indicated in Fig. 6.

[Back]

[Back]

Figure 4.

Electrical dc voltage as a function of the normalized dc current density and the microwave amplitude for all samples at T = 0.98T c and H m = 7.2 mT. (a–c) Experimental data. (d–f) Theoretical modeling results, as detailed in Discussion. The color scale is the same in all panels with the electrical voltage normalized to the first Shapiro step voltage V ₀. Vertical dashed lines mark mode-locking fringes. Arrows depict the dc loading capability of the vortex ratchets, as deduced from crossing the straight lines j d(j mw) at the point where the internal asymmetry of the pinning potential is effectively compensated by the extrinsic asymmetry of the WPP caused by the dc bias.

[Back]

[Back]

Figure 5.

(a) Shapiro steps in the CVC of microstrips S and A2 at jmw/jdmw=1.83, T = 0.98T c, and the fundamental matching field H m = 7.2 mT. Symbols: Experiment. Solid lines: Fits to Eq. (6). (b) Mode-locking “fringes” in the dc voltage with increasing mw amplitude at T = 0.98T c, H m = 7.2 mT, and a series of dc densities as indicated.

[Back]

[Back]

Figure 6.

(a–c) Reduction of the depinning frequency upon increasing the dc density as deduced from the mw power absorption data. The experimental data for the positive (●) and the negative (□) dc polarity are approximated by fits (solid lines of respective color) of the general form f d/f ₀ = [1 − (j/j d)k/l]m/n, with the exponents k, l, m, n as indicated. (d–f) Atomic force microscope images of the WPP landscapes milled by FIB on the surface of samples S, A1, and A2. When subjected to a combination of microwave and dc currents, fluxons act as movable coherent “sensors” probing the curvature of U(x): The ac component shakes the fluxons in the vicinity of their equilibrium points which are unequivocally determined by the local pinning force counterbalanced by the driving (Lorentz) force induced by the dc bias current. A positive dc bias makes the vortices to “probe” the gentle groove slope, whereas the steep groove slope is probed by vortices when they are subject to a negative dc bias. (g–i) The coordinate dependences of the pinning potentials U(x) deduced numerically using the approach outlined in Supplementary (solid lines) are plotted together with fits (dashed lines) to Eq. (2) and the cross-sectional AFM line scans across the grooves (symbols). The light red areas in panels (g–i) correspond to the left-hand WPP groove slopes in panels (a–c) probed at positive dc bias values (red solid symbols in panels (a–c)). The light blue areas in panels (g–i) correspond to the right-hand WPP groove slopes in panels (a–c) probed at negative dc bias values (blue open symbols in panels (a–c)).

[Back]

[Back]

Table 1.

Main parameters of the investigated Nb microstrips.

Parameter | Sample S | Sample A1 | Sample A2
d, nm | 40 | 56 | 70
T c, K | 8.66 | 8.85 | 8.94
f ₀, GHz | 5.72 | 3.95 | 3.12
jd+ at T = 0.3T c, MA/cm² | 0.75 | 0.7 | 0.52
jd− at T = 0.3T c, MA/cm² | 0.75 | 0.91 | 1.25
jd−/jd+, deduced from the CVC | 1 | 1.3 | 2.4
e, deduced asymmetry parameter | 0 | 0.13 | 0.56

[Back]