 interesting ways, which can be observed experimentally, using microscopic probes. For instance, both the order parameter, $\eta (\bf R)$, and magnetic field, $\bf B (R)$, are altered near the impurity. These perturbations therefore provide important insight into the nature of pinning. From a theoretical perspective, the calculation of the perturbations in $\eta (\bf R)$ and $\bf B (R)$ constitutes a non-trivial application of the Thuneberg theory.

In this work, we investigate the changes of the order parameter $\delta \eta (\bf R)$, and magnetic field $\delta \bf B(R)$, due to the impurity, when a vortex core is placed on a defect. Near the defect these quantities can be calculated analytically. Away from the defect, $\delta \eta (\bf R)$ decays exponentially over the correlation length scale, $\xi (T)$, while $\delta \bf B(R)$ decays exponentially over the penetration length scale $\lambda (T)$.

The paper is organized as follows. In Section II we review the appropriate GL theory, and set up the general equations for our problem. These equations, (\[eq:def1\]) and (\[eq:dea1\]), are the main results of this section. They are a pair of coupled, linear, partial differential equations for the impurity-induced changes in the order parameter and vector potential.

In Section III we consider in particular the spatial region close to the impurity, where $\delta \eta (\bf R)$ and $\delta \bf B(R)$ are largest. In this region we are able to derive explicit expressions for the perturbations—the simplicity of the result for $\delta \bf B(R)$ is somewhat remarkable. Eqs. (\[eq:dpsi2\]-\[eq:da2\]) contain our main results. In Section IV we discuss the results, and estimate the magnitude of the perturbation effects. In the appendix we calculate, numerically, the GL solutions for vortices in the absence of impurities.

GL Equations
============

General Formulation
-------------------

For simplicity, in this paper we consider a superconductor with a spin singlet, isotropic ([*i.e.*]{} $\hat{\bf k}$-independent) energy gap. The GL free energy, in terms of the complex order parameter $\eta 
(\bf R)$ and vector potential $\bf A (R)$, is given, in the absence of the impurity, by $$\begin{aligned}
\Omega_{\rm pure} & =  & \
interesting ways, which can be observed experimentally, using microscopic probes. For case, both the decree parameter, $ \eta (\bf R)$, and magnetic field, $ \bf B (R)$, are change near the impurity. These perturbations therefore supply important insight into the nature of trap. From a theoretical perspective, the calculation of the perturbations in $ \eta (\bf R)$ and $ \bf B (R)$ establish a non - trivial application of the Thuneberg hypothesis. 

 In this work, we investigate the changes of the order parameter $ \delta \eta (\bf R)$, and magnetic field $ \delta \bf B(R)$, due to the impurity, when a vortex kernel is placed on a defect. Near the defect these quantity can be calculated analytically. Away from the defect, $ \delta \eta (\bf R)$ disintegrate exponentially over the correlation length scale, $ \xi (T)$, while $ \delta \bf B(R)$ decays exponentially over the penetration length plate $ \lambda (T)$. 

 The paper is organized as follows. In Section   II we review the appropriate GL theory, and set up the general equation for our problem. These equations, (\[eq: def1\ ]) and (\[eq: dea1\ ]), are the main results of this section. They are a pair of coupled, linear, partial differential equations for the impurity - induce changes in the order argument and vector electric potential. 

 In Section III we consider in particular the spatial region close to the impurity, where $ \delta \eta (\bf R)$ and $ \delta \bf B(R)$ are big. In this region we are able to derive explicit expressions for the disruption — the simplicity of the result for $ \delta \bf B(R)$ is somewhat remarkable. Eqs.   (\[eq: dpsi2\]-\[eq: da2\ ]) hold our main results. In Section IV we discuss the results, and calculate the magnitude of the disruption consequence. In the appendix we calculate, numerically, the GL solutions for vortices in the absence of impurities. 

 GL Equations 
 = = = = = = = = = = = = 

 General Formulation 
 ------------------- 

 For simplicity, in this paper we consider a superconductor with a tailspin singlet, isotropic ([ * i.e. * ] { } $ \hat{\bf k}$-independent) energy col. The GL free energy, in term of the complex order parameter $ \eta 
 (\bf R)$ and vector potential $ \bf A (R)$, is given, in the absence of the impurity, by $ $ \begin{aligned } 
 \Omega_{\rm pure } & =   & \
 inheresting ways, which can be observed exkeeimentelly, usjng micruscopic probes. For instance, uoth the irder parameter, $\eta (\bf R)$, and mahnetic fueld, $\vf B (R)$, are altered near fme im'ucity. These pertorbations tharefore provida km'ortant insight into the nature of pynning. Grlm a theoreticwl ptrs[ectjne, the calculation of the pertudbationv in $\eta (\bf R)$ and $\bf B (R)$ constitutes a jon-tgivial application of the Thubebewt theory.

In tfis work, wt nnvestigate the changes of the order paramdter $\belta \eta (\bd E)$, ajg magnetic hield $\qelta \bf B(R)$, due to tve imputity, when a vovtex rore is placed on a defecv. Near the defect thgse quantidizs can be calculated qnqlytiwallf. Awxt ffom tie sefect, $\demta \eta (\bf R)$ decays ezponentially over tne bprrelation lsngth fcwle, $\xi (T)$, while $\delta \bf B(R)$ decays expontntiamly over the penetratiob length scale $\lambda (T)$.

The papqr is organized as follows. In Section II we review the eporokrlwgw HL theory, and set up the general equations fow oir problem. These equations, (\[rq:feg1\]) and (\[eq:dea1\]), arg the mcjn results of this sfction. Jhey aee a pair of voupled, linear, partial diffwrential equctiins for the impuricy-induced chcnges on thr order parameter and vzctor lotential.

In Section JKI we consider iv psrdicular ufe spatial region close to the impurigy, wnere $\dqlta \eta (\bv R)$ akg $\delta \bf B(R)$ are pargedt. In this rfgion we are able to derive explicit expressipnv fmr the pzrturbstions—the sim[licity of the result for $\aelta \bf B(G)$ is somexhat remarkafle. Eqs. (\[eq:dpsi2\]-\[ae:da2\]) contain mur main resylts. In Secgkon IV we disciss the rvsblts, and wstimate the magnibude ur the perturbatnub effects. In thr aopegdpx xe cajwulate, numerhcaluy, gne GL solutions nor voryices in the absence of jmpurities.

GL Equatooks
============

General Formulaeion
-------------------

For simplocity, in this papeg we ronsidxr a sipetconductor with a spin singlet, jsotropic ([*i.e.*]{} $\hat{\bf k}$-inde[endcnt) gnergy gap. Che GL free energy, in terms of the complxx order parameter $\eta 
(\vf R)$ and vector potgntlal $\bf A (R)$, ix givqn, in the absence of the impueity, by $$\begin{aligked}
\Omega_{\rm pure} & =  & \
interesting ways, which can be observed experimentally, probes. instance, both order parameter, $\eta $\bf (R)$, are altered the impurity. These therefore provide important insight into the of pinning. From a theoretical perspective, the calculation of the perturbations in $\eta R)$ and $\bf B (R)$ constitutes a non-trivial application of the Thuneberg theory. this we the of the order parameter $\delta \eta (\bf R)$, and magnetic field $\delta \bf B(R)$, due to impurity, when a vortex core is placed on defect. Near the defect quantities can be calculated analytically. from defect, $\delta (\bf decays over the correlation scale, $\xi (T)$, while $\delta \bf B(R)$ decays exponentially over the penetration length scale $\lambda (T)$. The is organized In Section we the GL theory, and the general equations for our problem. and (\[eq:dea1\]), are the main results of this They are pair of coupled, linear, partial differential for the impurity-induced changes in the order parameter vector potential. In Section III we consider in particular the spatial region close to the $\delta \eta (\bf R)$ $\delta \bf B(R)$ largest. this we able to explicit expressions for the perturbations—the simplicity of the result for $\delta B(R)$ is somewhat remarkable. Eqs. (\[eq:dpsi2\]-\[eq:da2\]) contain our main results. IV discuss the results, estimate the magnitude of perturbation In the appendix we the solutions the of GL Equations ============ General ------------------- For simplicity, in this we consider a superconductor ([*i.e.*]{} $\hat{\bf k}$-independent) energy gap. The GL free in terms of the complex order parameter (\bf R)$ and vector potential $\bf A (R)$, is given, in the of the $$\begin{aligned} \Omega_{\rm pure} & = & \
 interesting ways, which can be Observed exPerimEntAllY, uSing MicrOscopic probes. FOR insTance, both the order paramEter, $\eTa (\BF R)$, anD MaGnetiC field, $\bF b (R)$, ARE alTeReD neAr THe ImpurIty. these peRturbationS thErEfore provide IMpOrtant insiGht Into the naturE of PinninG. FRom A TheorEtiCal peRspectIVe, the cAlculatioN oF The perTUrbatioNS In $\Eta (\bF R)$ and $\bf B (R)$ constituTEs A Non-trivial applIcatioN oF ThE tHunEbeRg theory.

In ThIs worK, We invesTIgATE The CHanges of the orDer parameteR $\DelTa \eta (\bF R)$, And MAgnetiC fielD $\dELta \Bf B(R)$, due to thE impUrity, when A vorteX Core is pLAced on a Defect. neaR thE defECt ThEse QuANtiTIeS caN Be cAlculateD aNaLyticAlly. aWAY From The DefeCt, $\delTa \eta (\bf R)$ decayS exPoneNTiaLly ovEr the CorrElAtion Length Scale, $\Xi (t)$, while $\delta \bf B(R)$ DecaYs exponenTiaLlY ovEr The peNEtratiOn lEngTh scale $\Lambda (T)$.

tHe pApER IS oRganized as follows. IN SECTiOn II we reView thE ApPrOPriate GL ThEorY, and SET up thE genERaL equatioNs for oUR pRoBlem. TheSe EquatiOnS, (\[eq:Def1\]) And (\[eq:DEa1\]), arE the maIn resultS of thIS section. They arE A pair of coupleD, LiNEAr, PArtiAl dIfferential EquaTIons For tHE iMpuRIty-inDuced ChANgES in the order parameteR aNd vectOr potEntial.

In SectiOn III we conSIDEr in partIculAR tHE spatial region Close To the impurITy, where $\dElta \eTa (\bf R)$ and $\Delta \bf B(R)$ ARE largest. in tHis RegIon WE ArE able to derive EXPlicIt ExpressIonS for the PerTurBatIonS—tHe simplicIty of the ReSuLt FoR $\deLta \bf b(r)$ is somewHaT reMaRkaBle. EqS. (\[Eq:dpsi2\]-\[Eq:da2\]) cOntaIn OuR MaiN resultS. in sECtioN Iv wE disCusS tHe resUlts, ANd eStimate The magnitUde OF the PeRtUrbatioN effects. In the ApPendix we caLcUlaTe, numeRICally, the gL solutions for vortices iN The abseNce Of impUritIes.

GL EquaTioNs
============

GeneRal fOrmulaTion
-------------------

FoR simpLiCitY, IN this PAPeR we CoNsider a supERConDuctoR wIth a Spin sinGlet, isotropic ([*i.e.*]{} $\hat{\BF k}$-iNdependent) eneRgy Gap. THE gL FreE EnERgy, In TErmS OF the complex ordeR parameter $\EtA 
(\Bf r)$ and vector POteNtIal $\bf A (R)$, Is given, In the ABsence oF the impurIty, by $$\begiN{aLignED}
\omeGa_{\rm pure} & =  & \
 interesting ways, which c an be obse rvedexp eri me ntal ly,using microsco p ic p robes. For instance, b oth t he orde r  p arame ter, $\ e ta ( \bf  R )$ , a nd ma gneti c f ield, $ \bf B (R)$ , a re  altered nea r  t he impurit y.These pertur bat ions t he ref o re pr ovi de im portan t  insig ht into t he nature of pinn i n g.  Fro m a theoretical p e rs p ective, the ca lculat io n  o f the  pe rturbation sin $\ e ta (\bf R) $ a nd$ \bf B (R)$ co nstitutes a non -trivi al  ap p licati on of  t h e T huneberg th eory .

In thi s work ,  we inv e stigate  the c han ges  oft he  o rde rp ara m et er$ \de lta \eta  ( \b f R)$ , an d m a gnet icfiel d $\d elta \bf B(R) $,duet o t he im purit y, w he n a v ortexcoreis  placed on a de fect . Near th e d ef ect  t heseq uantit ies  ca n be ca lculate d  an al y t i ca lly. Away from the  d e f ec t, $\del ta \et a  ( \b f  R)$ dec ay s e xpon e n tiall y ov e rthe corr elatio n  l en gth sca le , $\xi  ( T)$ , w hile$ \del ta \bf  B(R)$ d ecays exponentiallyo ver the penet r at i o nl engt h s cale $\lamb da ( T )$.
The pa per is or ganiz ed as follows. In Section  I I we r eview  the appropri ate GL the o r y , and se t up th e  general equat ionsfor our pr o blem. Th ese e quations , (\[eq:d e f 1\]) and  (\ [eq :de a1\ ] ) ,are the mainr e sult sof this  se ction.The y a rea p ai r of coup led, lin ea r,  p ar tia l dif f erential  e qua ti ons  fort he imp urity -ind uc ed cha nges in th e orde rpa rame ter  a nd ve ctor pot ential. 

In Sect ion IIIwe  c onsider  in particula rthe spatia lreg ion cl o s e to the  impurity, where $\delt a  \eta ( \bf  R)$and$\delta \ bfB(R)$are larges t. Inthisre gio n we ar e ab leto  derive ex p l ici t exp re ssio ns forthe perturbations— t hesimplicity of  th e re s u lt  fo r  $ \ del ta \bf B (R)$ is somewha t remarkab le .  E qs. (\[eq: d psi 2\ ]-\[eq: da2\])conta i n our m ain resul ts. In Se ct ionI V  we  discuss t he resul ts, and e s timat e  t he ma gni tude o fthe  pert urbati o n e ffect s. Inth e appe ndixwe  calcula te, numerically, the GL  solut ionsfor  vortices  in the  absenceof i mpurities. 

G L E quati ons ===== ==== = == 

G e neral  For m ulation
- - -- --- - - -- --------

F o r sim plici ty, in thi s pa per we consider a superconductor  wit h a s pin sing le t, isotropic ( [*i .e . * ]{} $\ha t{ \bf k}$-ind ependent )e nergy  gap.The GL  free e n e rg y , in t erms  of  the comp lex  o r der par am et e r $\et a 
( \b f R)$and ve c torp o tential $\bf A ( R)$,i s  give n , i n the  a bsenceo f th e impurity , by $$\beg in{ali gned }
\Om ega_{\r mpure}& =  & \
 interesting_ways, which_can be observed experimentally,_using microscopic_probes._For instance,_both_the order parameter,_$\eta (\bf R)$,_and magnetic field, $\bf_B (R)$, are_altered_near the impurity. These perturbations therefore provide important insight into the nature of pinning._From_a theoretical_perspective,_the_calculation of the perturbations in_$\eta (\bf R)$ and $\bf_B (R)$_constitutes a non-trivial application of the Thuneberg theory.

In_this_work, we investigate_the changes of the order parameter $\delta \eta (\bf_R)$, and magnetic field $\delta \bf_B(R)$, due to_the_impurity,_when a vortex core_is placed on a defect. Near_the defect these quantities can be_calculated analytically. Away from the defect, $\delta_\eta (\bf R)$ decays exponentially over_the correlation length scale, $\xi_(T)$, while_$\delta \bf B(R)$ decays exponentially_over the penetration_length scale_$\lambda (T)$.

The paper_is organized as follows. In Section II_we review the_appropriate GL theory, and set up_the_general equations for_our_problem._These equations,_(\[eq:def1\]) and (\[eq:dea1\]),_are_the main_results_of this section. They are a_pair_of coupled, linear, partial differential equations for_the impurity-induced changes in_the_order parameter and vector_potential.

In Section III we consider_in particular the spatial region close_to the_impurity, where_$\delta \eta (\bf R)$ and $\delta \bf B(R)$ are largest. In_this region we are able to_derive explicit expressions for_the perturbations—the_simplicity_of the result_for_$\delta \bf_B(R)$ is somewhat remarkable. Eqs. (\[eq:dpsi2\]-\[eq:da2\]) contain our_main results._In Section IV we discuss the_results, and estimate the_magnitude_of the perturbation effects. In the_appendix we calculate, numerically, the GL_solutions for vortices in the_absence_of_impurities.

GL Equations
============

General Formulation
-------------------

For simplicity, in_this paper we consider a superconductor_with a spin_singlet, isotropic ([*i.e.*]{} $\hat{\bf k}$-independent) energy gap._The_GL free energy, in terms of_the_complex order parameter $\eta 
(\bf R)$_and_vector_potential $\bf A (R)$, is_given, in the absence of the_impurity, by $$\begin{aligned}
\Omega_{\rm pure} & =  & \