Abstract We describe by use of molecular dynamics simulations how the attraction, at short range, between the vortex cores aids the clustering of vortices. Our simulations suggest that the short range core attraction and ensuing vortex clustering is the mechanism responsible for the second peak in the magnetisation loops of superconductors. 

D K Jackson

Henrik Jeldtoft Jensen

N A Lindop

CiteSeerX

Ž .Physica C 332 2000 58–60
www.elsevier.nlrlocaterphysc
Simulations of clustering of vortices in superconductors: a
possible origin of the second peak
D.K. Jackson, N.A. Lindop 1, Henrik Jeldtoft Jensen )
Department of Mathematics, Imperial College, London SW7 2BZ, UK
Abstract
We describe by use of molecular dynamics simulations how the attraction, at short range, between the vortex cores aids
the clustering of vortices. Our simulations suggest that the short range core attraction and ensuing vortex clustering is the
mechanism responsible for the second peak in the magnetisation loops of superconductors. q 2000 Elsevier Science B.V. All
rights reserved.
Keywords: Magnetization; Vortex clustering; Molecular dynamics simulation
1. Introduction
There is still no consensus as to the explanation of
the second peak observed in the magnetisation loop
of high temperature superconductors and of some
w xlow temperature superconductors 1 . Experimen-
tally, one observes that the peak structure in YBCO
becomes more pronounced with decreasing tempera-
ture. From this, we expect that a zero temperature
mechanism must be at the heart of the origin of the
w xpeak 2 . This is our motivation for performing zero
temperature molecular dynamics simulations.
We have elsewhere briefly presented some of our
simulations in a discussion which focusses on the
relationship to experiments and more general conse-
w xquences of the vortex clustering 2 . In the present
) Corresponding author.
Ž .E-mail addresses: david@ic.ac.uk D.K. Jackson ,
Ž .h.jensen@ic.ac.uk H.J. Jensen .
1 Present address: Department of Chemistry, University of Ed-
inburgh, Edinburgh EH9 3JJ, UK
paper, we present a more detailed description of our
molecular dynamics simulation.
We consider an infinite slab of superconductor in
a magnetic field applied parallel to the slab’s face.
We choose the London limit in so far as we consider
discrete flux lines. The energy of interaction between
w xtwo segments of flux lines is 4 ,
f 20
U r sŽ .vv 2X2p l mŽ . 0
=
X Xyr rl yrrje e
™ ™ ™ ™
< < < <d l Pd l y d l d l1 2 1 2r r
1Ž .
™ ™< <where rs r yr is the distance between the two1 2
™ ™
™ ™segments d l and d l located at positions r and r ,1 2 1 2
X X' (respectively. l slr 1yb and j sjr 2 1ybŽ .
are the effective field dependent London penetration
depth and coherence length, respectively. bsBrBc 2
is the reduced magnetic induction relative to the
upper critical induction B .c 2
0921-4534r00r$ - see front matter q 2000 Elsevier Science B.V. All rights reserved.
Ž .PII: S0921-4534 99 00644-9
( )D.K. Jackson et al.rPhysica C 332 2000 58–60 59
For numerical convenience we treat the vortices
as straight parallel flux lines so that we can move
Ž .from a three-dimensional 3-D to a two-dimensional
Ž .2-D model. The interaction energy per unit length
becomes,
f 20 X XU r s K rrl yK rrj , 2Ž . Ž . Ž . Ž .vv 0 0X 22pl
w xwhere K is a modified Bessel function. Brandt 40
has shown this potential to reproduce to a good
approximation the elastic properties of a vortex lat-
tice found by Ginzburg–Landau theory and the mi-
croscopic theory of Gorkov. Its long-range repulsive
component reduces to the London potential at low
field. The attractive term is due to a change in the
condensation energy with the overlapping of vortex
cores.
Ž .It is immediately clear by comparing Eqs. 1 and
Ž .2 that if clustering is found for the 2-D case
clustering will also occur in three dimensions. This is
because the repulsive term of the interaction energy
Ž .in Eq. 1 depends on the relative orientation of the
two vortex segments and has a maximum when the
segments are parallel — as they are in the 2-D
situation. The attractive core-part of the interaction
energy does not depend on the orientation and will
accordingly be relatively more important for non-
parallel vortex lines. Clearly, much future research is
needed in order to understand quantitatively the rela-
tionship between our model and experiment. For
Ž .instance, the potential used, Eq. 2 , is still only an
approximation, though it does incorporate the core
attraction.
We note that for a system of evenly distributed
vortices the separation of the vortices is well charac-
terised by the average vortex separation a0
˚ ' Žs f rB s450 Ar B with the magnetic induc-( 0
tion measured in tesla and f denotes the flux0
.quantum . Clearly, at experimentally accessible
fields, a 4j since the coherence length, j , is0
˚typically of, say, 10 to 20 A. For this reason, simula-
tions have hitherto neglected the second term in Eq.
Ž . w x2 3 . A homogeneous vortex distribution is ex-
pected to arise in field-cooled experiments.
The situation is very different when magnetisation
loops are measured. Here, the external field is ramped
up and down forcing the vortices to move in and out
of the sample. In doing so, the vortices have to push
each other through the pinning potential present in-
side the superconducting sample. The shortest sepa-
ration encountered between a vortex being pushed in
Ž .or out of the sample and a vortex trapped at a
position in the pinning potential will be determined
by the mechanical force balance between the force
Ž .derived from Eq. 2 and the pinning force. It is easy
to estimate that physically realistic pinning forces
can be larger that the maximum of the vortex–vortex
force. This is a purely two-body argument. In reality,
many-body effects will come into play and tend to
lower the repulsive vortex force by smoothing.
2. Simulations
Ž .We use a deterministic zero temperature over
damped equation of motion. The pinning potential is
constructed by randomly positioning Gaussian-
shaped pinning wells of a concentration n , ampli-p
tude A , and range j . The amplitude is a fraction ofp
the superconducting condensation energy contained
per length in a cylinder of radius j .
Here, we show results from molecular dynamics
simulations for a system with periodic boundary
conditions of size 70j by 70j . A slice of the system
is left unpinned, and vortices are introduced here to
represent ramping of the applied field, an effective
w xcomputational method first used by Ref. 3 . The rest
of the system, representing the superconductor con-
tains the randomly distributed Gaussian pins. The
interaction force between the vortices is cut off at
half the smaller dimension of the system. This cutoff
is made smoothly to avoid false vortex configuration
w xthat may result from a sharp cutoff 5 .
As a measure of the magnetisation M we use the
gradient of the vortex density in the pinned region,
i.e., d Brd x. In the case of a linear vortex profile
Žinside the sample a good approximation at high
.induction , one has d Brd xAM. We plot this gradi-
ent as function of the average vortex–vortex density
in the pinned region, i.e., B. In this way, we do not
have to define the external magnetic field.
3. Simulation results
The inclusion of the attractive term in the vortex
interaction makes a considerable difference to the
( )D.K. Jackson et al.rPhysica C 332 2000 58–6060
Fig. 1. Magnetisation averaged over four different pinning configurations. On the left is a system with a vortex–vortex interaction
Ž Ž ..containing only a repulsive term, on the right, the attractive term is also used see Eq. 2 . Solid line is the magnetisation, dashed line is the
˚fraction of vortices within j of another vortex. Pin amplitude 0.3 of condensation energy. ks70.7. Fields calculated for js15 A.
magnetisation loop, most noticeably at high fields
Ž .see Fig. 1 The attractive term allows vortices to
cluster; a significant fraction of vortices sit within j
of another vortex, as Fig. 1 shows. The magnetisa-
tion results vary between different pinning configura-
tions. However, our computational method limits us
to small systems so this ensemble of results is aver-
aged to obtain a result more typical of the macro-
scopic scale.
4. Conclusion
We conclude that the second peak may require us
to consider situations where vortices come close
together. In such a situation, the magnetisation can-
not be described within a framework, say, such as
elastic medium theory, which considers long wave-
Žlength deformations possibly with small density of
.topological defects and which assume a separation
between vortices of order the average separation
a s f rB . The core attraction, which is most(0 0
often assumed to be of negligible importance, allows
for a clustering of vortices at favourable positions in
the pinning potential. We expect this clustering to
have significant dynamical consequences in mag-
netic relaxation experiments.
Finally, we remark that the clustering described
here is equivalent to large fluctuations in the local
superconducting current density.
Acknowledgements
We are grateful to Mario Nicodemi and Garry
Perkins for collaboration and discussions. DKJ and
NL are grateful to EPSRC and the Imperial College
Centre for High Temperature Superconductivity for
studentships. HJJ is supported by EPSRC.
References
w x Ž .1 L.F. Cohen, H.J. Jensen, Rep. Prog. Phys. 60 1997 1581.
w x2 D.K. Jackson, M. Nicodemi, G.K. Perkins, N.A. Lindop, H.J.
Jensen, cond-matr9908454.
w x3 C. Reichardt, C.J. Olson, J. Groth, S. Field, F. Nori, Phys.
Ž .Rev. B 52 1995 10441.
w x Ž .4 E.H. Brandt, Phys. Rev. B 34 1986 6514.
w x5 H. Fangohr, A. Price, S. Cox, P.A.J. de Groot, G.J. Daniell
Efficient methods for handling long-range forces in particle–
particle simulations, in press.


Simulations of clustering of vortices in superconductors: a possible origin of the second peak

http://citeseerx.ist.psu.edu/viewdoc/download;jsessionid=B3BE32B0D6E1740E7258983B31373784?doi=10.1.1.1043.5456&rep=rep1&type=pdf

Simulations of clustering of vortices in superconductors: a possible origin of the second peak

Abstract

Similar works

Full text

Available Versions

CiteSeerX