Mass of Abrikosov vortex in high-temperature superconductor YBa\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2Cu\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}3O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{7-\delta }$$\end{document}7-δ

For more than four decades, mass of Abrikosov vortices defied experimental observations. We demonstrate a method of its detection in high-temperature superconductors. Similarly to electrons, fluxons circulate in the direction given by the magnetic field, causing circular dichroism. We report the magneto-transmittance of a nearly optimally doped thin YBa\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2Cu\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}3O\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{7-\delta }$$\end{document}7-δ film, measured using circularly polarized submillimeter waves. The circular dichroism emerges in the superconducting state and increases with dropping temperature. Our results confirm the dominant role of quasiparticle states in the vortex core and yield the diagonal fluxon mass of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2.2 \times 10^8$$\end{document}2.2×108 electron masses per centimeter at 45 K and zero-frequency limit, and even larger off-diagonal mass of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.9 \times 10^8 m_e$$\end{document}4.9×108me/cm.


Far-infrared measurements of the fluxon mass
Magneto-transmittance of circularly polarized light. We designed and developed a unique far-infrared (FIR) transmission experiment capable of probing the circular dichroism (see Fig. 1). The breakthrough that has eventually allowed us to conduct this research is a custom-made retarder 22 inserted in the optical path near the FIR-laser output aperture; it delays the horizontal polarization relative to the vertical one by an adjustable phase difference. A mutual phase shift of ±π/2 converts the THz beam of equal vertical and horizontal polarization components into the circularly polarized state. Since the phase delay introduced by the retarder is inversely proportional to the wavelength of the incoming light, each laser line requires a separate adjustment. We measured the transmittance of the sample for several laser lines at wavelengths 119, 163, 312, 419, and 433 µm 23 , corresponding to the terahertz circular frequencies ω of 15.9, 11.6, 6.05, 4.50, and 4.35 ×10 12 rad/s, respectively. Our setup shown in Fig. 1a enables fast flips between clockwise ( + ) and anti-clockwise (−) circular polarizations; consequently, we probed the transmittance T + and T − of both polarizations under identical conditions. The transmittance, i.e., the fraction of laser energy transmitted through the sample, is evaluated as the bolometer-to-pyrodetector signal ratio, effectively eliminating any possible time instability in the laser power. Identical profiles of both signals confirm that the transmittance is measured in the linear regime.
Our experimental protocol is as follows: In experimental runs, we apply a magnetic field at a temperature well above T c . As the temperature drops, vortices freeze into a regular Abrikosov lattice. The applied magnetic field B is kept constant since any change in B would result in inhomogeneous patterns of the vortex density. The temperature T is swept down and up at a steady sweep rate of 2.5 K/min, the instantaneous temperature of the sample being recorded together with the transmittance. Care is taken to avoid hysteresis of the down and up sweeps. The continuous laser beam is split by a linear wire-grid polarizer; the reflected part is monitored using a pyroelectric detector to keep trace of unavoidable power fluctuations, whereas the transmitted part proceeds toward the sample. The retarder converts the light from linear to circular polarization. The propagation of the circularly polarized beam and the magnetic field are perpendicular to the film surface, as detailed in panel (b). (b) Vortices in the film (gray circles) control the transmittance of the sample via the following mechanism: The electric field of the laser light drives the supercurrent. The Magnus force accelerates vortices in the direction perpendicular to the supercurrent; in reaction, the vortex motion affects the supercurrent and, thus, the transmittance. In the sketch, the electric field in the sample, as well as the vortices, rotate clockwise. If the light frequency is close to the cyclotron frequency of vortices, the motion of vortices is resonantly enhanced, leading to the observed dichroism. The extent of the cyclotron motion is strongly exaggerated; in fact, the fluxon circulates on a radius of less than 10 −12 m at the strongest laser line. (c) Transmittance of the YBa 2 Cu 3 O 7−δ superconducting sample, normalized to the normal-state transmittance T N at 100 K and plotted for two circular polarizations versus temperature. The dichroism is clearly visible below 70 K in a magnetic field of 10 T; in a zero field, no dichroism appears. (d) Transmittance ratio T + /T − measured in several applied magnetic fields plotted versus temperature. The data were obtained using a 312 µ m laser line ( 6.1 × 10 12 rad/s). Above the critical temperature, the dichroism is absent, showing that the normal-state Hall component is negligible. www.nature.com/scientificreports/ Figure 1c displays a typical temperature behavior of the transmittance observed in the external magnetic field of 10 T using a 312 µ m laser line. The results obtained with different laser lines are similar. As expected, the transmittance measured in zero field does not show any dichroism, in absence of Abrikosov vortices to induce an asymmetry. In nonzero fields, however, the low-temperature circular dichroism is clearly observed and can be attributed to the formation of vortices threading the sample. The effect is enhanced in higher applied magnetic fields, thanks to the growing areal density of Abrikosov vortices (see Fig. 1d). We focus on the low-temperature region, T < 50 K, where the thermal quasiparticles play a negligible role and the response of the system is fully governed by the motion of fluxons.
Sample specification by time-domain terahertz spectroscopy. We chose the most common high-T c material YBa 2 Cu 3 O 7−δ in the form of a thin film with a thickness of L = 107 nm, and with CuO 2 planes parallel to the surface. The sample was prepared at National Chiao Tung University (Taiwan) using a pulsed laser deposition method from a stoichiometric target on a lanthanum aluminate substrate oriented in the (100) plane. The substrate dimensions are 10 × 10 × 0.5 mm 3 . Several measurements were performed to establish sample properties. Figure 2 summarizes some relevant results.
The critical temperature of the film, T c = 87.6 K, was determined from a dedicated measurement of dc resistivity ρ (Fig. 2a). For our slightly underdoped sample, this T c corresponds to the hole density 24 n 0 = 1.68 × 10 27 /m 3 . The linear slope of ρ extrapolated to zero temperature shows negligible residual resistivity, so the relaxation time is inversely proportional to temperature.
Additional film properties were established in a separate experiment using standard time-domain THz spectroscopy. Broadband linearly polarized THz pulses (0.3-2.5 THz) were generated by exciting an interdigitated LT-GaAs emitter with a Ti:sapphire femtosecond laser beam at 800 nm . 25 We measured the complex conductivity σ = σ ′ + iσ ′′ for frequencies ω/2π in the range 0.5-2 THz at temperatures from 4 to 100 K (Fig. 2b,d,e). At the zero magnetic field and below T c , the two-fluid model fit confirms the dominant London contribution, σ 0 ≈ ine 2 /(mω) , with a temperature-dependent condensate density n = n 0 (1 − T 4 /T 4 c ) , the hole mass m = 3.3 m e , and the elementary charge e. Comparing the conductivities at temperatures of 4 K and T N = 100 K in Fig. 2d, we found the relaxation time τ N = 5 × 10 −14 s. Conductivities at all temperatures from 4 to 100 K are consistent with τ = τ N T N /T.
Our sample is moderately clean. Its purity is given by the lifetime measured on the energy scale 19 as is the Boltzmann constant and is the reduced Planck constant. The Fermi energy E F = 2 k 2 F /2m depends on the hole doping. For the hole density of our sample, the Fermi surface is cylindrical rather than spherical, and the Fermi momentum follows from the 2D density of holes in the CuO 2 plane n 2D = n 0 c/2 = k 2 F /(2π) , where c = 11.68 Å is the YBaCuO lattice parameter in the z-direction. The resulting Fermi energy E F = 71 meV yields the value of k 2 B T 2 c τ/ E F ∼ 0.1 , corresponding to the moderately clean sample. According to Kopnin and Vinokur 19 , the moderately clean d-wave superconductor behaves as the s-wave one. In the absence of reliable formulas for angular frequency ω 0 of quasiparticles bounded in the vortex core in the d-wave superconductors, we used To complete the sample characteristics, we assume the pinning of vortices by layer imperfections, for example, the surface roughness. Figure 2e shows the conductivity in the magnetic field of 7 T applied perpendicularly to the film 28,29 . It was interpreted either with the theory specified below or with the model used by Parks 30 , both indicating that the vortex pinning is rather weak with the Labusch coefficient κ ≈ 2 × 10 5 N/m 2 .
Theoretical prediction. Our aim is to compare experimental values of T + /T − with a theoretical model.
Using the above sample parameters, we can evaluate the film conductivity σ ± from which the transmittance T ± results. The theoretical prediction shown in Fig. 3 is based on the Yeh formalism 31 and covers interferences in the weakly birefringent substrate. For the purpose of discussion, we refer to an approximation T + /T − = |σ − | 2 /|σ + | 2 , which differs from the exact theory by less than 4% as shown in the Supplementary Information . 32 Kopnin and Vinokur 33 provide the theoretical conductivities σ ± derived under very general conditions. Our study allows for two simplifications. First, we focus on temperatures about 45 K, where we observe a large dichroic signal. At such low temperatures, extended quasiparticles are very dilute so we can neglect their contribution to the electric current J as well as their effect on vortices. Second, the cyclotron frequency ω c = eB/m is small on the scale of the quasiparticle lifetime, ω c τ ≪ 1 , for all experimental magnetic fields, which simplifies dynamics of quasiparticles in the vortex core. Under these conditions, the equation of motion for a fluxon of unitary length takes the form of Newton's law 33 with the time derivative of momentum p and the force F from the interaction of the vortex core with the crystal lattice. Unlike Kopnin and Vinokur, we include the pinning force with the Labusch parameter κ and vortex displacement u related to its velocity as u = v . The Magnus force, given by the vector product of the magnetic field direction z = B/B and the vortex velocity related to the condensate current, covers a force by which the flowing condensate acts on the fluxon.
In Kopnin's model, the fluxon momentum p is a total momentum of quasiparticles in the vortex core which rotate about the vortex axis at an angular frequency ω 0 . For our moderately clean sample, the fluxon momentum at the low-temperature limit depends on its velocity as 19  www.nature.com/scientificreports/ The diagonal mass µ is complex at finite frequencies, which reflects a delay between a change of the vortex velocity and a change of the total momentum of quasiparticles in its core. The off-diagonal mass µ ⊥ describes a property common in anisotropic systems that the velocity of an excitation is not parallel with its momentum. Quasiparticles disturbed by the vortex motion and the action of FIR light eventually lose their momentum in collisions with impurities and phonons. Via these collisions, the crystal lattice acts on the fluxon by force F . With the collision integral approximated by a single relaxation time τ , the force and the momentum are simply related by F = −p/τ . The longitudinal force −(µ � /τ )v is a vortex friction. The transversal force (µ ⊥ /τ )[v × z] reduces the Magnus force giving the Kopnin-Kravtsov force 20,21 in the dc limit.
The electric field in the film 33

Discussion
Experiment versus theoretical prediction. With the full set of sample parameters established from the time-domain spectroscopy, the simplified Kopnin-Vinokur conductivity (4) furnishes us with the theoretical prediction of the circular dichroism. Figure 3 compares this prediction and the experimentally observed dichroism for several values of the applied magnetic field and the THz-laser wavelength. The differences between theory and experiment are smaller than the experimental errors.
With increasing wavelength and field strength, the transmittance ratio T + /T − gradually deviates from unity. The observed trends can be understood in terms of Eq. (3). The field dependence arises from the dominant London contribution 1/σ 0 complemented by the Josephson-type resistivity, which is linear in B . The variation with wavelength has a similar cause: for lower frequencies, the London resistivity 1/σ 0 ≈ −iωm/(ne 2 ) is smaller so the Josephson part becomes dominant.
Vortex mass from experiment. Figure 3 documents that the theory of Kopnin and Vinokur is relevant for the THz dynamics of vortices. Since the observed frequency dependence of magneto-transmission agrees with the theoretical one, even with no adjustable parameters, the extrapolation of our results to low frequencies is justified. Based on this, we used their theory to find the vortex mass from our data.
While the sample parameters n, σ 0 , and τ are sound, κ and ω 0 are less clear. The Labusch parameter κ has a minor effect on the dichroism; therefore, we kept the value κ = 2 × 10 5 N/m 2 . The angular frequency ω 0 was established by the least-square fit of T + /T − data. The best-fit value ω 0 = 4.3 × 10 12 rad/s was very close to 4.4 × 10 12 rad/s estimated above. We believe that such close agreement of observed and estimated frequency is fortuitous.
At THz frequencies, where the circular dichroism was found, both diagonal and off-diagonal masses are complex. They become real in the low-frequency limit, as apparent from Eq. (2). Using the experimentally established values of ω 0 = 4.3 × 10 12 rad/s and other sample parameters, we evaluated the zero-frequency components of the fluxon mass. In YBa 2 Cu 3 O 7−δ at 45 K, the diagonal mass µ amounts to 2.2 × 10 8 m e /cm, while the off-diagonal mass µ ⊥ is more than twice larger, 4.9 × 10 8 m e /cm.
In summary, we have developed a reliable experimental method to measure the circular dichroism of superconducting films threaded by Abrikosov vortices. To interpret our data in terms of vortex dynamics, we have established all the essential material parameters from independent time-domain THz spectroscopy and dcresistivity measurements. The observed dichroism is in good agreement with the theory of Kopnin and Vinokur based on the circular motion of quasiparticles in the vortex core. Their angular frequency was experimentally determined and used to extrapolate the vortex mass from THz frequencies to low-frequency motion.