Vortices in superconductors with a columnar defect: Finite size effects

https://doi.org/10.1016/j.physc.2005.02.007Get rights and content

Abstract

In the present work we investigate the behavior of a vortex in a long superconducting cylinder near to a columnar defect at the center. The derivations of the local magnetic field distribution and the Gibbs free energy will be carried out for a cylinder and a cavity of arbitrary sizes. From the general expressions, it is considered two particular limits: one in which the radius of the cavity is very small but the radius of the superconducting cylinder is kept finite; and one in which the radius of the superconducting cylinder is taken very large (infinite) but the radius of the cavity is kept finite. In both cases the maximum number of vortices which are allowed in the cavity is determined. In addition, the surface barrier field for flux entrance into the cavity is calculated.

Introduction

The analysis of vortex motion in the presence of defects is one of the most important topic in the superconductivity research field. The relevance of this theme is based on the fact that the behavior of type-II superconductors in the presence of an applied magnetic field can be described by the motion and pinning of vortices in the material.

In the last decades a number of magnetic anomalies, such as the Wohlleben effect, fish tail anomaly and the jumps on the magnetic response of mesoscopic samples, were reported relating the occurrence of defects and the granular conformance of both bulk samples and thin films [1], [2]. Simultaneously, several initiatives using the most different theoretical approaches have been carried out to describe granular superconductors and their characteristic magnetic behavior [3], [4]. Today it is well known that all these features and anomalies are straight linked to the vortex dynamics.

Recently, the advances in the fabrication process related to nanotechnology made possible the production of mesoscopic rings [5] and superconducting arrays with columnar defects [6]. Experiments using these nanoscopic samples showed that the connectivity and vorticity are quite different than those for the macroscopic ones. The physics revealed by those experiments has contributed to a major understanding of type-II granular superconductors. This has motivated both experimental and theoretical physicists to investigate the magnetic properties of either small superconductors or bulk samples with nanosized defects. The properties of the vortex lattice, for instance, in this superconductors of confined geometries change radically with respect to the bulk superconductor ones.

One of the open questions within the framework of mesoscopic superconducting materials is the validity and complete understanding of the simple rule proposed by Mkrtchyan and Shmidt [7]. This rule indicates that the maximum number of vortices that enter a columnar defect is ns = r/2ξ, where r is the radius of the cylindrical cavity and ξ is the superconducting coherence length. Although this rule has been questioned, it is a reasonable estimate of the cavity occupation [8], [9].

In this paper the work of Refs. [7], [10] is extended for the case of a long superconducting wire, of arbitrary size, with a columnar defect. The applied magnetic field is parallel to the cylinder axis and the cylindrical cavity axis coincides with the direction of the applied field. First we determine the local magnetic field with appropriate boundary conditions and further, we calculate the Gibbs free energy of an ensemble of vortices near the cavity with n vortices inside it. On using the Gibbs free energy, the force acting on a vortex near the cavity and near to the external surface of the superconductor is calculated. From the expression of the force, we find a criterion for the maximum number of vortices allowed in the cavity. It is shown that if the size of the cavity is not too small, the saturation number differs substantially from the classical rule of Mkrtchyan and Shmidt [7]. The surface barrier field for flux penetration into the cavity is also determined.

Section snippets

The magnetic field

The starting point of our study is the London equation for the local magnetic field of a very strong type-II superconductor for which the Ginzburg–Landau parameter κ=λξ1. For our purposes, this equation can be more conveniently written in cylindrical coordinates as-λ22hρ2+1ρhρ+1ρ22hϕ2+h=Φ0i=1Nδ(ρ-ρi),where ρi is the position of the i-vortex line, Φ0 is the quantum flux, and N is the number of vortices outside the columnar defect.

Consider a long superconducting wire under an applied

The Gibbs free energy

In order to proceed we need to evaluate the Gibbs free energy (per unit length) which is given byG=F-ABH4π,where F is the London free energy (the Helmholtz free energy in the thermodynamic context) and B is the average induction; here A is the area of the cylinder cross section. In the London approximation, the free energy is given byF=18π02πraλ2ρ2hϕ2+λ2hρ2+h2ρdρdϕ.By using integration by parts, we can writeF=Φ08πi=1Nh(ρi,ϕi)+λ2aH8π02πhρρ=adϕ-λ2rH08π02πhρρ=rdϕ.Next, we compute the

Force

We will be most interested in the particular case N = 1. From Eq. (20), we haveG=Φ028πG(ρ1)+Φ0H4πΔ0(r,a)[Δ0(ρ1,a)+Δ0(r,ρ1)]+λrH0(ρ1)2δ0(r,a)Δ0(r,a)H0(ρ1)2-H-λaH24δ0(a,r)Δ0(r,a)+λ2H221Δ0(r,a)-Φ0H4π.Now the Green’s function is given byG(ρ1)=12πλ2K0(ξ/λ)+12πλ2m=-Km(a/λ)Im(ρ1/λ)Δm(ρ1,r)+Km(ρ1/λ)Im(r/λ)Δm(a,ρ1)1Δm(r,a),and the magnetic field in the cavity isH0(ρ1)=Φ02πλrn+Δ0(ρ1,a)Δ0(r,a)+2πλ2HΦ0Δ0(r,a)r2λ+δ0(r,a)Δ0(r,a).The force acting on a single vortex can be found by taking the derivative of

Summary and conclusions

In summary, we have derived the local field and the free energy of an ensemble of vortices around a columnar defect in a superconducting wire for both the cavity and the superconductor of arbitrary sizes. We also evaluated the force near to the cavity and external surfaces of the superconductor. It has also been found that for a large superconductor, the cavity saturates at a larger number of vortices as its size increases. However, not linearly in r/2ξ as predicted in [7]. In fact, in [8], [9]

References (12)

  • W.A.C. Passos et al.

    Physica C

    (2001)
  • Y. Fouad et al.

    Physica C

    (1994)
  • R. de Luca

    Physica C

    (2004)
  • D.A. Luzhbin et al.

    Phys. Rev. B

    (2004)
  • X. Zhang et al.

    Phys. Rev. B

    (1997)
    S. Pedersen et al.

    Phys. Rev. B

    (2001)
  • V.V. Moshchalkov et al.

    Phys. Rev. B

    (1996)
There are more references available in the full text version of this article.
View full text