Next Article in Journal
Soliton Solutions of Mathematical Physics Models Using the Exponential Function Technique
Next Article in Special Issue
Vortex Dynamics of Charge Carriers in the Quasi-Relativistic Graphene Model: High-Energy k · p Approximation
Previous Article in Journal
Modulational Instability, Inter-Component Asymmetry, and Formation of Quantum Droplets in One-Dimensional Binary Bose Gases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relation of Superconducting Pairing Symmetry and Non-Magnetic Impurity Effects in Vortex States

1
Department of Physics, Okayama University, Okayama 700-8530, Japan
2
Research Institute for Interdisciplinary Science, Okayama University, Okayama 700-8530, Japan
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(1), 175; https://doi.org/10.3390/sym12010175
Submission received: 30 December 2019 / Revised: 14 January 2020 / Accepted: 15 January 2020 / Published: 18 January 2020
(This article belongs to the Special Issue Vortex, Topology and Singularity in Quantum Systems)

Abstract

:
Non-magnetic impurity scattering effects on the vortex core states are theoretically studied to clarify the contributions from the sign-change of the pairing function in anisotropic superconductors. The vortex states are calculated by the Eilenberger theory in superconductors with p x -wave pairing symmetry, as well as the corresponding anisotropic s-wave symmetry. From the spatial structure of the pair potential and the local electronic states around a vortex, we examine the differences between anisotropic superconductors with and without sign-change of the pairing function, and estimate how twofold symmetric vortex core images change with increasing the impurity scattering rate both in the Born and the unitary limits. We found that twofold symmetric vortex core image of zero-energy local density of states changes the orientation of the twofold symmetry with increasing the scattering rate when the sign change occurs in the pairing function. Without the sign change, the vortex core shape reduces to circular one with approaching dirty cases. These results of the impurity effects are valuable for identifying the pairing symmetry by observation of the vortex core image by the STM observation.

1. Introduction

Superconductivity is caused by the long-range order of the pair potential, which is wave function of Cooper pair. In type II superconductors [1,2,3], magnetic fields penetrate into the superconductor as quantized flux lines, and form a flux line lattice under applied fields between upper and lower critical fields. Around a flux line with a vortex of supercurrent, the complex value pair potential has the phase winding 2 π around the vortex center, and the amplitude of the pair potential vanishes at the vortex center. At the vortex core where the pair potential is suppressed, low-energy bound states called Caroli-de Gennes Matricon states appear [4,5,6]. The local electronic structure of the bound state at the vortex core was experimentally observed by the scanning tunneling microscope (STM) observation by Hess et al. [7,8,9]. The STM observations of the vortex core image have been performed to study exotic properties of superconductivity in many superconductors [10,11,12,13], including high- T c cuprate superconductors [14,15,16,17,18,19], iron-based superconductors [20,21,22,23,24], and topological superconductors [24,25].
In unconventional superconductors with pairing mechanism other than the conventional electron–phonon interaction, anisotropic pairing symmetries such as p-, d-, and f-wave symmetries are realized rather than the isotropic s-wave. Therefore, to understand the pairing mechanism of unconventional superconductivity, we have to identify the pairing symmetry including the node positions of the pairing function. For the purpose, we can obtain important feature of the anisotropic superconductivity from the local density of states (LDOS) around a vortex in type-II superconductors. For example, in isotropic s-wave superconductors, the LDOS shows circular vortex core image, whereas in d-wave superconductors with line nodes, we see a fourfold symmetric vortex core image [26,27,28,29]. Reflecting the anisotropic superconducting gap, the star-shaped vortex core images are observed by STM experiments in NbSe 2 [8,30,31] and YNi 2 B 2 C [10,11,32]. In the direction where the LDOS extends far from the vortex center, pairing function of the superconductivity has nodes or gap-minimum on the Fermi surface. Recently, in quasi-one-dimensional nickel-compounds such as NiBi 3 , STM experiments observed twofold symmetric vortex core image [13], which may reflect the Fermi surface structure and anisotropic superconducting gap. Twofold symmetric vortex core image was observed also in a nematic superconductor FeSe [20,21,22].
An example for pairing function φ ( θ ) of anisotropic superconductivity with vertical line nodes is φ p x ( θ ) = ( 2 k x / k F = ( 2 cos θ for the p x -wave symmetry pairing. This is schematically shown in Figure 1a. Here, for simplicity we consider a two-dimensional cylindrical Fermi surface k = ( k x , k y , k z ) = ( k F cos θ , k F sin θ , k z ) . In the p x -wave symmetry pairing, nodes of the pairing function appear at θ = ± π / 2 . As the zero-energy LDOS N ( E = 0 , r ) around a vortex extends toward the node directions, vortex core image becomes twofold symmetric shape, which is expected to be observed by STM. However, in the clean limit the LDOS structure is reproduced even if we consider anisotropic s-wave pairing functions ϕ | p x | ( θ ) = | ϕ p x ( θ ) | in Figure 1b. As zero-energy LDOS is determined by the amplitude | φ ( θ ) | and independent of the sign of φ ( θ ) in the clean limit, it is difficult to distinguish the line node of φ p x ( θ ) and the deep-minimum of φ | p x | ( θ ) . On the other hand, it is know that non-magnetic impurity scattering effect is sensitive to the pairing symmetry. In uniform states, although isotropic s-wave pairing superconductivity is not affected by the impurity effect as called Anderson theorem [33,34], anisotropic s-wave superconductivity is weakly suppressed by the impurity scattering effect [35,36]. As the impurity effect easily suppresses unconventional p and d wave superconductivities [37,38], the strength to the impurity scattering is used as a clue to distinguish the pairing symmetry of superconductors. Therefore, to find a clue to identify the sign change of the pairing function φ ( θ ) , we have to study the impurity effect also in the vortex states in anisotropic superconductors. As other type of impurity effects, local bound states may appear at the impurity site [39,40], which are not covered in this work.
The purpose of this work is to theoretically study the non-magnetic impurity scattering effects on the spatial structure of vortex states. Especially, we focus the contribution of the sign change of the pairing function on the Fermi surface, comparing the results of p x -wave paring symmetry φ p x ( θ ) with sign change, and the anisotropic s-wave φ | p x | ( θ ) without sign change. We evaluate the spatial structure of the pair potential and the LDOS in the vortex state by quantitative calculation of selfconsistent Eilenberger theory in the vortex lattice state [28,29,41,42,43,44] and study the dependence on the scattering rate 1 / τ 0 of non-magnetic impurities both in the Born and the unitary limits. Around the twofold symmetric vortex core, the anisotropy ratio is easily evaluated comparing the core radius of the vortex core image along the x and the y directions. From the results of the vortex structure, we discuss contributions from the sign change of φ ( θ ) . The study about the anisotropic s-wave pairing can be applied to the twofold symmetric vortex core such as in NbNi 3 [13] and FeSe [20,21,22], which have deep minimum in the pairing function.
This paper is organized as follows. After this introduction, we explain our formulation in Section 2. In Section 3, we see the impurity effects in the uniform states. In Section 4 we study impurity effects on the pair potential and the LDOS in the vortex states for the cases of φ p x ( θ ) and φ | p x | ( θ ) . The last section is devoted to discussion and conclusion.

2. Formulation

2.1. Eilenberger Theory with Non-Magnetic Impurity Scatterings

Our numerical calculations of the vortex state are based on the Eilenberger theory [28,29,41,42,43,44]. Eilenberger equation is derived from Gor’kov equation in the superconducting state under the assumption Δ 0 / E F 1 , which is satisfied in many superconductors, for the superconducting gap Δ 0 and Fermi energy E F . The Eilenberger theory is suitable to treat inhomogeneous superconducting states, such as surface, interface, and vortex states. There, pair potential of superconductivity shows spatial variation in the length scale of the coherence length, which is longer than the atomic scale. To include contributions of non-magnetic impurities, we consider self-energies of the impurity scattering, as explained below. In addition to the Born limit [41,42], we cover also the unitary limit of the impurity scatterings by the t-matrix approximation [45,46,47,48,49,50].
Eilenberger equation is given by
i v · g ^ + i ω ˜ n τ ^ 3 Δ ^ Σ ^ , g ^ = 0 ^
in Nambu matrix form [45,50]. ω ˜ n = ω n + i v · A with Matsubara frequency ω n , and the vector potential A . v = v F / v F with Fermi velocity v F and v F 2 = v F 2 k . k indicates the Fermi surface average. The quasi-classical Green’s functions g ^ and the pair potential Δ ^ are, respectively,
g ^ ( ω n , k , r ) = i π g i f i f g , Δ ^ ( k , r ) = 0 Δ ( r ) φ ( k ) Δ ( r ) φ * ( k ) 0 .
In our formulation using the Eilenberger unit, length, temperature, and magnetic field are, respectively, measured in unit of ξ 0 , T c , and B 0 . Here, ξ 0 = v F / 2 π k B T c and B 0 = ϕ 0 / 2 π ξ 0 2 with the flux quantum ϕ 0 . T c is superconducting transition temperature in the clean limit at a zero magnetic field. The energy E, pair potential Δ and ω n are in unit of π k B T c .
In the t-matrix approximation for the non-magnetic impurity scattering [45,46,47,48,49,50], self-energy Σ ^ ( ω n , k , r ) = n s t ^ ( ω n , k , k , r ) is given by
t ^ ( ω n , k , k , r ) = u ^ ( ω n , k , k ) + N 0 u ^ ( ω n , k , k ) g ^ ( ω n , k , r ) t ^ ( ω n , k , k , r ) k .
N 0 is the density of states (DOS) at the Fermi energy in the normal state, and n s the number density of impurities. When impurity scattering is assumed to be s-wave scattering, impurity potential is a constant as u ^ = u 0 1 ^ , and t ^ ( ω n , r ) = u ^ + N 0 u ^ g ^ ( ω n , k , r ) k t ^ ( ω n , r ) . Therefore, as t ^ ( ω n , r ) = [ 1 N 0 u ^ g ^ ( ω n , k , r ) k ] 1 u ^ , we obtain the self-energy Σ ^ as
Σ ^ = n s 1 ^ + i π N 0 u 0 g k i f k i f k g k 1 u 0 = 1 τ 1 π N 0 u 0 1 ^ + i g k f k f k i g k
where
1 τ = 1 / τ 0 cos 2 δ 0 + ( g k 2 + f k f k ) sin 2 δ 0
with
1 τ 0 = n s π N 0 u 0 2 1 + ( π N 0 u 0 ) 2 , cos 2 δ 0 = 1 1 + ( π N 0 u 0 ) 2 , sin 2 δ 0 = ( π N 0 u 0 ) 2 1 + ( π N 0 u 0 ) 2 ,
and δ 0 = tan 1 ( π N 0 u 0 ) . τ 0 is in unit of 2 π k B T c / . In the Born limit of weak impurity scattering potential u 0 , δ 0 0 . In the unitary limit of strong scattering potential, δ 0 π / 2 .
From Equation (4) and [ 1 ^ , g ^ ] = 0 , Eilenberger Equation (1) is written as
i v · g ^ + i ( ω ˜ n + G ) ( Δ φ + F ) Δ φ * + F ( i ω ˜ n + G ) , g ^ = 0 ^
with components of the self-energy,
F ( ω n , r ) = 1 τ f k , F ( ω n , r ) = 1 τ f k , G ( ω n , r ) = 1 τ g k .
From matrix components of Equation (7), we obtain Eilenberger equations
ω n + G + v · + i A f = Δ φ + F g , ω n + G v · i A f = Δ * φ * + F g .
For numerical analyses of Eilenberger equation, it is more convenient to employ functions a ( ω n , k , r ) and b ( ω n , k , r ) defined by
f = 2 a 1 + a b , f = 2 b 1 + a b , g = 1 a b 1 + a b .
Substituting these into Eilenberger Equation (9), we obtain the Riccati differential equations [26,51],
v · a = Δ φ + F Δ * φ * + F a 2 ω n + G + i v · A a ,
v · b = Δ * φ * + F Δ φ + F b 2 ω n + G + i v · A b .
Unlike the original Equation (9), the Riccati equations for a and b are decoupled. The Riccati Equations (11) and (12) are then, respectively, solved by numerical integration along the trajectory parallel to the vector v , under given Δ , A , G, F, and F .

2.2. Methods of Numerical Calculations in the Vortex Lattice

We calculate the spatial structure of vortices in the vortex lattice state by selfconsistent Eilenberger theory [28,29,41,42,43,44,51,52]. We determine the spatial structure of vortices to be selfconsistent with the quasi-classical Green’s functions of electronic states, solving the Riccati differential equations. We introduce contributions of the impurity scatterings by self-energies from non-magnetic s-wave impurity scatterings [43,44,45,46,47,48,49,50]. This method of Eilenberger theory can be applied to all temperature and magnetic field range in the vortex state, and appropriately capture contributions of vortex core and inter-vortex interaction. In many previous researches, local electronic states in the vortex states were studied in the clean limit. The researches of the impurity effect in the vortex states were done in pairing symmetries with isotropic gap, such as isotropic s-wave and chiral p-wave. Although there were some studies of impurity effect in the vortex state in the d x 2 y 2 -wave pairing symmetry [44,47], we need systematic study of the impurity effect about the impurity scattering rate dependence and of the identification of contributions from the sign-change in the pairing function of anisotropic superconductors with line nodes.
In our calculation, magnetic fields B are applied z axis direction, so that B = ( 0 , 0 , B ) . The vector potential is given by A ( r ) = 1 2 B × r + a ( r ) in the symmetric gauge. For simplicity, we consider a two-dimensional cylindrical Fermi surface with a Fermi velocity v F = ( v x , v y , 0 ) ( cos θ , sin θ , 0 ) at a Fermi wave number k = ( k x , k y , k z ) = ( k F cos θ , k F sin θ , k z ) . We consider two cases of the pairing function, φ ( θ ) = φ p x ( θ ) and φ | p x | ( θ ) .
As the first step in our selfconsistent calculations of the vortex states, we calculate the quasiclassical Green’s functions in the Eilenberger theory, under given Δ , A , G, F, and F . To obtained the solution a ( ω n , k , r ) , we numerically integrate Equation (11) along the trajectory from r v s 0 to r [26,51,52]. As demonstrated in Figure 8 of [52], choosing enough long length s 0 of these trajectories we obtain unique solution independent of s 0 and initial values at r v s 0 . The solution b ( ω n , k , r ) is obtained in the similar way by integrating Equation (12) along the trajectory from r + v s 0 to r . From the solutions a and b, we obtain quasiclassical Green’s functions g, f, and f in Equation (10) at r .
As an initial state, we start our calculations from the Ablikosov lattice solution,
Δ ( r ) = 2 a y a x 1 4 p = exp π a y a x y + y 0 a y + p 2 2 π i p x 0 a x + ζ 2 p + y 0 a y + p × exp i π x y a x a y ,
and a ( r ) = 0 of the vector potential A ( r ) = 1 2 B × r + a ( r ) in the symmetric gauge. Unit vectors of the vortex lattice are given by r 1 = ( a x , 0 ) and r 1 = ( a x , ζ a y ) with ζ = 1 2 and B a x a y = ϕ 0 . We set a x / a y = 3 / 2 in the triangular vortex lattice. r 0 = ( x 0 , y 0 ) = 1 2 ( r 1 + r 2 ) , so that one of vortex centers is located at 1 2 ( r 1 + r 2 ) + r 0 = 0 . A unit cell in our calculation is set to be u 1 ( r 1 r 2 ) + u 2 r 2 with 0.5 u i < 0.5 ( i = 1 , 2 ).
Using the quasiclassical Green’s functions g, f, and f obtained from Equations (10)–(12), we estimate new values of Δ ( r ) and A ( r ) by the gap equation
Δ ( r ) = g 0 N 0 T 0 < ω n ω cut φ * f + f * k
and the current equation
× × A = 2 T κ 2 0 < ω n v Im g k .
We set
( g 0 N 0 ) 1 = ln T + 2 T 0 < ω n ω cut ω n 1
in Equation (14), so that the value of the cut-off ω cut does not seriously affect the numerical results. In the present calculation, ω cut = 20 k B T c . We set the Ginzburg–Landau (GL) parameter κ = 30 assuming typical type-II superconductors. For the large κ , since variation of internal fields due to a ( r ) is very small compared with the applied field as is shown in Ref. [44], contributions of small a ( r ) do not seriously affect on results of the pair potential and the LDOS around a vortex. The self-energies G, F and F are obtained from Equations (5) and (8) using the quasiclassical Green’s functions. We calculate Δ ( r ) , a ( r ) , and self-energies at mesh points r in a unit cell of the vortex lattice, and estimate them at arbitrary points by interpolation of the values on the mesh points. Once we obtain Δ ( r ) and a ( r ) in a unit cell, we can know Δ ( r ) and a ( r ) in other unit cells by the lattice translation R = m r 1 + n r 2 (m and n are integers) as
Δ ( r + R ) = Δ ( r ) e i χ ( r , R ) , a ( r + R ) = a ( r ) ,
where
χ ( r , R ) = π ϕ 0 ( B × R ) · ( r + 2 r 0 ) π m n
in the symmetric gauge.
First, we solve Equations (9)–(15) at T = 0.5 T c , and obtain selfconsistent solutions of Δ ( r ) , A ( r ) , and quasi-classical Green’s functions. We perform the calculations for some values of 1 / τ 0 in the Born and the unitary limits, in addition to the clean limit 1 / τ 0 = 0 . Using the renewed pair potential, vector potential, and self-energies, we solve the Riccati Equations (11) and (12) again to obtain new quasi-classical Green’s functions. We iterate these calculation steps, until a sufficiently self-consistent solution is obtained.
Next, we study local electronic states. Using the selfconsistently obtained Δ ( r ) and A ( r ) in the calculation of ω n , we solve Riccati Equations (11) and (12) with i ω n E + i η to obtain g ( i ω n E + i η , k , r ) , f ( i ω n E + i η , k , r ) , and f ( i ω n E + i η , k , r ) as a function of real energy E. η is an infinitesimal constant. Then, we calculate self-energies in Equation (8). These calculations of Riccati equations and the self-energies are iterated, until a selfconsistent solution is obtained. Using the selfconsistently obtained g, the LDOS is given by N ( E , r ) = N ( E , k , r ) k with k -resolved LDOS
N ( E , k , r ) = N 0 g ( i ω n E + i η , k , r ) .

3. Impurity Effects in Uniform States

Before studying the vortex states, we evaluate the strength of the superconductivity to the non-magnetic impurity scattering in uniform states without vortices at a zero magnetic field. Figure 2a shows the 1 / τ 0 -dependence of uniform Δ at T / T c = 0.5 to see how the impurity scattering suppresses p x -wave and anisotropic s-wave superconductivities. For the p x -wave pairing symmetry, Δ is slightly smaller in the unitary limit compared with the Born limit. In both limits Δ vanishes at 1 / τ 0 0.19 . On the other hand, by the Anderson theorem, Δ is not changed by the non-magnetic impurity scattering in isotropic s-wave symmetry. Even in the anisotropic s-wave pairing symmetry φ | p x | ( θ ) , Δ shows very weak suppression as is seen in Figure 2a. Therefore, s-wave superconductors are strong to the non-magnetic impurity scattering. Figure 2b shows Δ in wider range until higher 1 / τ 0 . There we see that Δ for φ | p x | ( θ ) survives until 1 / τ 0 ~22. The drop near 1 / τ 0 ~22 is related to ω cut . We see that Δ is smaller for larger ω cut when 1 / τ 0 is large. Our study is performed at 1 / τ 0 1 , where the ω cut -dependence of Δ is negligible. The Born and the unitary limits have almost the same 1 / τ 0 dependence for φ | p x | ( θ ) . The difference between two limits comes from L g k 2 + f k f k in Equation (5). For example, in isotropic s-wave pairing, since L = 1 from the relation g 2 + f f = 1 of the quasiclassical Green’s functions in Equation (10), 1 / τ = 1 / τ 0 in both limits.
Impurity effect on the DOS spectrum N ( E ) in uniform states for the p x wave pairing φ p x ( θ ) is presented in Figure 3a for the Born limit and in Figure 3b for the unitary limit. N ( E ) for the corresponding anisotropic s-wave φ | p x | ( θ ) is shown in Figure 3c for the Born limit, and in Figure 3d for the unitary limit. The horizontal axis is E / Δ with Δ for each 1 / τ 0 . Both φ p x ( θ ) and φ | p x | ( θ ) have the same N ( E ) in the clean limit (not shown in the figure), as the DOS is determined only by the amplitude of the pairing function, | φ ( θ ) | . We see a peak of N ( E ) at the gap edge E / Δ = 2 , as the maximum of | φ p x ( θ ) | is 2 . The peak is smeared with increasing the impurity scattering rate 1 / τ 0 both in the Born and unitary limits. The smearing is seen in all cases of Figure 3. In the presence of line nodes, it is know that N ( E ) E near E = 0 in the clean limit 1 / τ 0 = 0 . For the p x -wave symmetry in Figure 3a in the Born limit, N ( E ) = 0 even when 1 / τ 0 becomes finite. With increasing 1 / τ 0 , N ( E ) at low E ( 0 ) is enhanced. In the unitary limit in Figure 3b, finite DOS appears at E = 0 [53,54]. Flat N ( E ) near E = 0 becomes larger with increasing 1 / τ 0 . On the other hand, for anisotropic s-wave symmetry in Figure 3c,d, the Born and the unitary limits show the similar behaviors. Even in the unitary limit of the impurity scattering, the zero-energy states do not appear, as N ( E = 0 ) = 0 . Rather, at finite 1 / τ 0 gap opens, since N ( E ) = 0 at finite E near E = 0 . With increasing 1 / τ 0 , gap amplitude becomes larger.
The reason for the finite gap is because the gap structure is determined by Δ φ + F as is seen in Eilenberger Equations (9). In uniform states, the quasiclassical anomalous Green’s function f has the same symmetry as the pairing function φ ( θ ) . By the Fermi surface average, the anomalous self-energy F = 0 for φ ( θ ) = φ p x ( θ ) . However, for φ ( θ ) = φ | p x | ( θ ) , F 0 and finite gap are induced due to the finite F, which has isotropic s-wave symmetry on the Fermi surface. However, as the anomalous f and F show different behaviors at the vortex core from those in the uniform states, we need to study the impurity effect around the vortex core by detailed calculations of the spatial structure in the vortex state.

4. Impurity Effects in Vortex States for p x -Wave Pairing φ p x ( θ ) and Anisotropic s-Wave φ | p x | ( θ )

Figure 4a presents calculated spatial structure of the pair potential’s amplitude | Δ ( r ) | within a unit cell of the vortex lattice in x y plane for the p x -wave pairing symmetry φ p x ( θ ) at a low field B = 0.01 B 0 . The amplitude | Δ ( r ) | decreases with approaching the vortex center. The vortex core shape shows twofold symmetric shape, reflecting the twofold symmetric pairing function φ p x ( θ ) on the Fermi surface. This vortex core shape indicates the relation ξ x > ξ y , where ξ x ( ξ y ) is the coherence length relating to the vortex core radius in x (y) direction. This anisotropy of the coherence length is consistent to the estimate ξ x v x 2 | φ p x | 2 k 1 / 2 and ξ y v y 2 | φ p x | 2 k 1 / 2 in the GL theory. | φ p x ( θ ) | is larger (smaller) where v x 2 ( v y 2 ) is larger on the Fermi surface. Even when the impurity effect is stronger, vortex core shape keeps similar twofold symmetric anisotropy both in the Born and the unitary limits. The suppression by the impurity scattering is stronger in the unitary limit compared with the Born limit also around the vortex core. These behaviors are also seen from the comparison of the profile between x and y directions shown in the lower panels of Figure 4a. In each panel, | Δ ( r ) | and the slope around the vortex core are smaller along the x axis, compared to those along the y axis. When we compare the Born (solid lines) and the unitary (dashed lines) limits, the slope of | Δ ( r ) | at the vortex center is almost the same in the both limits. Approaching outside region of the vortex core, deviation of | Δ ( r ) | between two limits becomes larger. There, | Δ ( r ) | in the Born limit is larger compared with the unitary limit.
Figure 4b presents | Δ ( r ) | for anisotropic s-wave pairing φ | p x | ( θ ) . When the impurity scattering is weak with small 1 / τ 0 , the vortex core shape shows similar twofold symmetric shape as in Figure 4a. However, approaching the dirty case with larger 1 / τ 0 , the vortex core becomes circular shape. This is because contributions of the self-energy F with isotropic s-wave symmetry becomes dominant in the effective pair potential Δ φ + F . Differences between the Born and the unitary limits are small in this anisotropic s-wave pairing φ | p x | ( θ ) . These behaviors are also seen from the comparison of the profile between x and y directions shown in the lower panels of Figure 4b. In each panel, as lines for the Born limit overlap with those for the unitary limit, the x and y dependences of Δ ( r ) are same in both limits. Although | Δ ( r ) | is smaller along the x-direction, the deviation between x and y directions becomes small with increasing 1 / τ 0 to be circular vortex core shape.
From the fitting of the slope of Δ ( r ) at the vortex center, we define the core radius ξ x along the x-direction as Δ ( x , y = 0 ) = Δ max x / ξ x near x = 0 . Δ max is maximum of | Δ ( r ) | within a unit cell of the vortex lattice. Similarly we define the core radius ξ y along the y direction as Δ ( x = 0 , y ) = Δ max y / ξ y near y = 0 . We present 1 / τ 0 dependence of ξ x and ξ y in Figure 5a. Although each ξ increases as a function of 1 / τ 0 , and ξ x is longer than ξ y . For p x -wave pairing symmetry, ξ x and ξ y rapidly increase toward the critical point 1 / τ 0 0.19 . For anisotropic s-wave pairing symmetry, ξ x and ξ y gradually increase in this range of 1 / τ 0 , as the critical point 1 / τ 0 22 is far. The anisotropy ratio ξ y / ξ x is presented in Figure 5b as a function of 1 / ξ 0 . For the p x -wave, since the ratio ξ y / ξ x 0.65 in the range 0 1 / τ 0 < 0.19 of the superconductivity, twofold symmetric vortex core shape is unchanged by the impurity scattering, as seen Figure 4a. On the other hand, for anisotropic s-wave φ | p x | ( θ ) , the ratio ξ y / ξ x gradually increases towards 1 changing to circular vortex core shape seen in Figure 4b.
The zero-energy LDOS N ( E = 0 , r ) around a vortex core is presented in Figure 6a for the p x -wave pairing symmetry. N ( E = 0 , r ) has a peak at the vortex center, and extends toward the y-axis, which is the node direction of the pairing function φ p x ( θ ) on the Fermi surface. In the vortex lattice, N ( E = 0 , r ) is slightly suppressed on the line connecting nearest neighbor vortices [28,29]. This suppression is seen in N ( E = 0 , r ) on the y axis at 1 / τ 0 = 0.01 .
Figure 6a shows how the local electronic states around a vortex are changed by the non-magnetic impurity scattering effects in the Born and the unitary limits, respectively. The sharp structure of N ( E = 0 , r ) around the vortex core is smeared with increasing 1 / τ 0 , and the smearing effect is greater in the unitary limit than in the Born limit. With increasing the impurity effect, orientation of the twofold symmetric vortex core shape changes. At larger 1 / τ 0 , N ( E = 0 , r ) extends toward the x direction. This is similar anisotropy to that of | Δ ( r ) | in Figure 4a. In the unitary limit, we see enhancement of N ( E = 0 , r ) outside of vortex core by the impurity effect, which occurs even in uniform states as shown in Figure 3b. These behaviors are also seen in lower panels of Figure 4b. There, we see that N ( E = 0 , r ) at larger 1 / τ 0 is larger on the x-axis, compared with the y-axis direction, because twofold symmetric vortex core image extends toward the x-direction. In the unitary limit, peak height at the vortex center is lower, and N ( E = 0 , r ) is higher outside of the vortex core, compared with the Born limit.
On the other hand, impurity effect on the zero-energy LDOS N ( E = 0 , r ) is presented in Figure 6b for anisotropic s-wave pairing φ | p x | ( θ ) . When the impurity effect is weak with small 1 / τ 0 , N ( E = 0 , r ) shows similar anisotropic core image to that in Figure 6a, as the gap anisotropy is same as | φ p x ( θ ) | = | φ | p x | ( θ ) | . With increasing 1 / τ 0 of the impurity effect, peak structure of N ( E = 0 , r ) around the vortex core is smeared, and the vortex core image of N ( E = 0 , r ) transforms to circular shape. In this process, the extension of N ( E = 0 , r ) toward the node direction along the y axis is weakened, because contributions of the isotropic self-energy F become dominant. This changes to circular core shape is consistent to that of the pair potential in Figure 4b. The impurity effect in the unitary limit is qualitatively same as in the Born limit for the anisotropic s-wave pairing symmetry. These behaviors are also seen in lower panels of Figure 6b. Although peak height at the vortex center is higher in the unitary limit, the difference between the Born and the unitary limits becomes smaller with increasing 1 / τ 0 . Thus, N ( E = 0 , r ) has almost the same spatial structure in the two limits at 1 / τ 0 = 1 .
To see the 1 / τ 0 dependence quantitatively, we present the peak height N ( E = 0 , r = 0 ) at the vortex center as a function of 1 / τ 0 in Figure 7a. For the p x wave pairing φ p x ( θ ) in the range 0 1 / τ 0 0.15 , with increasing 1 / τ 0 , the peak height is suppressed toward normal state value 1. It is slightly smaller in the unitary limit than in the Born limit. On the other hand, for the anisotropic s-wave pairing φ | p x | ( θ ) in the range 0 1 / τ 0 1 , the suppression effect of the peak height by the impurity scattering is greater in the Born limit than in the unitary limit. Although the peak height is rapidly suppressed as in the p x -wave pairing in the Born limit, it decreases slowly in the unitary limit. These behavior in the anisotropic s-wave is consistent to results of isotropic s-wave pairing [47].
To quantitatively estimate the twofold symmetric vortex core shape in the zero-energy LDOS, the vortex core radii ξ x and ξ y are, respectively, determined as a radius of half width by N ( E = 0 , x = ξ x , y ) = 0.5 N ( E = 0 , x = y = 0 ) on the x axis and N ( E = 0 , x = 0 , y = ξ y ) = 0.5 N ( E = 0 , x = y = 0 ) on the y axis. The 1 / τ 0 dependences of ξ x and ξ y are presented in Figure 7b. There we see clear differences between the Born and the unitary limits, as the peak height is different between two limits as shown in Figure 7a. For p x -wave pairing symmetry, ξ x and ξ y rapidly increase toward the critical point 1 / τ 0 0.19 . The increase is greater in the unitary limit because N ( E = 0 , r ) largely appear outside of the vortex. For anisotropic s-wave pairing symmetry, ξ x and ξ y gradually increase in this range of 1 / τ 0 . The increase is smaller in the unitary limit than in the Born limit.
The anisotropy ratio ξ y / ξ x from N ( E = 0 , r ) is presented in Figure 7c as a function of 1 / ξ 0 . In the clean limit, ξ y / ξ x 1.7 , which is different anisotropy from the order parameter. For the p x -wave φ p x ( θ ) , the ratio ξ y / ξ x rapidly decreases as a function of 1 / τ 0 to ξ y / ξ x 0.5 , which is the same orientation of twofold symmetry as in Δ ( r ) . On the other hand, for anisotropic s-wave φ | p x | ( θ ) , the ratio decrease until ξ y / ξ x 1.0 to be circular vortex core image. The decrease of ξ y / ξ x is slower in the unitary limit, compared with the Born limit.
Last, we study the spatial structure of the LDOS N ( E , r ) at finite energy E = 0.01 , 0.1, 0.2, and 0.5. It is presented in Figure 8 for the p x -wave pairing φ p x ( θ ) , and in Figure 9 for the anisotropic s-wave pairing φ | p x | ( θ ) . These figures focus the vortex core region. At low energy E = 0.01 , spatial structure of the LDOS is almost the same as the zero-energy LDOS N ( E = 0 , r ) at 1 / τ 0 = 0.01 and 0.1 presented in Figure 6. At 1 / τ 0 = 0.05 in the Born limit, N ( E , r ) shows similar spatial structure both for the p x -wave φ p x ( θ ) in Figure 8 and the anisotropic s-wave φ | p x | ( θ ) in Figure 9. There, twofold symmetric vortex core image of the LDOS at E = 0.01 extends along the y axis, which is node-direction of the pairing function φ p x ( θ ) . At E = 0.1 and higher E, the vortex core shape is seen to be rotated 90 , and it extends to the x-direction. At E = 0.1 and 0.2, N ( E , r ) has large intensity near the vortex core, but it is suppressed at the vortex center. At higher energy E = 0.5 , N ( E , r ) is rather smaller in the vortex core region, and has peak on the y axis far from the vortex center. These behaviors are also seen in the upper panels of Figure 8b and Figure 9b. With increasing E, the peak position is shifted from the vortex center r = 0 at E = 0.01 to finite r. The distance between the peak position and the vortex center is longer at higher E. At 1 / τ 0 = 0.01 , N ( E , r ) in the unitary limit (dashed lines) has almost the same spatial structure as in the Born limit (solid lines) in Figure 8 and Figure 9.
The spatial pattern of the LDOS and the E-dependence at 1 / τ 0 = 0.01 are qualitatively explained by picture of the quasiparticle trajectory [27,31], which is schematically presented in Figure 10 for the p x -wave pairing symmetry φ p x ( θ ) . This is the same also for φ | p x | ( θ ) . As schematically presented in Figure 10a, the k -resolved LDOS N ( E , k , r ) for k with the velocity v = ( cos θ , sin θ , 0 ) is mainly distributed on the line with distance r ( E , θ ) from the vortex center. Roughly speaking, r is given by the relation E = | Δ ( r ) φ p x ( θ ) | . If we assume | Δ ( r ) | = Δ max r / ξ for simplicity, we obtain r = ξ ( E / Δ max ) / | φ p x ( θ ) | . When a quasiparticle is propagating toward an antinode direction θ = 0 or π , r takes the minimum value ξ E / Δ max . When propagating toward a node direction θ = ± π / 2 , r . On a line with angle α from the x-axis as shown in Figure 10a, distance s between the main distribution of N ( E , k , r ) and the vortex center is given by
s ( E , α , θ ) = r | sin ( θ α ) | = ξ E / Δ max | φ p x ( θ ) sin ( θ α ) | .
The LDOS N ( E , r ) is given by the sum of N ( E , k , r ) distribution in the range 0 θ < 2 π . The peak of the LDOS appears at the inner edge of the LDOS distribution, where distance s on a line of angle α takes the minimum s min ( E , α ) as a function of θ . For the p-wave pairing symmetry φ p x ( θ ) , s min = 2 ξ ( E / Δ max ) / ( 1 sin α ) at θ = ( α + π / 2 ) / 2 . In Figure 10b,c, s min on a line of angle α is plotted for 0 α < 2 π , as quasiparticle trajectories. On the trajectories, first, in-coming quasiparticles is pointing to the direction θ = π / 2 or π / 2 , where φ p x ( θ ) = 0 , far from the vortex center. On going around the vortex core, when the direction θ of the quasiparticle is increasing, the distance between the trajectories and the vortex center is narrower since | φ p x ( θ ) | becomes larger. When the direction of the trajectory is θ = 0 or π , the distance is the minimum at y = ± ξ E / Δ max . Then, trajectories go away from the vortex center towards the direction of θ = π / 2 or 3 π / 2 , with increasing θ . For small E in Figure 10b, tails of the trajectories extend toward ± y directions. This pattern of the trajectories corresponds to twofold symmetric vortex core image extending to the y-direction in the spatial structure of N ( E , r ) at E = 0.01 in Figure 8. On the other hand, for larger E in Figure 10c, shape of trajectories within the crossing points of two trajectories gives twofold symmetric vortex core image extending to the x-direction in N ( E , r ) at larger E in Figure 8, as tails of the trajectories far from the vortex are smeared.
When we compare the LDOS N ( E , r ) in the clean case, 1 / τ 0 = 0.01 , and the dirty case, 1 / τ 0 = 0.1 , for the p x wave symmetry φ p x ( θ ) in Figure 8, the spatial pattern of the quasiparticle trajectories is smeared in the dirty case. There, the trajectories are not seen far from the vortex. Around the vortex core, the distance between the trajectories and the vortex center increases, because E / Δ becomes larger due to the suppression of Δ by 1 / τ 0 , seen in Figure 2a. From lower panels of Figure 8, we see that the peak height becomes lower and the LDOS outside of vortex core is increased when 1 / τ 0 = 0.1 , compared with the case of 1 / τ 0 = 0.01 . These tendencies are greater in the unitary limit than in the Born limit. On the other hand, for anisotropic s-wave φ | p x | ( θ ) in Figure 9, we see weak suppression of the quasiparticle trajectory in the case of 1 / τ 0 = 0.1 , since anisotropic s-wave pairing is stronger to the impurity effect, compared with the p x wave pairing. From lower panels of Figure 9, we see that the peak height at the quasiparticle trajectories becomes lower at 1 / τ 0 = 0.1 , compared with that of 1 / τ 0 = 0.01 . The large differences between the Born and the unitary limits are only seen near the vortex core.

5. Discussion and Conclusions

To clarify contributions from the sign-change of the pairing function φ ( θ ) in the non-magnetic impurity scattering effects on the vortex state, we studied the twofold symmetric vortex core structure of the pair potential | Δ ( r ) | and the LDOS N ( E , r ) in anisotropic superconductors with p x -wave pairing symmetry φ p x ( θ ) and anisotropic s-wave pairing φ | p x | ( θ ) . The LDOS is accessible by the STM observation. In the STM experiments, twofold symmetric vortex core was observed in NiBi 3 [13] and FeSe [20,21,22], which may be anisotropic s-wave superconductors. If the impurity effects are experimentally studied in these superconductors, the theoretical results in this work may be examined by the STM observation. If a chiral p-wave superconductor exists, the p x -wave pairing superconductivity may be realized under uniaxial pressure, as tried in Sr 2 RuO 4 [55]. We note that twofold symmetric vortex core shape may also come from the contributions of anisotropic Fermi surface, in addition to anisotropic pairing symmetry. The study about the contribution of anisotropic Fermi surface belongs to future study.
As anisotropic superconductors with higher symmetry of the pairing, d x 2 y 2 -wave pairing is realized in high T c cuprate superconductors [56,57], and some of heavy fermion superconductors [58,59]. As d x 2 y 2 -wave pairing has four vertical line node on the Fermi surface, vortex core becomes four-fold symmetric shape [26,27]. As the corresponding anisotropic s-wave superconductors, four-fold symmetric vortex core image was observed in YNi 2 B 2 C by STM [10,11,32]. Therefore, detailed studies about the impurity effect of the vortex states are expected also in these superconductors. The twofold or four-fold symmetric vortex core structure may be related to various physical quantities in the anisotropic superconductors under magnetic fields, such as vortex lattice configuration [60] and angle-resolved specific heat measurement [58,59]. Therefore, studies about the impurity effects on these physical quantities are also interesting.
In summary, non-magnetic impurity scattering effect on the vortex core states are studied by the Eilenberger theory in superconductors with p x -wave pairing symmetry, as well as the corresponding anisotropic s-wave symmetry. From the comparison of the pair potential and the local electronic states around a vortex, we examine the differences between anisotropic superconductors with and without sign-change of the pairing function, and find how twofold symmetric vortex core changes with increasing the impurity scattering rate both in the Born and the unitary limits. For example, we found that twofold symmetric vortex core image of zero-energy LDOS changes the orientation of the twofold symmetry with increasing the scattering rate when the sign change occurs in the pairing function. Without the sign change, the vortex core shape reduces to circular one with approaching dirty cases. These results clarify the contributions from the sign-change of the pairing function in anisotropic superconductors, which are helpful to identify the pairing symmetry in unconventional superconductors, in relation to the STM vortex core observation.

Author Contributions

Data curation, Y.S. and T.U.; Investigation, Y.S., T.U., H.A. and M.I.; Visualization, Y.S.; Writing-original draft, M.I.; Writing-review & editing, Y.S. and H.A.; Project administration, M.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by JSPS KAKENHI Grant Number 17K05542 and 19H00657.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abrikosov, A.A. Fundamentals of the Theory of Metals; Elsevier: New York, NY, USA, 1988; Chapter 18. [Google Scholar]
  2. de Gennes, P.G. Superconductivity of Metals and Alloys; Benjamin: New York, NY, USA, 1966; Chapter 3. [Google Scholar]
  3. Fetter, A.L.; Hohenberg, P.C. Theory of Type II Superconductors. In Superconductivity; Parks, R.D., Ed.; Marcel Dekker: New York, NY, USA, 1969; pp. 817–923. [Google Scholar]
  4. Caroli, C.; de Gennes, P.G.; Matricon, J. Bound Fermion states on a vortex line in a type II superconductor. Phys. Lett. 1964, 9, 307–309. [Google Scholar] [CrossRef]
  5. Gygi, F.; Schlüter, M. Electronic tunneling into an isolated vortex in a clean type-II superconductor. Phys. Rev. B 1990, 41, 822–825. [Google Scholar] [CrossRef] [PubMed]
  6. Gygi, F.; Schlüter, M. Self-consistent electronic structure of a vortex line in a type-II superconductor. Phys. Rev. B 1991, 43, 7609–7621. [Google Scholar] [CrossRef] [PubMed]
  7. Hess, H.F.; Robinson, R.B.; Dynes, R.C.; Valles, J.M.; Waszczak, J.V. Scanning-Tunneling-Microscope Observation of the Abrikosov Flux Lattice and the Density of States near and inside a Fluxoid. Phys. Rev. Lett. 1989, 62, 214–216. [Google Scholar] [CrossRef]
  8. Hess, H.F.; Robinson, R.B.; Waszczak, J.V. Vortex-Core Structure Observed with a Scanning Tunneling Microscope. Phys. Rev. Lett. 1990, 64, 2711–2714. [Google Scholar] [CrossRef]
  9. Hess, H.F.; Robinson, R.B.; Waszczak, J.V. Scanning-Tunneling-Microscopy Study of Distortion and Instability of Inclined Flux-Line-Lattice Structures in the Anisotropic Superconductor 2H-NbSe2. Phys. Rev. Lett. 1992, 69, 2138–2141. [Google Scholar] [CrossRef]
  10. Nishimori, H.; Uchiyama, K.; Kaneko, S.; Tokura, A.; Takeya, H.; Hirata, K.; Nishida, N. First Observation of the Fourfold-symmetric and Quantum Regime Vortex Core in YNi2B2C by Scanning Tunneling Microscopy and Spectroscopy. J. Phys. Soc. Jpn. 2004, 73, 3247–3250. [Google Scholar] [CrossRef]
  11. Kaneko, S.; Matsuba, K.; Hafiz, M.; Yamasaki, K.; Kalizaki, E.; Nishida, N.; Takeya, H.; Hirata, K.; Kawakami, T.; Mizushima, T.; et al. Quantum Limiting Behaviors of a Vortex Core in an Anisotropic Gap Superconductor. J. Phys. Soc. Jpn. 2012, 81, 063701. [Google Scholar] [CrossRef] [Green Version]
  12. Eskildsen, M.R.; Kugler, M.; Tanaka, S.; Jun, J.; Kazakov, S.M.; Karpinski, J.; Fischer, Ø. Vortex Imaging in the π Band of Magnesium Diboride. Phys. Rev. Lett. 2002, 89, 187003. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, W.-L.; Zhang, Y.-M.; Lv, Y.-F.; Ding, H.; Wang, L.; Li, W.; He, K.; Song, C.-L.; Ma, X.-C.; Xue, Q.-K. Anisotropic superconductivity and elongated vortices with unusual bound states in quasi-one-dimensional nickel-bismuth compounds. Phys. Rev. B 2018, 97, 134524. [Google Scholar] [CrossRef] [Green Version]
  14. Maggio-Aprile, I.; Renner, C.; Erb, A.; Walker, E.; Fischer, Ø. Direct Vortex Lattice Imaging and Tunneling Spectroscopy of Flux Lines on YBa2Cu3O7−δ. Phys. Rev. Lett. 1995, 75, 2754–2757. [Google Scholar] [CrossRef] [PubMed]
  15. Fischer, Ø.; Kugler, M.; Maggio-Aprile, I.; Berthod, C.; Renner, C. Scanning tunneling spectroscopy of high-temperature superconductors. Rev. Mod. Phys. 2007, 79, 353–419. [Google Scholar] [CrossRef] [Green Version]
  16. Berthod, C.; Maggio-Aprile, I.; Bruér, J.; Erb, A.; Renner, C. Observation of Caroli-de Gennes-Matricon Vortex States in YBa2Cu3O7−δ. Phys. Rev. Lett. 2017, 119, 237001. [Google Scholar] [CrossRef] [Green Version]
  17. Renner, C.; Revaz, B.; Kadowaki, K.; Maggio-Aprile, I.; Fischer, Ø. Observation of the Low Temperature Pseudogap in the Vortex Cores of Bi2Sr2CaCu2O8−δ. Phys. Rev. Lett. 1998, 80, 3606–3609. [Google Scholar] [CrossRef] [Green Version]
  18. Pan, S.H.; Hudson, E.W.; Gupta, A.K.; Ng, K.-W.; Eisaki, H.; Uchida, S.; Davis, J.C. STM Studies of the Electronic Structure of Vortex Cores in Bi2Sr2CaCu2O8−δ. Phys. Rev. Lett. 2000, 85, 1536–1539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Yoshizawa, S.; Koseki, T.; Matsuba, K.; Mochiku, T.; Hirata, K.; Nishida, N. High-resolution scanning tunneling spectroscopy of vortex cores in inhomogeneous electronic states of Bi2Sr2CaCu2Ox. J. Phys. Soc. Jpn. 2013, 82, 083706. [Google Scholar] [CrossRef] [Green Version]
  20. Song, C.-L.; Wang, Y.-L.; Cheng, P.; Jiang, Y.-P.; Li, W.; Zhang, T.; Li, Z.; He, K.; Wang, L.; Jia, J.-F.; et al. Direct Observation of Nodes and Twofold Symmetry in FeSe Superconductor. Science 2011, 332, 1410–1413. [Google Scholar] [CrossRef] [Green Version]
  21. Watashige, T.; Tsutsumi, Y.; Hanaguri, T.; Kohsaka, Y.; Kasahara, S.; Furusaki, A.; Sigrist, M.; Meingast, C.; Wolf, T.; Löhneysen, H.V.; et al. Evidence for Time-Reversal Symmetry Breaking of the Superconducting State near Twin-Boundary Interfaces in FeSe Revealed by Scanning Tunneling Spectroscopy. Phys. Rev. X 2015, 5, 031022. [Google Scholar] [CrossRef] [Green Version]
  22. Hanaguri, T.; Kasahara, S.; Böker, J.; Eremin, I.; Shibauchi, T.; Matsuda, Y. Quantum Vortex Core and Missing Pseudogap in the Multiband BCS-BEC Crossover Superconductor FeSe. Phys. Rev. Lett. 2019, 122, 077001. [Google Scholar] [CrossRef] [Green Version]
  23. Hanaguri, T.; Kitagawa, K.; Matsubayashi, K.; Mazaki, Y.; Uwatoko, Y.; Takagi, H. Scanning tunneling microscopy/spectroscopy of vortices in LiFeAs. Phys. Rev. B 2012, 85, 214505. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, D.; Kong, L.; Fan, P.; Chen, H.; Zhu, S.; Liu, W.; Cao, L.; Sun, Y.; Du, S.; Schneeloch, J.; et al. Evidence for Majorana bound states in an iron-based superconductor. Science 2018, 362, 333–335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Chen, M.; Chen, X.; Yang, H.; Du, Z.; Wen, H.-H. Superconductivity with twofold symmetry in Bi2Te3/FeTe0.55Se0.45 heterostructures. Sci. Adv. 2018, 4, eaat1084. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Schopohl, N.; Maki, K. Quasiparticle spectrum around a vortex line in a d-wave superconductor. Phys. Rev. B 1995, 52, 490–493. [Google Scholar] [CrossRef]
  27. Ichioka, M.; Hayashi, N.; Enomoto, N.; Machida, K. Vortex structure in d-wave superconductors. Phys. Rev. B 1996, 53, 15316–15326. [Google Scholar] [CrossRef] [Green Version]
  28. Ichioka, M.; Hasegawa, A.; Machida, K. Vortex lattice effects on low-energy excitations in d-wave and s-wave superconductors. Phys. Rev. B 1999, 59, 184–187. [Google Scholar] [CrossRef]
  29. Ichioka, M.; Hasegawa, A.; Machida, K. Field dependence of the vortex structure in d-wave and s-wave superconductors. Phys. Rev. B 1999, 59, 8902–8916. [Google Scholar] [CrossRef] [Green Version]
  30. Hayashi, N.; Ichioka, M.; Machida, K. Star-Shaped Local Density of States around Vortices in a Type-II Superconductor. Phys. Rev. Lett. 1996, 77, 4074–4077. [Google Scholar] [CrossRef] [Green Version]
  31. Hayashi, N.; Ichioka, M.; Machida, K. Effects of gap anisotropy upon the electronic structure around a superconducting vortex. Phys. Rev. B 1997, 56, 9052–9063. [Google Scholar] [CrossRef] [Green Version]
  32. Nagai, Y.; Kato, Y.; Hayashi, N.; Yamauchi, K.; Harima, H. Calculated positions of point nodes in the gap structure of the borocarbide superconductor YNi2B2C. Phys. Rev. B 2007, 76, 214514. [Google Scholar] [CrossRef] [Green Version]
  33. Anderson, P.W. Theory of dirty superconductors. J. Phys. Chem. Solids 1959, 11, 26–30. [Google Scholar] [CrossRef]
  34. Abrikosov, A.A.; Gor’kov, L.P. On the Theory of Superconducting Alloys I. The Electrodynamics of Alloys at Absolute Zero. Zh. Eksp. Teor. Fiz. 1958, 35, 1558–1571, [Sov. Phys. JETP 1959, 8, 1090–1098]. [Google Scholar]
  35. Fulde, P. Effects of Gap Anisotropy in Superconductors Containing Paramagnetic Impurities. Phys. Rev. 1965, 139, A726–A728. [Google Scholar] [CrossRef]
  36. Yamashita, T.; Takenaka, T.; Tokiwa, Y.; Wilcox, J.A.; Mizukami, Y.; Terazawa, D.; Kasahara, Y.; Kittaka, S.; Sakakibara, T.; Konczykowski, M.; et al. Fully gapped superconductivity with no sign change in the prototypical heavy-fermion CeCu2Si2. Sci. Adv. 2017, 3, e1601667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Larkin, A.I. Vector Pairing in Superconductors of Small Dimensions. Pis’ma Zh. Éksp. Teor. Fiz. 1965, 2, 205–209, [JETP Lett. 1965, 2, 130–132]. [Google Scholar]
  38. Abrikosov, A.A.; Gor’kov, L.P. Contribution to the theory of superconducting alloys with paramagnetic impurities. Zh. Eksp. Teor. Fiz. 1960, 39, 1781–1796, [Sov. Phys. JETP 1961, 12, 1243–1253]. [Google Scholar]
  39. Balatsky, A.V.; Vekhter, I.; Zhu, J.-X. Impurity-induced states in conventional and unconventional superconductors. Rev. Mod. Phys. 2006, 78, 373–433. [Google Scholar] [CrossRef] [Green Version]
  40. Ptok, A.; Kapcia, K.J. Probe-type of superconductivity by impurity in materials with short coherence length: The s-wave and η-wave phases study. Supercond. Sci. Technol. 2015, 28, 045022. [Google Scholar] [CrossRef] [Green Version]
  41. Eilenberger, G. Transformation of Gorkov’s equation for type II superconductors into transport-like equations. Z. Phys. 1968, 214, 195–213. [Google Scholar] [CrossRef]
  42. Klein, U. Microscopic calculations on the vortex state of type II superconductors. J. Low Temp. Phys. 1987, 69, 1–37. [Google Scholar] [CrossRef]
  43. Miranović, P.; Ichioka, M.; Machida, K. Effects of nonmagnetic scatterers on the local density of states around a vortex in s-wave superconductors. Phys. Rev. B 2004, 70, 104510. [Google Scholar] [CrossRef] [Green Version]
  44. Tanaka, K.K.; Ichioka, M.; Nakai, N.; Machida, K. Knight shift spectrum in vortex states in s- and d-wave superconductors on the basis of Eilenberger theory. Phys. Rev. B 2014, 89, 174504. [Google Scholar] [CrossRef] [Green Version]
  45. Thuneberg, E.V.; Kurkijärvi, J.; Rainer, D. Elementary-flux-pinning potential in type-II superconductors. Phys. Rev. B 1984, 29, 3913–3923. [Google Scholar] [CrossRef]
  46. Kato, Y. Phase-Sensitive Impurity Effects in Vortex Core of Moderately Clean Chiral Superconductors. J. Phys. Soc. Jpn. 2000, 69, 3378–3386. [Google Scholar] [CrossRef] [Green Version]
  47. Eschrig, M.; Rainer, D.; Sauls, J.A. Vortex Core Structure and Dynamics in Layered Superconductors. In Vortices in Unconventional Superconductors and Superfluids; Huebener, R.P., Schopohl, N., Volovik, G.E., Eds.; Springer: Heidelberg, Germany, 2002; pp. 175–293. [Google Scholar]
  48. Hayashi, N.; Kato, Y. Elementary vortex pinning potential in a chiral p-wave superconductor. Phys. Rev. B 2002, 66, 132511. [Google Scholar] [CrossRef] [Green Version]
  49. Hayashi, N.; Kato, Y.; Sigrist, M. Impurity effect on Kramer-Pesch core shrinkage in s-wave vortex and chiral p-wave vortex. J. Low Temp. Phys. 2005, 139, 79–96. [Google Scholar] [CrossRef] [Green Version]
  50. Sauls, J.A.; Eschrig, M. Vortices in chiral, spin-triplet superconductors and superfluids. New J. Phys. 2009, 11, 075008. [Google Scholar] [CrossRef]
  51. Miranović, P.; Machida, K. Thermodynamics and magnetic field profiles in low-κ type-II superconductors. Phys. Rev. B 2003, 67, 092506. [Google Scholar] [CrossRef] [Green Version]
  52. Ichioka, M.; Kogan, V.G.; Schmalian, J. Locking of length scales in two-band superconductors. Phys. Rev. B 2017, 95, 064512. [Google Scholar] [CrossRef] [Green Version]
  53. Hirschfeld, P.J.; Wölfle, P.; Einzel, D. Consequences of resonant impurity scattering in anisotropic superconductors: Thermal and spin relaxation properties. Phys. Rev. B 1988, 37, 83–97. [Google Scholar] [CrossRef]
  54. Kübert, C.; Hirschfeld, P.J. Vortex contribution to specific heat of dirty d-wave superconductors: Breakdown of scaling. Solid State Commun. 1998, 105, 459–463. [Google Scholar] [CrossRef] [Green Version]
  55. Hicks, C.W.; Brodsky, D.O.; Yelland, E.A.; Gibbs, A.S.; Bruin, J.A.N.; Barber, M.E.; Edkins, S.D.; Nishimura, K.; Yonezawa, S.; Maeno, Y.; et al. Strong Increase of Tc of Sr2RuO4 Under Both Tensile and Compressive Strain. Science 2014, 344, 283–285. [Google Scholar] [CrossRef] [PubMed]
  56. Wollman, D.A.; Van Harlingen, D.J.; Lee, W.C.; Ginsberg, D.M.; Leggett, A.J. Experimental Determination of the Superconducting Pairing State in YBCO from the Phase Coherence of YBCO-Pb dc SQUIDs. Phys. Rev. Lett. 1993, 71, 2134–2137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Tsuei, C.C.; Kirtley, J.R.; Chi, C.C.; Yu-Jahnes, L.S.; Gupta, A.; Shaw, T.; Sun, J.Z.; Ketchen, M.B. Pairing Symmetry and Flux Quantization in a Tricrystal Superconducting Ring of YBa2Cu3O7−δ. Phys. Rev. Lett. 1994, 73, 593–596. [Google Scholar] [CrossRef] [PubMed]
  58. An, K.; Sakakibara, T.; Settai, R.; Onuki, Y.; Hiragi, M.; Ichioka, M.; Machida, K. Sign Reversal of Field-Angle Resolved Heat Capacity Oscillations in a Heavy Fermion Superconductor CeCoIn5 and dx2−y2 Pairing Symmetry. Phys. Rev. Lett. 2010, 104, 037002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Sakakibara, T.; Yamada, A.; Custers, J.; Yano, K.; Tayama, T.; Aoki, H.; Machida, K. Nodal Structures of Heavy Fermion Superconductors Probed by the Specific-Heat Measurements in Magnetic Fields. J. Phys. Soc. Jpn. 2007, 76, 051004. [Google Scholar] [CrossRef]
  60. Suzuki, K.M.; Inoue, K.; Miranović, P.; Ichioka, M.; Machida, K. Generic First-Order Orientation Transition of Vortex Lattices in Type II Superconductors. J. Phys. Soc. Jpn. 2010, 79, 013702. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Some pairing functions φ ( θ ) are schematically presented on a circular Fermi surface in ( k x , k y ) plane. (a) φ p x ( θ ) = 2 cos θ for the p x -wave symmetry. (b) φ | p x | ( θ ) = | φ p x ( θ ) | = 2 | cos θ | for anisotropic s-wave symmetry which has the same amplitude as φ p x ( θ ) .
Figure 1. Some pairing functions φ ( θ ) are schematically presented on a circular Fermi surface in ( k x , k y ) plane. (a) φ p x ( θ ) = 2 cos θ for the p x -wave symmetry. (b) φ | p x | ( θ ) = | φ p x ( θ ) | = 2 | cos θ | for anisotropic s-wave symmetry which has the same amplitude as φ p x ( θ ) .
Symmetry 12 00175 g001
Figure 2. (a) Amplitude of the pair potential Δ as a function of 1 / τ 0 in uniform states at a zero field. Δ is plotted in the range 0 1 / τ 0 0.2 for p x -wave pairing φ ( θ ) = φ p x ( θ ) in the Born and the unitary limits, with Δ for anisotropic s-wave paring φ ( θ ) = φ | p x | ( θ ) in the Born limit. T / T c = 0.5 . (b) Δ for anisotropic s-wave paring φ ( θ ) = φ | p x | ( θ ) is plotted in wider range 0 1 / τ 0 25 in the Born limit. Δ in the unitary limit shows almost the same 1 / τ 0 -dependence.
Figure 2. (a) Amplitude of the pair potential Δ as a function of 1 / τ 0 in uniform states at a zero field. Δ is plotted in the range 0 1 / τ 0 0.2 for p x -wave pairing φ ( θ ) = φ p x ( θ ) in the Born and the unitary limits, with Δ for anisotropic s-wave paring φ ( θ ) = φ | p x | ( θ ) in the Born limit. T / T c = 0.5 . (b) Δ for anisotropic s-wave paring φ ( θ ) = φ | p x | ( θ ) is plotted in wider range 0 1 / τ 0 25 in the Born limit. Δ in the unitary limit shows almost the same 1 / τ 0 -dependence.
Symmetry 12 00175 g002
Figure 3. DOS spectrum N ( E ) in uniform states at a zero field. (a) N ( E ) for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . 1 / τ 0 = 0.01 , 0.05, 0.1, and 0.15 in the Born limit. (b) The same as (a), but in the unitary limit. (c) N ( E ) for anisotropic s-wave paring φ ( θ ) = φ | p x | ( θ ) . 1 / τ 0 = 0.01 , 0.1, 0.6, and 1.0 in the Born limit. (d) The same as (c), but in the unitary limit.
Figure 3. DOS spectrum N ( E ) in uniform states at a zero field. (a) N ( E ) for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . 1 / τ 0 = 0.01 , 0.05, 0.1, and 0.15 in the Born limit. (b) The same as (a), but in the unitary limit. (c) N ( E ) for anisotropic s-wave paring φ ( θ ) = φ | p x | ( θ ) . 1 / τ 0 = 0.01 , 0.1, 0.6, and 1.0 in the Born limit. (d) The same as (c), but in the unitary limit.
Symmetry 12 00175 g003
Figure 4. (a) Spatial structure of the pair potential | Δ ( r ) | around a vortex core within a unit cell of the vortex lattice in x y plane for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . 1 / τ 0 = 0.01 , 0.05, 0.10, and 0.15. The upper (middle) panels are for the Born (unitary) limit. Lower panels show profile of | Δ ( r ) | along the x (red lines) and the y axes (blue lines). The solid (dashed) lines are for the Born (unitary) limit. (b) The same as panel (a), but for anisotropic s wave pairing symmetry φ ( θ ) = φ | p x | ( θ ) . 1 / τ 0 = 0.01 , 0.1, 0.6, and 1.0.
Figure 4. (a) Spatial structure of the pair potential | Δ ( r ) | around a vortex core within a unit cell of the vortex lattice in x y plane for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . 1 / τ 0 = 0.01 , 0.05, 0.10, and 0.15. The upper (middle) panels are for the Born (unitary) limit. Lower panels show profile of | Δ ( r ) | along the x (red lines) and the y axes (blue lines). The solid (dashed) lines are for the Born (unitary) limit. (b) The same as panel (a), but for anisotropic s wave pairing symmetry φ ( θ ) = φ | p x | ( θ ) . 1 / τ 0 = 0.01 , 0.1, 0.6, and 1.0.
Symmetry 12 00175 g004
Figure 5. (a) Impurity scattering rate 1 / τ 0 dependence of vortex core radius ξ x and ξ y , which are determined from the slope of | Δ ( r ) | at the vortex center, for φ p x ( θ ) in the range 0 1 / τ 0 0.15 and φ | p x | ( θ ) in the range 0 1 / τ 0 1 . The Born (solid lines) and the unitary (dashed lines) limit shows almost identical 1 / τ 0 dependence. (b) The anisotropy ratio ξ y / ξ x as a function of 1 / τ 0 , from ξ x and ξ y in panel (a).
Figure 5. (a) Impurity scattering rate 1 / τ 0 dependence of vortex core radius ξ x and ξ y , which are determined from the slope of | Δ ( r ) | at the vortex center, for φ p x ( θ ) in the range 0 1 / τ 0 0.15 and φ | p x | ( θ ) in the range 0 1 / τ 0 1 . The Born (solid lines) and the unitary (dashed lines) limit shows almost identical 1 / τ 0 dependence. (b) The anisotropy ratio ξ y / ξ x as a function of 1 / τ 0 , from ξ x and ξ y in panel (a).
Symmetry 12 00175 g005
Figure 6. (a) Spatial structure of zero-energy LDOS N ( E = 0 , r ) around a vortex core in x y plane for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . 1 / τ 0 = 0.01 , 0.05, 0.10, and 0.15. The upper (middle) panels are for the Born (unitary) limit. Lower panels show profile of N ( E = 0 , r ) along the x (red lines) and the y axes (blue lines). The solid (dashed) lines are for the Born (unitary) limit. (b) The same as panel (a), but for anisotropic s wave pairing symmetry φ ( θ ) = φ | p x | ( θ ) . 1 / τ 0 = 0.01 , 0.1, 0.6, and 1.0.
Figure 6. (a) Spatial structure of zero-energy LDOS N ( E = 0 , r ) around a vortex core in x y plane for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . 1 / τ 0 = 0.01 , 0.05, 0.10, and 0.15. The upper (middle) panels are for the Born (unitary) limit. Lower panels show profile of N ( E = 0 , r ) along the x (red lines) and the y axes (blue lines). The solid (dashed) lines are for the Born (unitary) limit. (b) The same as panel (a), but for anisotropic s wave pairing symmetry φ ( θ ) = φ | p x | ( θ ) . 1 / τ 0 = 0.01 , 0.1, 0.6, and 1.0.
Symmetry 12 00175 g006
Figure 7. (a) Peak height N ( E = 0 , r = 0 ) at the vortex center as a function of impurity scattering rate 1 / τ 0 for φ p x ( θ ) in the range 0 1 / τ 0 0.15 and φ | p x | ( θ ) in the range 0 1 / τ 0 1 . The solid (dashed) lines are for the Born (unitary) limit. (b) 1 / τ 0 dependence of vortex core radius ξ x and ξ y , which are determined from the zero energy LDOS N ( E = 0 , r ) around a vortex core as described in text, for φ p x ( θ ) and φ | p x | ( θ ) . Solid (dashed) lines are for the Born (unitary) limit. (c) The anisotropy ratio ξ y / ξ x as a function of 1 / τ 0 , from ξ x and ξ y in panel (b).
Figure 7. (a) Peak height N ( E = 0 , r = 0 ) at the vortex center as a function of impurity scattering rate 1 / τ 0 for φ p x ( θ ) in the range 0 1 / τ 0 0.15 and φ | p x | ( θ ) in the range 0 1 / τ 0 1 . The solid (dashed) lines are for the Born (unitary) limit. (b) 1 / τ 0 dependence of vortex core radius ξ x and ξ y , which are determined from the zero energy LDOS N ( E = 0 , r ) around a vortex core as described in text, for φ p x ( θ ) and φ | p x | ( θ ) . Solid (dashed) lines are for the Born (unitary) limit. (c) The anisotropy ratio ξ y / ξ x as a function of 1 / τ 0 , from ξ x and ξ y in panel (b).
Symmetry 12 00175 g007
Figure 8. (a) Spatial structure of LDOS N ( E , r ) around a vortex core in x y plane for the p x -wave pairing φ ( θ ) = φ p x ( θ ) at E = 0.01 , 0.1, 0.2, and 0.5. 1 / τ 0 = 0.01 (upper panels), 0.10 (middle panels) in the Born limit, and 1 / τ 0 = 0.1 in the unitary limit (lower panels). (b) Profile of N ( E , r ) along the x (red lines) and the y axes (blue lines) for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . The solid (dashed) lines are for the Born (unitary) limit. 1 / τ 0 = 0.01 (upper panels) and 0.1 (lower panels).
Figure 8. (a) Spatial structure of LDOS N ( E , r ) around a vortex core in x y plane for the p x -wave pairing φ ( θ ) = φ p x ( θ ) at E = 0.01 , 0.1, 0.2, and 0.5. 1 / τ 0 = 0.01 (upper panels), 0.10 (middle panels) in the Born limit, and 1 / τ 0 = 0.1 in the unitary limit (lower panels). (b) Profile of N ( E , r ) along the x (red lines) and the y axes (blue lines) for the p x -wave pairing φ ( θ ) = φ p x ( θ ) . The solid (dashed) lines are for the Born (unitary) limit. 1 / τ 0 = 0.01 (upper panels) and 0.1 (lower panels).
Symmetry 12 00175 g008
Figure 9. The same as Figure 8, but for anisotropic s wave pairing symmetry φ | p x | ( θ ) .
Figure 9. The same as Figure 8, but for anisotropic s wave pairing symmetry φ | p x | ( θ ) .
Symmetry 12 00175 g009
Figure 10. (a) Bold line schematically shows main distribution of the k -resolved LDOS N ( E , k , r ) . Vortex center is located at O(0,0). See text for details. (b) Quasiparticle trajectories at low E for the p x -wave pairing symmetry φ p x ( θ ) are schematically presented. (c) The same as panel (b), but at higher E.
Figure 10. (a) Bold line schematically shows main distribution of the k -resolved LDOS N ( E , k , r ) . Vortex center is located at O(0,0). See text for details. (b) Quasiparticle trajectories at low E for the p x -wave pairing symmetry φ p x ( θ ) are schematically presented. (c) The same as panel (b), but at higher E.
Symmetry 12 00175 g010

Share and Cite

MDPI and ACS Style

Sera, Y.; Ueda, T.; Adachi, H.; Ichioka, M. Relation of Superconducting Pairing Symmetry and Non-Magnetic Impurity Effects in Vortex States. Symmetry 2020, 12, 175. https://doi.org/10.3390/sym12010175

AMA Style

Sera Y, Ueda T, Adachi H, Ichioka M. Relation of Superconducting Pairing Symmetry and Non-Magnetic Impurity Effects in Vortex States. Symmetry. 2020; 12(1):175. https://doi.org/10.3390/sym12010175

Chicago/Turabian Style

Sera, Yasuaki, Takahiro Ueda, Hiroto Adachi, and Masanori Ichioka. 2020. "Relation of Superconducting Pairing Symmetry and Non-Magnetic Impurity Effects in Vortex States" Symmetry 12, no. 1: 175. https://doi.org/10.3390/sym12010175

APA Style

Sera, Y., Ueda, T., Adachi, H., & Ichioka, M. (2020). Relation of Superconducting Pairing Symmetry and Non-Magnetic Impurity Effects in Vortex States. Symmetry, 12(1), 175. https://doi.org/10.3390/sym12010175

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop