Analytical Charge Density Profile of Vortex Core in Weak-Coupling Superconductor

기술 요약: 약결합 초전도체 내 보텍스 코어의 분석적 전하 밀도 프로파일

문제 정의
제2종 초전도체의 보텍스 코어(vortex core) 내 전자 구조가 이산적인 준입자 결합 상태(Caroli–de Gennes–Matricon 또는 CdGM 상태)를 수용한다는 사실은 오래전부터 알려져 왔다. 자기 일관적(self-consistent) Bogoliubov–de Gennes (BdG) 계산과 포아송 방정식을 결합한 수치 해석을 통해, 보텍스 코어 내부의 전하 밀도가 $\pi/k_F$ 주기로 부호가 교대하는 진동을 보인다는 점이 수치적으로 입증되었으나, 이러한 현상에 대한 단순하고 분석적인 설명은 오랫동안 부족했다. 이전의 수치 연구들(예: Machida와 Koyama)은 스크리닝된 전하 밀도가 부호를 바꾼다는 것을 관찰했으나, 이 부호 반전의 메커니즘—구체적으로 결합 상태의 전하 분포와 정전기적 스크리닝 사이의 상호작용—이 제1원리로부터 폐쇄형(closed form)으로 유도되지는 않았다. 난제는 BdG 방정식의 복잡성에 있는데, 여기서 규격화 상수(normalization constant)는 일반적으로 모드에 따라 달라지며 파동함수는 결합된 비선형 미분 방정식들을 포함하므로, 분석적인 합산(summation)이 매우 어렵다.

방법론
저자는 보텍스 코어의 핵심 물리량을 유지하면서 가장 경제적인 모델을 사용하는 기초적이고 설명적인 접근법을 채택한다:

1. 모델 시스템: 계단형 갭 프로파일($r < \xi$일 때 $\Delta(r) = 0$, 그 외에는 $\Delta(r) = \Delta_\infty$)을 갖는 약결합 $s$-파 초전도체 내의 단일 보텍스 선. 스펙트럼은 표준적인 CdGM 결과인 $E_\mu = \mu \Delta^2/E_F$를 사용한다.
2. 파동함수 구성:
   • 내부 영역 ($r < \xi$): 갭이 사라지므로 BdG 방정식이 디커플링(decoupling)된다. 저자는 WKB/BKJT 안사츠(ansatz) 대신 정확한 베셀 함수($J_{\mu \mp 1/2}$) 해를 사용한다. 이는 위상적 회전(topological winding)으로 인해 보텍스 선 근처에서 WKB/BKJT 방식이 정의되지 않기 때문이다.
   • 외부 영역 ($r > \xi$): BKJT 엔벨로프(envelope)를 선형화하며, 계수는 코어 경계에서의 매칭 조건에 의해 결정된다.
3. 규격화 분석: 핵심 단계는 규격화 상수 $C_\mu$를 계산하는 것이다. 저자는 고전적으로 허용되는 창(window)과 감쇠율이 변함에도 불구하고, 대부분의 CdGM 래더(ladder)에서 규격화 상수가 각 각운동량 양자수 $\mu$에 대해 거의 독립적임을 입증한다. 이를 통해 $C_\mu$를 모드 합산 시 상수 전계 인자 $\langle C^2 \rangle$로 취급할 수 있다.
4. 모드 합산:
   • 코어 내부: 전하 밀도 합은 부분 베셀-파르스발 합(Bessel–Parseval sum)으로 축소된다. 베셀 함수의 완비성 항등식(completeness identity)을 적용하고, 위상적 제약에 의해 $\nu=0$ 채널이 제외됨을 이용하여, 이 합은 폐쇄형으로 수렴한다.
   • 코어 외부: 합은 잔여 리플(ripple)을 포함하는 감쇠하는 $1/r$ 엔벨로프를 생성한다.
5. 정전기학: 결합 상태의 전하 밀도는 포아송 방정식의 소스 항(source term) 역할을 한다. 저자는 운동량 공간에서 스크리닝된 포아송 방정식을 풀며, 매끄러운 배경 성분과 $2k_F$ 진동 성분을 분리하여 취급한다.

주요 기여 및 결과

• 폐쇄형 내부 밀도: 저자는 코어 내부의 결합 상태 전하 밀도에 대한 분석적 표현식을 유도한다: $\delta n(r < \xi) \propto 1 - J_0^2(k_F r)$. 이 결과는 경계 조건 $n(0)=0$을 만족하며($\nu=0$ 베셀 채널의 제외로 인해), 밀도가 매끄러운 배경과 $2k_F$ 파수를 갖는 프라이델(Friedel) 유사 진동으로 구성됨을 보여준다.
• 스크리닝 메커니즘: 저자는 왜 스크리닝된 전하 밀도가 부호를 교대하는지 밝혀낸다. 결합 상태의 전하는 적분값이 0이 아니므로, 전체적인 전하 중성을 달성하기 위해서는 확장된 상태(extended states)가 이를 보상해야 한다.
  • 토마스-페르미(Thomas–Fermi) 스크리닝은 장파장(낮은 운동량 전달)에서 매우 효율적이어서, 전하의 매끄러운 배경 성분을 완전히 중화시킨다.
  • 그러나 스크리닝은 큰 운동량 전달($q \approx 2k_F$)에서는 비효적으로 작용한다. 결과적으로, 결합 상태 전하의 $2k_F$ 진동 성분은 스크리닝 과정을 거친 후에도 살아남는다.
• 생존하는 전하에 대한 분석적 표현식: 최종적인 스크리닝된 전하 밀도는 다음과 같다:
  $$\delta n_{\text{total}}(r) \simeq -\langle C^2 \rangle \frac{4k_F^2}{4k_F^2 + k_{TF}^2} \frac{\sin(2k_F r)}{\pi k_F r}$$
  이 식은 생존하는 진동이 부호가 교대하는 $2k_F$ 변조임을 보여준다. 진폭은 $4k_F^2/(4k_F^2 + k_{TF}^2)$라는 인자에 의해 감소하며, 이 값은 일반적인 금속의 경우 $1/2$에서 $3/4$ 사이에 위치한다.
• 검증: 분석적 결과는 전체 CdGM 모드 합을 사용한 직접적인 포아송 방정식의 수치 해와 대조하여 검증되었다. 분석적 폐쇄형은 수치 데이터와 1% 이내로 일치하였다(잔여 오차는 규격화 상수의 모드 평균 근사에 기인한다).

의의 및 주장
본 논문은 수치 시뮬레이션에서 관찰된 부호 변화 전하 밀도 프로파일에 대한 기초적인 설명을 제공한다고 주장한다. 주요 의의는 다음과 같다:

1. 메커니즘의 명확성: 저자는 이 진동이 미세한 효과가 아니라, 전하 중성과 큰 운동량 전달($2k_F$)에서의 토마스-페르미 스크리닝의 비효율성의 직접적인 결과임을 확인하였다.
2. 분석적 벤치마크: 저자는 갭 프로파일이 경직된 약결합 한계에서 향후 자기 일관적 수치 계산을 검증할 수 있는 폐쇄형 분석적 벤치마크(식 18, 22, 28)를 제공한다.
3. 불순물 산란과의 구별: 저자는 주기성($2k_F$)은 불순물 주변의 프라이델 진동과 동일하지만, 메커니즘은 반대임을 명확히 한다. 불순물 산란에서는 소스가 구조가 없으며 진동은 응답 함수의 비해석성(non-analyticity)에서 발생하는 반면, 보텍스 코어에서는 소스 자체가 $2k_F$에서 피크를 가지며, 응답(스크리닝)은 단지 매끄러운 성분을 필터링할 뿐이다.

저자는 자신의 목적이 예측보다는 설명에 있으며, 강결합 영역이나 복잡한 갭 대칭성에 대한 정량적 세부 사항을 해결하려는 것이 아님을 명시적으로 밝히면서도, 이를 자연스러운 확장 과제로 언급하였다. 본 연구는 약결합 초전도체의 보텍스 코어에서 전하 밀도의 "리플(ripples)"이 견고한 특징임을 확인시켜 준다.