Abstract: Depending on the Ginzburg-Landau parameter kappa, superconductors can either be fully diamagnetic if kappa < 1/root 2 (type I superconductors) or allow magnetic flux to penetrate through Abrikosov vortices if kappa > 1/root 2 (type II superconductors; refs 1,2). At the Bogomolny critical point, kappa = kappa(c) = 1/root 2, a state that is infinitely degenerate with respect to vortex spatial configurations arises(3,4). Despite in-depth investigations of conventional type I and type II superconductors, a thorough understanding of the magnetic behaviour in the near-Bogomolny critical regime at kappa similar to kappa(c) remains lacking. Here we report that in confined systems the critical regime expands over a finite interval of kappa forming a critical superconducting state. We show that in this state, in a sample with dimensions comparable to the vortex core size, vortices merge into a multi-quanta droplet, which undergoes Rayleigh instability(5) on increasing kappa and decays by emitting single vortices. Superconducting vortices realize Nielsen-Olesen singular solutions of the Abelian Higgs model, which is pervasive in phenomena ranging from quantum electrodynamics to cosmology(6-9). Our study of the transient dynamics of Abrikosov-Nielsen-Olesen vortices in systems with boundaries promises access to non-trivial effects in quantum field theory by means of bench-top laboratory experiments.

Abstract: Considering a noncentrosymmetric pinning texture composed of a square array of triangular holes, the magnetic flux penetration and expulsion are investigated experimentally and theoretically. A direct visualization of the magnetic landscape obtained using a magneto-optical technique on a Nb film is complemented by a multiscale numerical modeling. This combined approach allows the magnetic flux dynamics to be identified from the single flux quantum limit up to the macroscopic electromagnetic response. Within the theoretical framework provided by time-dependent Ginzburg-Landau simulations, an estimation of the in-plane current anisotropy is obtained and its dependence with the radius of the curvature of hole vertices is addressed. These simulations show that current crowding plays an important role in channeling the flux motion, favoring hole-to-hole flux hopping rather than promoting interstitial flux displacement in between the holes. The resulting anisotropy of the critical current density gives rise to a distinct pattern of discontinuity lines for increasing and decreasing applied magnetic fields, in sharp contrast to the invariable patterns reported for centrosymmetric pinning potentials. This observation is partially accounted for by the rectification effect, as demonstrated by finite-element modeling. At low temperatures, where magnetic field penetration is dominated by thermomagnetic instabilities, highly directional magnetic flux avalanches with a fingerlike shape are observed to propagate along the easy axis of the pinning potential. This morphology is reproduced by numerical simulations. Our findings demonstrate that anisotropic pinning landscapes and, in particular, ratchet potentials produce subtle modifications to the critical state field profile that are reflected in the distribution of discontinuity lines.

Abstract: Bogoliubov-de Gennes theory is used to investigate the effect of the size of a superconducting square on the vortex states in the quantum confinement regime. When the superconducting coherence length is comparable to the Fermi wavelength, the shape resonances of the superconducting order parameter have strong influence on the vortex configuration. Several unconventional vortex states, including asymmetric ones, giant-multivortex combinations, and states comprising giant antivortices, were found as ground states and their stability was found to be very sensitive on the value of k(F)xi(0), the size of the sample W, and the magnetic flux Phi. By increasing the temperature and/or enlarging the size of the sample, quantum confinement is suppressed and the conventional mesoscopic vortex states as predicted by the Ginzburg-Laudau (GL) theory are recovered. However, contrary to the GL results we found that the states containing symmetry-induced vortex-antivortex pairs are stable over the whole temperature range. It turns out that the inhomogeneous order parameter induced by quantum confinement favors vortex-antivortex molecules, as well as giant vortices with a rich structure in the vortex core-unattainable in the GL domain.

Abstract: Based on the Gor'kov formalism for a clean s-wave superconductor, we develop an extended version of the single-band Ginzburg-Landau (GL) theory by means of a systematic expansion in the deviation from the critical temperature T(c), i.e., tau = 1 – T/T(c). We calculate different contributions to the order parameter and the magnetic field: the leading contributions (proportional to tau(1/2) in the order parameter and. t in the magnetic field) are controlled by the standard GL theory, while the next-to-leading terms (proportional to tau(3/2) in the gap and proportional to tau(2) in the magnetic field) constitute the extended GL (EGL) approach. We derive the free-energy functional for the extended formalism and the corresponding expression for the current density. To illustrate the usefulness of our formalism, we calculate, in a semianalytical form, the temperature-dependent correction to the GL parameter at which the surface energy becomes zero, and analytically, the temperature dependence of the thermodynamic critical field. We demonstrate that the EGL formalism is not just a mathematical extension to the theory: variations of both the gap and the thermodynamic critical field with temperature calculated within the EGL theory are found in very good agreement with the full BCS results down to low temperatures, which dramatically improves the applicability of the formalism compared to its standard predecessor.

Abstract: An experimental determination of the mean vortex velocity in superconductors mostly relies on the measurement of flux-flow resistance with magnetic field, temperature, or driving current. In the present work we introduce a method combining conventional transport measurements and a frequency-tuned flashing pinning potential to obtain reliable estimates of the vortex velocity. The proposed device is characterized using the time-dependent Ginzburg-Landau formalism, where the velocimetry method exploits the resonances in mean vortex dissipation when temporal commensuration occurs between the vortex crossings and the flashing potential. We discuss the sensitivity of the proposed technique on applied current, temperature and heat diffusion, as well as the vortex core deformations during fast motion.

Abstract: After the discovery of magnetism in monolayer CrI3, the magnetic properties of different 2D materials from the chromium-trihalide family are intuitively assumed to be similar, yielding magnetic anisotropy from the spin-orbit coupling on halide ligands. Here we reveal significant differences between the CrI3 and CrBr3 magnetic monolayers in their magnetic anisotropy, resulting Curie temperature, hysteresis in external magnetic field, and evolution of magnetism with strain, all predominantly attributed to distinctly different interplay of atomic contributions to spin-orbit coupling in two materials.

Abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

Abstract: A thorough understanding of the skyrmion motion through nanotracks is a prerequisite to realize the full potential of spintronic applications like the skyrmion racetrack memory. One of the challenges is to place the data, i.e. skyrmions, on discrete fixed positions, e.g. below a read or write head. In the domain-wall racetrack memory, one proposed solution to this problem was patterning the nanotrack with notches. Following this approach, this paper reports on the skyrmion mobility through a nanotrack with periodic notches (constrictions) made using variations in the chiral Dzyaloshinskii-Moriya interaction. We observe that such notches induce a coupling between the mobility and the skyrmion breathing mode, which manifests itself as velocity-dependent oscillations of the skyrmion diameter and plateaus in which the velocity is independent of the driving force. Despite the fact that domain walls are far more rigid objects than skyrmions, we were able to perform an analogous study and, surprisingly, found even larger plateaus of constant velocity. For both systems it is straightforward to tune the velocity at these plateaus by changing the design of the notched nanotrack geometry, e.g. by varying the distance between the notches. Therefore, the notch-induced coupling between the excited modes and the mobility could offer a strategy to stabilize the velocity against unwanted perturbations in racetrack-like applications. In the last part of the paper we focus on the low-current mobility regimes, whose very rich dynamics at nonzero temperatures are very similar to the operating principle of recently developed probabilistic logic devices. This proves that the mobility of nanomagnetic structures through a periodically modulated track is not only interesting from a fundamental point of view, but has a future in many spintronic applications.

Abstract: A superconducting diode is an electronic device that conducts supercurrent and exhibits zero resistance primarily for one direction of applied current. Such a dissipationless diode is a desirable unit for constructing electronic circuits with ultralow power consumption. However, realizing a superconducting diode is fundamentally and technologically challenging, as it usually requires a material structure without a centre of inversion, which is scarce among superconducting materials. Here, we demonstrate a superconducting diode achieved in a conventional superconducting film patterned with a conformal array of nanoscale holes, which breaks the spatial inversion symmetry. We showcase the superconducting diode effect through switchable and reversible rectification signals, which can be three orders of magnitude larger than that from a flux-quantum diode. The introduction of conformal potential landscapes for creating a superconducting diode is thereby proven as a convenient, tunable, yet vastly advantageous tool for superconducting electronics. This could be readily applicable to any superconducting materials, including cuprates and iron-based superconductors that have higher transition temperatures and are desirable in device applications. A superconducting diode is dissipationless and desirable for electronic circuits with ultralow power consumption, yet it remains challenging to realize it. Here, the authors achieve a superconducting diode in a conventional superconducting film patterned with a conformal array of nanoscale holes.

Abstract: A Raman scattering study covering the entire substitution range of the FeSe1-xSx solid solution is presented. Data were taken as a function of sulfur concentration x for 0 <= x <= 1, of temperature and of scattering symmetry. All types of excitations including phonons, spins, and charges are analyzed in detail. It is observed that the energy and width of the iron-related B-1g phonon mode vary continuously across the entire range of sulfur substitution. The A(1g) chalcogenide mode disappears above x = 0.23 and reappears at a much higher energy for x = 0.69. In a similar way the spectral features appearing at finite doping in A(1g) symmetry vary discontinuously. The magnetic excitation centered at approximately 500 cm(-1) disappears above x = 0.23 where the A(1g) lattice excitations exhibit a discontinuous change in energy. The low-energy mode associated with fluctuations displays maximal intensity at the nematostructural transition and thus tracks the phase boundary.

Abstract: A numerical approach to GinzburgLandau (GL) theory is demonstrated and we review its applications to several examples of current interest in the research on superconductivity. This analysis also shows the applicability of the two-dimensional approach to thin superconductors and the re-defined effective GL parameter κ. For two-gap superconductors, the conveniently written GL equations directly show that the magnetic behavior of the sample depends not just on the GL parameter of two bands, but also on the ratio of respective coherence lengths.

Abstract: A Josephson phase shift can be induced in a Josephson junction by a strategically nearby pinned Abrikosov vortex (AV). For an asymmetric distribution of an imprinted phase along the junction (controlled by the position of the AV) such a simple system is capable of rectification of ac current in a broad and tunable frequency range. The resulting rectified voltage is a consequence of the directed motion of a Josephson antivortex which forms a pair with the AV when at local equilibrium. The proposed realization of the ratchet potential by an imprinted phase is more efficient than the asymmetric geometry of the junction itself, is easily realizable experimentally, and provides rectification even in the absence of an applied magnetic field.

Abstract: <script type='text/javascript'>document.write(unpmarked('We present an advanced method to study spin fluctuations in superconductors quantitatively and entirely from first principles. This method can be generally applied to materials where electron-phonon coupling and spin fluctuations coexist. We employ it here to examine the recently synthesized superconductor iron tetraboride (FeB4) with experimental T-c similar to 2.4 K [H. Gou et al., Phys. Rev. Lett, 111, 157002 (2013)]. We prove that FeB4 is particularly prone to ferromagnetic spin fluctuations due to the presence of iron, resulting in a large Stoner interaction strength, I = 1.5 eV, as calculated from first principles. The other important factor is its Fermi surface that consists of three separate sheets, among which two are nested ellipsoids. The resulting susceptibility has a ferromagnetic peak around q = 0, from which we calculated the repulsive interaction between Cooper pair electrons using the random phase approximation. Subsequently, we combined the electron-phonon interaction calculated from first principles with the spin fluctuation interaction in fully anisotropic Eliashberg theory calculations. We show that the resulting superconducting gap spectrum is conventional, yet very strongly depleted due to coupling to the spin fluctuations. The critical temperature decreases from T-c = 41 K, if they are not taken into account, to T-c = 1.7 K, in good agreement with the experimental value.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('Spin-triplet chiral p-wave superconductivity is typically described by a two-component order parameter, and as such is prone to unique emergent effects when compared to the standard single-component superconductors. Here we present the equilibrium phase diagram for small mesoscopic chiral p-wave superconducting disks in the presence of magnetic field, obtained by solving the microscopic Bogoliubov-de Gennes equations self-consistently. In the ultrasmall limit, the cylindrically symmetric giant-vortex states form the ground state of the system. However, with increasing sample size, the cylindrical symmetry is broken as the two components of the order parameter segregate into domains, and the number of fragmented domain walls between them characterizes the resulting states. Such domain walls are topological defects unique for the p-wave order, and constitute a dominant phase in the mesoscopic regime. Moreover, we find two possible types of domain walls, identified by their chirality-dependent interaction with the edge states.'));